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EXECUTIVE  SUMMARY 


In  this  work,  the  effect  of  processing  on  the  origin  and  the  evolution  of  interfaces  in 
metallic/intermetallic  matrix  composites  were  studied.  Novel  computa,tiond  methods  using 
nonlinear  finite  element  method  combined  with  experimental  work 

push-out  testing)  were  used  in  this  investigation.  Typical  Titanium  based  MMCs  and  IMCs 
Lh  as  SiC/Ti-6Al-4V,  SiC/Ti-15V-3Sn-3Al  and  SiC/Ti-24Al-llNb  systems  were  used  as 
model  material  in  this  study.  MMCs  and  IMCs  are  significantly  infiuenced  by  the  presence 
of  large  levels  of  residual  stresses  arising  from  processing  as  well  as  use  temperature  and 
time  This  is  quite  different  from  that  of  other  composite  systems  such  as  polymeric  and 
ceramic  matrix  composite  systems.  Thus  great  caution  needs  to  be  exercised  when  results 
are  extended  from  these  systems  to  MMCs.  A  comprehensive  study  of  the  origin  and  role  of 
residual  stresses  on  the  interface  behavior  was  carried  out  in  this  work. 

Interfaces  in  composites  play  a  very  vital  role  in  determining  the  strength,  stiffoess 
and  failure  characteristics  of  composite  systems.  Hence  in  this  work,  the  thermoinechan- 
ical  behavior  of  interfaces  was  carefully  analyzed.  The  strength  and  toughness  behaviour 
of  composites  were  modeled  using  energy  based  criterion.  The  shear  strength  and  fixture 
toughness  of  interfaces  were  evaluated  from  push-out  test  results  in  conjunction  with  non¬ 
linear  finite  element  analysis  incorporating  equivalent  domain  integral  methods.  This  report 
contains  details  of  the  results  in  the  areas  of  elevated  temperature  effects  evolution  of  inter¬ 
face  chemistry,  effect  of  residual  stress  on  fracture  behavior,  the  effect  of  debonding  ^d  fiber 
fracture  on  the  evolution  of  damage  in  composites  and  the  prediction  of  composite  behavior 

given  the  interfacial  properties. 
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Chapter  1 

Overview  of  the  Research  Effort 


This  research  was  focused  on  the  effect  of  processing  on  the  evolution  of  interfaces  in  metal 
matrix  composites  (MMCs).  The  effort  was  carried  out  in  two  different  heirarchical  levels 
(micromechancical  and  macroscopic)  using  both  experimental  and  computational  methods. 
Significant  work  was  carried  out  in  understanding  the  origin  of  interfacial  mechanical  prop¬ 
erties  and  hence  was  performed  at  the  micromechanics  level.  The  report  consists  of  about 
ten  chapters  spanning  different  aspects  of  the  problem.  While  the  first  eight  attempts  to 
characterize  the  interface  in  terms  of  its  thermomechanical  properties,  the  last  two  uses  the 
concept  developed  in  the  earlier  chapters  to  design  composites  for  practical  applications. 


1.1  Introduction 

Metallic  and  Intermrtallic  Matrix  Composites  (MMCs  and  IMCs)  are  being  considered  for 
crucial  applications  in  the  aerospace  industry  because  of  their  high  specific  stiffness  and 
strength  even  at  elevated  temperatures  (>  600°C).  A  critical  issue  in  the  success^  applica¬ 
tion  of  these  composites  is  the  behavior  of  fiber-matrix  interfaces.  Push-out  test  is  emerging 
as  an  important  experimental  tool  for  characterizing  the  interfacial  behavior  of  this  class  of 
composites,  mainly  after  its  widespread  use  in  Ceramic  Matrix  Composites  (CMCs). 

When  MMCs  and  IMCs  are  cooled  after  processing  from  a  stress  free  temperature,  residual 
stresses  arise  due  to  the  difference  in  the  coefficients  of  thermal  expansion  (CTE)  of  fiber 
and  matrix  materials.  These  stresses  reach  considerable  levels  because,  in  typical  titanium 
based  MMCs  and  IMCs  such  as  SiC/Ti-6Al-4V,  SiC/Ti-15V-3Sn-3Al,  SiC/Ti-24Al-llNb 
the  mismatch  in  the  CTE  is  quite  significant  {amatrix  >  ^Oifiber)  and  also  the  temperature 
differential  involved  in  the  cooling  process  is  high  (about  800  °C).  These  stresses  cause 
matrix  cracking  near  the  interface  on  processing  and  also  affect  the  interfacial  response  of 
these  composites  under  service  conditions. 
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1.2  Elevated  Temperature  Interfacial  Behavior  of  IVIlVICs 

The  thermomechanical  behavior  of  fiber-matrix  interfaces  play  an  iihportant  role  in  the  suc¬ 
cessful  application  of  this  class  of  composites.  Push-out  test  is  emerging  as  an  important 
experimental  tool  for  characterizing  the  interfacial  behavior  of  MMCs  and  IMCs.  In  this 
study,  single  fiber  push-out  test  is  modeled  using  the  finite  element  method,  with  the  objec¬ 
tive  of  studying  the  interface  failure  process,  and  extracting  interfacial  properties  from  the 
experimental  test  results.  Earlier  studies  by  the  author’s  group  emphasize  the  significance 
of  processing  induced  residual  stresses  in  titanium  based  composites  and  their  effects  on 
push-out  test  results.  In  the  present  work,  the  developed  methodology  is  used  to  study  the 
interfacial  behavior  during  push-out  tests  at  elevated  temperatures.  An  attempt  is  made  in 
Chapter  2  to  predict  interfacial  shear  strengths  at  elevated  temperatures,  by  correlating  the 
numerical  simulations  and  the  experimental  results. 

In  this  chapter,  it  is  shown  that  residual  shear  stresses  aflFect  the  peak  push-out  load 
as  well  as  the  failure  initiation  load  at  elevated  temperatures.  The  average  shear  strength 

_  P^gj.j'KdV)  evaluated  from  the  peak  push-out  load  does  not  reflect  the  actual  shear 
strSigth  of  the  interface,  since  it  is  influenced  by  the  presence  of  processing  induced  residual 
stresses  and  other  geometric  and  testing  variables.  Numerical  predictions  show  that  for 
SiC/Ti-15-3  the  intrinsic  interfacial  shear  strength  (which  is  a  material  property)  is  about 
four  times  at  room  temperature  and  drops  to  2.0  at  700°C.  The  support  distance  in  the 
experimental  set-up  has  to  be  controlled  properly  at  elevated  temperatures  to  minimize  the. 
specimen  bending  effects. 


1.3  Effect  of  Interface  Chemistry  on  the  Fracture  Prop¬ 
erties 


In  Chapter  3,  the  effect  of  fiber-matrix  interfacial  reactions  on  the  fracture  properties  of 
the  interface  is  studied  using  experimental  characterization  and  computational  modeling 
techniques.  The  objective  of  this  study  is  to  establish  a  link  between  the  evolution  of  the 
interfacial  chemistry  and  the  resulting  mechanical  properties.  SCS-6/Timetal21s  composite 
is  chosen  as  the  candidate  material  system.  The  composite  specimens  are  exposed  to  tem¬ 
peratures  as  high  as  927'’C  for  extended  hours.  The  diffusion  of  elements  across  the  interface 
is  investigated  through  metallurgical  techniques.  Fiber  push-out  is  used  to  characterize  the 
mechanical  properties  of  the  interface.  A  novel  computational  method  is  used  to  simulate 
the  propagation  of  interfacial  cracks  during  the  tests.  The  fracture  toughness  of  the  interface 
is  evaluated  from  the  experimental  data  using  this  model. 

The  evolution  of  interface  in  SCS-6/Timetal  21s  composite  system  is  shown  in  Chapter 
3  by  means  of  heat  treatment  experiments,  with  the  objective  of  establishing  a  link  between 
the  evolution  of  chemistry  and  geometry  of  the  reaction  zone  and  the  interfacial  properties. 
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•  Microstructural  studies  reveal  that  reaction  growth  is  due  to  transformation  of  both 
coating  and  the  matrbc.  The  growth  of  the  reaction  zone  into  the  matrbc  is  much  more 
than  the  growth  towards  the  fiber,  as  carbon  is  able  to  diffuse  much  farther  into  the 
matrix  than  the  diffusion  of  Titanium  into  the  fiber. 

•  Reaction  zone  sizes  increase  rapidly  with  time  at  927°  C.  At  lower  temperatures  such 
as  700°  C  and  450°  C  no  significant  increase  in  reaction  zone  size  is  observed  for  this 

system. 

•  Push-out  tests  conducted  on  heat  treated  specimens  reveal  that  the  effect  of  thermal 
exposure  on  interfacial  properties  is  predominantly  due  to  exposure  temperature.  Time 
of  exposure  is  found  to  have  only  a  secondary  effect. 

•  Fracture  toughness  values  remain  more  or  less  constaiit  for  normal  use  temperatures 
in  the  range  of  450°  to  700°  C.  The  values  however,  increase  when  the  exposure  tem¬ 
perature  is  raised  to  927°  C. 

•  Fracture  surfaces  of  the  specimens  exposed  at  927°  C  show  that  debonding  took  place 
in  the  reaction  zone  whereas  in  as-processed  specimens  and  those  heat  treated  at  the 
lower  temperature  ranges  debonding  initiated  in  the  coating.  This  is  because  the  weak 
region  in  the  coating  is  replaced  by  the  stronger  interfacial  region. 


1.4  Effect  of  residual  stress  on  interfacial  fracture  be¬ 
havior 


A  fracture  mechanics  approach  is  used  in  Chapter  4  to  examine  the  interfacial  debonding 
process  in  MMCs  and  CMCs  during  a  fiber  push-out  test.  The  Equivalent  Domain  Integral 
(EDI)  method  is  implemented  in  a  finite  element  code  and  is  used  to  compute  the  strain 
energy  release  rates  for  the  interface  crack.  The  cooling  process  from  the  composite  con¬ 
solidation  temperature,  specimen  preparation  for  the  push-out  test  and  the  actual  testing 
are  included  in  the  finite  element  simulation.  A  strain  energy  based  debonding  criterion  is 
used  to  predict  the  interfacial  behavior.  The  experimentally  observed  phenomenon  of  bot¬ 
tom  debonding  in  MMCs  is  explained  from  the  energy  release  rate  variation  for  the  loading 
and  support  end  cracks.  It  is  shown  that  processing  induced  residual  stresses  significantly 
affect  the  initiation  and  propagation  of  interface  cracks.  The  advantage  of  the  EDI  method 
over  conventional  methods  for  modeling  interface  crack  propagation,  by  eliminating  the  need 
for  singular  elements  and  thus  remeshing  with  crack  advance  is  demonstrated  through  the 
simulation  of  the  push-out  test. 

Strain  energy  release  rates  for  the  fiber-matrix  interface  due  to  processing  induced 
residual  stresses  are  determined  by  the  above  method.  These  rates  serve  as  a  lower  bound 
for  the  Mode  II  fracture  toughness  of  the  composite  systems  examined.  The  phenomenon 
of  bottom  (support  end)  debonding  observed  in  MMC  systems  during  the  push-out  test  is 
explained  by  the  study  of  the  variation  of  energy  release  rates  for  the  loading  and  support  end 
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cracks  upto  debond  initiation.  The  effect  of  residual  stresses  on  crack  initiation  is  studied  by 
altering  the  residual  stress  state  by  plastic  straining.  It  is  found  that  crack  initiation  switches 
from  the  support  to  the  loading  end  with  reduction  in  residual  stresses  with  increase  in  plastic 
strain.  For  the  specimen  geometry  and  material  system  studied,  this  switching  occurs  when 
the  Gir  values  drop  to  values  of  the  order  of  AJ/w?'  from  the  as- received  value  of  around 
AQJ/rr?.  It  is  shown  that  the  EDI  method  can  be  used  to  simulate  interface  crack  propagation 
without  the  need  for  remeshing  with  crack  advance.  Evaluation  of  strain  energy  release  rates 
through  the  use  of  J-integral  for  cracks  along  bimaterial  interfaces  in  the  presence  of  friction 
is  still  a  “mathematical  problem”  and  care  needs  to  be  exercised  in  the  interpretation  of  the 
numerical  predictions. 


1,5  Evaluation  of  fracture  toughness  using  thin-slice 
push-out  Tests 


Push-out  test  has  emerged  as  an  important  experimental  tool  for  characterizing  the  interface 
behavior  of  MMCs  mainly  after  its  widespread  use  in  the  study  of  interfaces  in  ceramic 
matrix  composites  (CMC).  The  presence  of  high  levels  of  residual  stresses  in  MMCs  alters 
the  mechanics  of  the  failure  process,  making  the  interpretation  of  test  results  and  evaluation 
of  interfacial  properties  very  difficult.  Moreover  the  relatively  high  bond  strength  of  MMCs 
and  residual  clamping  stresses  at  the  interface  limit  the  specimen  thickness  for  push-out 
experiments.  In  the  thin  slice  push-out  tests  it  is  observed  that  interface  failure  initiates 
from  the  bottom  of  the  specimen  and  not  from  the  loading  end  as  in  the  case  of  CMCs.  In 
Chapter  5,  fracture  toughness  (Gic)  is  considered  as  a  more  appropriate  quantitative  measure 
of  the  bond  at  the  interface.  The  interfacial  fracture  energies  of  metal  matrix  composites 
are  not  currently  available  in  the  literature  because  of  the  difficulty  in  evaluating  these  from 
the  experimental  results. 

In  'Lis  chapter,  a  numerical  method  based  on  the  equivalent  domain  integral  technique 
developed  previously  is  used  to  simulate  the  fracture  process  during  a  thin  slice  push-out  test 
and  determine  the  critical  interfacial  fracture  toughness  Gic-  Experimental  correlations  are 
made  with  push-out  test  results  of  ceramic  matrix  composites  to  validate  the  Gic  predictions 
using  the  EDI  method.  This  approach  is  then  used  to  evaluate  the  Gic  of  some  as  processed 
SCS-6/Titanium  matrix  composites  from  the  experimental  data. 


1.6  Effect  of  Fiber  Fracture  and  Matrix  Yielding  on 
Load  Sharing 


Chapter  6  addresses  the  issue  of  stress  re-distribution  in  the  presence  of  a  fiber  fracture  in 
a  composite  lamina.  Stress  elevation  in  the  fiber  adjacent  to  the  broken  fiber  is  the  focus  of 
this  study.  The  stress  concentration  effects  in  the  vicinity  of  the  fiber  break  and  its  influence 
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on  the  neighboring  intact  fiber  is  analyzed  using  finite  element  method  as  a  function  of 
fiber  volume  fraction.  The  role  of  the  inelastic  behavior  of  the  matrix  in  causing  the  stress 
elevation  is  studied.  It  is  found  that  the  state  of  stress  in  the  intact  fiber  is  significantly 
affected  by  the  propagation  of  plastic  front  due  to  matrix  yielding.  The  local  stress  field 
is  affected  also  if  the  fiber  break  is  followed  by  interfacial  debonding.  Titanium  matrix 
reinforced  with  continuous  fibers  of  silicon  carbide  is  chosen  as  the  metal  matrix  composite 
system  for  this  study.  Experimental  comparisons  are  made  with  tencile  tests  conducted  using 
a  single  ply  lamina. 


1.7  Effect  of  Fiber  Fracture  and  Matrix  Yielding  on 
Load  Sharing 

A  new  approach  for  modeling  the  behavior  laminated  composite  structures  using  computa¬ 
tional  methods  is  presented  in  Chapter  7  considering  the  damage  evolution  at  the  microme¬ 
chanical  level.  Micromechanical  models  are  developed  to  predict  the  stress-strain  response 
of  a  composite  lamina  explicitly  accounting  for  the  local  damage  mechanisms  such  as  fiber 
fracture  and  interfacial  bonding.  The  model  is  applied  to  metal  matrix  composites  and  hence 
the  inelastic  constitutive  behavior  of  the  matrix  phase  is  included.  The  stochastic  variation 
of  the  fiber  properties  is  incorporated  in  this  simulation  using  the  two-parameter  weibull 
model.  The  effect  of  fiber  volume  fraction  and  the  properties  of  the  fiber,  matrbc  and  inter¬ 
face  on  the  damage  evolution  is  studied  using  this  approach.  A  constitutive  damage  tensor 
for  the  composite  lamina  is  developed  from  the  micromechanical  models  which  can  be  input 
into  laminate  structural  analysis  codes. 


1.8  Effect  of  Interfacial  Properties  on  the  Fiber  Frag¬ 
mentation  Process 


Fiber  FVagmentation  Test  or  Single  Fiber  Composite  Test  (SFC)  has  been  widely  used  to 
characterize  the  interfacial  behavior  in  composites.  Though  it  is  accepted  that  SFC  gives 
useful  information  about  the  interfacial  bond  quality  and  the  fiber  strength,  there  is  no 
clear  consensus  on  how  to  interpret  the  data  and  quantify  the  interfacial  properties.  In 
this  paper,  the  fiber  fragmentation  test  is  modeled  with  the  objective  of  studying  the  effect 
of  interfacial  bond  strength  and  post-debonding  frictional  effects  on  the  progress  of  the 
fragmentation  process.  The  load  transfer  at  the  fiber-matrix  interface  is  modeled  using 
the  shear-lag  principles  and  the  statistical  variation  of  fiber  strength  is  incorporated*  using 
Weibull  theory.  A  parametric  study  of  the  variation  of  interfacial  shear  properties  on  the 
mean  fragment  length  and  debond  progression  with  applied  strain  is  conducted.  It  is  found 
that  in  the  case  of  weak  interfaces  saturation  of  fiber  fragments  occurs  at  much  lower  applied 
strain  values  and  is  usually  associated  with  extensive  debonding.  The  numerical  simulation 
shows  that  the  fragmentation  results  are  very  much  sensitive  to  the  frictional  effects  as  it 
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controls  the  debonded  length.  The  results  from  the  model  axe  correlated  with  experimental 
data  to  validate  the  model  and  for  extracting  the  interfacial  properties  from  the  test  data. 


1.9  Summary  and  Conclusion 


The  final  chapter  summarizes  the  goal  of  the  research,  technical  approach  taken,  and  some 
of  the  significant  findings.  Future  research  direction  to  pursue  this  work  is  also  suggested. 
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Chapter  2 

Elevated  Temperature  Interfacial 
Behavior  of  MMCs:  A 
Computational  Study 


2.1  Introduction 


Metallic  and  Intermetallic  Matrix  Composites  (MMCs  and  IMCs)  are  being  considered  for 
crucial  applications  in  the  aerospace  industry  because  of  their  high  specific  stiffiiess  and 
strength  even  at  elevated  temperatures  (>  600°C).  A  critical  issue  in  the  successful  applica¬ 
tion  of  these  composites  is  the  behavior  of  fiber-matrix  interfaces.  Push-out  test  is  emerging 
as  an  important  experimental  tool  for  characterizing  the  interfacial  behavior  of  this  class  of 
composites,  mainly  after  its  widespread  use  in  Ceramic  Matrix  Composites  (CMCs). 


When  MMCs  and  IMCs  are  cooled  after  processing  from  a  stress  free  temperature,  residual 
stresses  arise  due  to  the  difference  in  the  coefficients  of  thermal  expansion  (CTE)  of  fiber 
and  matrix  materials.  These  stresses  reach  considerable  levels  because,  in  typical  titanium 
based  MMCs  and  IMCs  such  as  SiC/Ti-6Al-4V,  SiC/Ti-15V-3Sn-3Al,  SiC/Ti-24Al-llNb 
the  mismatch  in  the  CTE  is  quite  significant  {amatrix  >  "^ctfiber)  and  also  the  temperature 
differential  involved  in  the  cooling  process  is  high  (about  800  "C)  [3].  These  stresses  cause 
matrix  cracking  near  the  interface  on  processing  and  also  affect  the  interfacial  response  of 
these  composites  under  service  conditions. 


Push-out  tests  involve  application  of  compressive  load  through  a  flat  indentor  to  .push  a 
fiber  out  of  a  composite  specimen.  A  schematic  of  the  experimental  set-up  and  a  typical 
push-out  curve  are  shown  in  Fig.  2.1.  The  cross  head  displacement  is  plotted  along  the 
X-axis.  Applied  load,  P  or  average  interfacial  shear  stress,  r  is  plotted  along  the  y-axis.  t  is 
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evaluated  from  P  using  the  relationship 


where  d  is  the  diameter  of  the  fiber,  and  L  is  the  length  of  the  specimen.  Eldridge  et  al.  [9] 
conducted  thin-slice  push-out  tests  of  SiC/Ti-15-3  and  SiC/Ti-24-11  composites  (Fig.  2.2) 
at  various  temperatures.  In  this  figure,  the  average  shear  stress  evaluated  from  equation  (1) 
is  plotted  against  the  cross  head  displacement.  Pmax  (or  r^)  and  Pjr  (or  t^)  are  the  two 
parameters  reported  from  the  push-out  tests. 


The  higher  bond  strength  of  MMCs  and  IMCs  compared  to  CMCs  limit  the  specimen  thick¬ 
ness  for  push-out  tests  in  MMCs  to  a  value  less  than  1.0  mm  (with  thickness  less  than  0.5  mm 
being  the  most  common),  due  to  experimental  problems  such  as  fiber  damage  during  load¬ 
ing  or  failure  of  the  indentor  [5].  Previous  work  by  Ananth  and  Chandra  [3]  indicates  that 
processing  induced  residual  stresses  considerably  influence  the  debonding  sequence  during 
push-out  tests  at  these  thicknesses  and  makes  the  interpretation  of  the  test  results  difficult. 
Nonlinear  finite  element  method  was  used  in  that  work  to  model  the  propagation  of  interface 
failure  during  push-out  tests.  They  showed  that  the  intrinsic  shear  strength  of  the  interface 
can  be  very  much  different  than  from  the  experiment;  and,  developed  a  methodology  to 
evaluate  the  interfacial  properties  by  correlating  the  numerical  predictions  with  the  experi¬ 
mental  data.  In  this  work,  the  comprehensive  analysis  procedure  developed  by  the  authors 
[3,  7]  is  used  to  study  the  interfacial  behavior  during  push-out  tests  at  elevated  temperatures. 


2.2  Modeling  the  single  fiber  push-out  test 


The  single  fiber  push-out  test  is  analyzed  using  axis’>-mmetric  cylindrical  models.  A  two  phase 
unit  cell  consisting  of  a  fiber  and  matrix  is  used  with  periodic  boundary  conditions  at  the 
outer  radius  (rv  =  constant)  to  represent  the  composite.  The  interface  is  modeled  as  a  region 
of  negligible  thickness  [7,  8,  9,  26]  possessing  the  required  stiffness  and  strength.  A  stress 
based  criterion  for  debonding,  and  frictional  resistance  based  criterion  for  interfacial  sliding 
are  used  to  capture  debonding  and  sliding  which  are  presumed  to  occur  in  a  chronological 
sequence  with  increasing  load.  Debonding  is  postulated  to  occur  under  the  combined  action 
of  normal  tensile  stress  which  could  cause  mode  I  failure  and  shear  stress  which  could  cause 
mode  II  failure  at  the  interface.  A  quadratic  stress  based  failure  theory 


is  applied  [25,  26]  where  cr/  is  the  interfacial  strength  in  the  normal  direction  resisting  craok 
opening,  and  is  the  effective  shear  strength  of  the  interface.  In  the  above  equation  the 
effect  of  cTr  is  included  only  when  it  is  positive.  Subsequent  to  the  debonding  process, 
interfacial  frictional  sliding  occurs  in  the  mating  debonded  surfaces  according  to  Coulomb’s 
law. 
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Any  contact  pressure  at  the  interface  will  result  in  an  increase  in  the  value  of  t^,  which  is 
then  evaluated  as 

=  r*  +  HP  with  p  =  —Or  if  Cr  <  0  (2-3) 

p  =  0  if  cTr  >  0 

where  t*  is  the  shear  strength  of  the  interface  (the  shear  stress  required  to  break  the  bond), 
p  is  the  contact  pressure  at  the  interface,  and  p,  is  the  coefficient  of  friction,  a}  is  taken  to 
be  equal  to  be  of  the  same  order  as  r*  [26]  in  the  absence  of  any  experimental  data. 


2.3  Numerical  Implementation 


The  problem  is  analyzed  by  assuming  axial  symmetry.  The  interface  is  modeled  using  a 
contact-friction  formulation  [7]  and  the  fiber-matrbc  bonding  is  simulated  by  connecting  the 
nodes  at  the  interface  by  means  of  high  stiffness  springs.  MARC  finite  element  software  is 
used  to  perform  the  analysis  along  with  external  subroutines.  The  stresses  at  the  interface  are 
monitored  in  these  subroutines  and  the  stiffness  of  the  springs  is  dropped  to  a  negligible  value 
once  the  failure  criterion  is  satisfied,  indicating  the  onset  of  frictional  sliding.  Coulomb’s 
law  is  used  to  model  the  frictional  behavior.  The  coefficient  of  friction  p  is  approximately 
calculated  from  the  relation  P  =  where  is  the  average  frictional  stress  from  the  push- 
out  experiment,  and  Pres  is  the  average  residual  contact  pressure  (at  the  interface)  from  the 
numerical  results. 


Elastic  constitutive  behavior  is  assumed  for  the  fiber  phase.  Matrix  is  assumed  to  be  a 
rate  independent  elastic-plastic  material.  The  temperature  dependency  of  the  elastic  and 
inelastic  properties  of  the  constituent  phases  are  included  in  the  analysis  in  a  piece- wise  linear 
manner,  for  different  titanium  based' matrices  and  silicon  carbide  (SCS-6)  fiber  [3,  21,  13). 

To  simulate  the  thin-slice  push-out  test,  the  analysis  is  done  in  three  steps.  The  boundary 
conditions  for  these  steps  are  shown  in  Fig.  4.5. 

1.  In  the  first  step  (Fig.  4.5a  )  the  cooling  process  after  composite  consolidation  at  high 
temperature  is  modeled.  A  reference  temperature  Tref  is  assumed  above  which  the 
composite  is  stress  free. 

2.  The  second  step  (Fig.  4.5b)  models  the  specimen  preparation  process.  This  is  sim¬ 
ulated  by  removing  the  existing  tying  constraints  and  boundary  conditions  from  the 
top  and  bottom  faces  and  allowing  the  stresses  in  the  specimen  to  reach  equilibrium, 
symmetrical  about  a  plane  passing  through  the  center  of  specimen  thickness.  This 
process  results  in  shear  stresses  at  the  interface  due  to  differential  axial  residual  strains 
between  the  fiber  and  matrix. 
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3.  The  push-out  load  is  applied  by  prescribing  displacements  to  a  rigid  punch,  till  the 
fiber  is  pushed  through.  The  boundary  conditions  used  in  this  step  are  shown  in  Fig. 
4.5(c).  The  clearance  between  the  supporting  hole  and  the  fiber  diameter  is  made  as 
smaJl  as  possible  to  minimize  bending. 


The  interface  failure  criterion  is  checked  during  steps  2  and  3  to  detect  the  possibility  of 
debonding.  Hence,  if  the  residual  shear  stresses  are  high  and  the  bonding  at  the  interface  is 
weak,  then  debonding  can  occur,  even  before  the  application  of  punch  load. 


2.4  Results  and  Discussion 


The  push-out  test  is  numerically  modeled  with  SiC/Ti-15-3  MMC  as  the  model  material 
system.  The  silicon  carbide  fiber  used  in  this  study  is  the  SCS-6  type  (continuous  fibers) 
with  140/im  fiber  diameter.  The  reference  temperature  at  which  the  composite  is  stress  free 
is  taken  to  be  [13].  Fiber  volume  fraction  of  35%,  specimen  thickness  of  0.45  mm  and 

an  indentor  of  diameter  are  used  so  that  the  numerical  predictions  can  be  compared 

with  the  experimental  results  of  Eldridge  et  al.  [9]. 


2.4.1  Interface  failure  at  room  temperature 

The  3-step  numerical  procedure  outlined  in  the  previous  section  is  used  to  simulate  the 
interfacial  failure  process  during  push-out  test.  Residual  stress  state  in  the  bulk  composite 
after  the  cooling  process  is  shown  in  Fig.  2.4  (step  1).  It  can  be  seen  that  axial  compressive 
stresses  {o^)  of  the  order  of  750  MPa  is  produced  in  the  fiber  and  a  corresponding  tensile 
stress  is  produced  in  the  matrix  due  to  the  thermal  mismatch  between  the  fiber  and  the 
matrix.  The  radial  stress  (a^)  is  compressive  at  the  interface  and  acts  as  a  clamping  force 
on  the  fiber.  Step  2  which  models  the  specimen  preparation  process  results  in  considerable 
amount  of  shear  stresses  at  the  ends  (Fig.  2.5)  due  to  the  removal  of  boundary  constraints. 
The  effect  of  this  shear  stress  is  to  push  the  fiber  (which  is  under  axial  compression)  out  of 
the  matrix  at  both  the  ends.  Adequate  modeling  of  the  end  effects  through  mesh  refinement 
reveals  that  the  end  points  of  the  specimen  debond  due  to  the  combined  effect  of  the  high 
levels  of  shear  stresses  and  radial  normal  stresses  (resulting  from  the  local  bending  of  the 
matrix).  When  an  axial  load  is  applied  to  the  fiber  under  these  conditions,  debonding  starts 
from  the  far  end  (or  support  end)  and  progresses  towards  the  punch  end  (Fig.  7.15).  Fig. 
7.15a  shows  the  residual  shear  stress  state  and  (b),  (c),  (d)  and  (e)  are  the  stress  states  at 
increasing  punch  displacements.  It  can  also  be  seen  from  Fig.  7.15  that  the  fiber  protrudes 
outside  in  the  axial  direction  at  the  support  end,  as  the  residual  axial  compressive  stresses 
are  relieved  in  the  debonded  region.  This  is  contrary  to  the  normally  observed  debonding 
behavior  in  CMCs  where  the  interface  failure  always  initiates  at  the  loading  end.  A  detailed 
study  of  this  phenomena  can  be  found  in  reference  [8,  3]. 
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2.4.2  Elevated  temperature  push-out  behavior 

Residual  stresses  decrease  Avith  increase  in  temperature.  This  results  in  a  decrease  in  the 
axial  strain  mismatch  between  the  fiber  and  matrix  phases  thus  reducing  the  residual  shear 
stress  at  the  interface  as  shown  in  Fig.  2.7.  The  push-out  test  is  numerically  simulated  at 
23°C  and  400°C  assuming  a  constant  value  of  interfaeial  shear  strength  (r*  =  400.0  MPa)  in 
this  temperature  raiige,  and  the  results  are  shown  in  Fig.  2.8.  We  can  see  from  this  figure 
and  from  Fig.  2.9  that  the  peak  load  Prnax  increases  with  temperature.  This  is  because 
the  reduction  in  shear  stresses  at  elevated  temperatures  requires  a  larger  load  for  failure 
initiation  from  the  support  end.  Fig.  2.9  also  shows  the  plot  of  debond  initiation  load  Pi 
with  temperature.  The  debonding  sequence  indicates  that  the  failure  under  these  conditions 
of  temperature  and  thickness,  still  initiates  from  the  support  end.  It  can  be  seen  from  Fig. 
2.9  that  the  difference  between  and  Pi  decreases  with  increase  in  temperature.  This 
can  be  attributed  to  the  reduction  in  residual  shear  stresses  at  the  loading  end  at  elevated 
temperatures  which  causes  the  resistance  to  further  debonding  to  decrease  as  the  crack 
propagates  towards  the  loading  end. 


Apart  from  the  shear  stresses,  radial  clamping  stresses  at  the  interface  also  decrease  with 
increase  in  temperature  as  shown  in  Fig.  2.10.  The  effect  of  this  drop  can  be  seen  from -the 
significant  reduction  in  Tfr  from  23‘’C  to  400°C  in  Fig.  2.8.  The  coefficient  of  friction  /i  is 
taken  as  0.25  in  these  calculations. 


2.4.3  Evaluation  of  interfacial  properties 


The  experimental  data  of  at  different  temperatures  for  SiC/Ti-15-3  composite  is  available 
from  Fig.  2.2  [9].  Tut  is  only  an  average  value  and  is  not  a  measure  of  the  actual  shear 
strength  of  the  int-iiace  [7].  The  intrinsic  shear  strength  of  the  interface  which  is  a  material 
property,  could  be  higher  than  the  which  is  computed  from  the  peak  load  from  push- 
out  test,  using  equation  (1).  Thus  one  of  our  main  objective  is  to  establish  a  quantitative 
relationship  between  the  experimental  value  and  the  intrinsic  shear  strength  r*,  for 
MMCs  and  IMCs.  In  the  present  study,  the  average  interfacial  shear  stress  corresponding 
to  the  peak  load  r^p,  is  evaluated  from  finite  element  simulation  for  different  values  of  t* 
(Fig.  2.11).  This  serves  as  a  calibration  curve  from  which  the  material  shear  strength  r*  can 
be  evaluated.  We  can  now  define  an  Interface  Strength  Factor,  Ks 


(2.4) 


where  r*  is  the  intrinsic  shear  strength  of  the  interface  and  r^-p  is  the  average  interfacial 
shear  stress  corresponding  to  peak  push-out  load  Pmax  from  the  experiment.  Kg  is  found 
to  be  in  the  range  of  4.0  from  room  temperature  to  400°C  and  drops  to  about  2.0  at  about 
700“C. 
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2.4.4  Support  distance  and  its  effect  on  push-out  test  results  at 
different  temperatures 

Kallas  et  al.  [5]  found  that  in  niobiura/sapphire  MMC  system  specimen  bending  under  push- 
out  load  is  a  problem,  especially  at  very  low  aspect  ratios  L/h,  where  L  is  the  thickness  and 
h  is  the  distance  between  the  bottom  supports  in  Fig.  2.1.  This  results  in  tensile  normal 
stresses  at  the  bottom  surface  which  causes  interface  failure  to  initiate  from  the  bottom  even 
in  the  absence  of  residual  shear  stress.  To  determine  the  effect  of  support  distance  in  the 
material  system  under  consideration,  the  numerical  simulation  is  repeated  at  different  L/h 
ratios  ranging  from  2.2  to  3.3.  The  results  indicate  that  the  peak  push-out  load  reduces 
with  reduction  in  L/h  ratio;  but  the  effect  is  significant  only  at  higher  temperatures.  The 
variation  of  the  peak  push-out  load  is  found  to  be  about  3%  at  temperatures  less  than  400°C; 
but,  increases  to  15%  at  700°C.  This  is  due  to  the  reduction  in  the  mechanical  properties  of 
the  material  at  elevated  temperatures. 


2.5  Conclusions 


Residual  shear  stresses  affect  the  peak  push-out  load  as  well  as  the  failure  initiation  load 
at  elevated  temperatures.  The  average  shear  strength  iTacp  —  Pmax/'^dL)  evaluated  from 
the  peak  push-out  load  does  not  reflect  the  actual  shear  strength  of  the  interface,  since  it  is 
influenced  by  the  presence  of  processing  induced  residual  stresses  and  other  geometric  and 
testing  variables.  Numerical  predictions  show  that  for  SiC/Ti-15-3  the  intrinsic  interfacial 
shear  strength  (which  is  a  material  property)  is  about  four  times  at  room  temperature 
and  drops  to  2.0  at  TOO^C.  The  support  distance  in  the  experimental  set-up  has  to  be 
controlled  properly  at  elevated  temperatures  to  minimize  the  specimen  bending  effects. 
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Interfacial  Shear  Stress,  MPa 


Figure  2.1:  Typical  push-out  test  curve  with  the  schematic  of  the  experimental  set-up 


(a)  SiCn’i-24-11  (b)  SiC/Ti-15-3 

Figure  2.2:  Push-out  experimental  results  for  SiC/Ti-15-3  and  SiC/Ti-24-11  composites 
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(a)  (b)  (c) 


Figure  2.3:  Boundary  conditions  used  in  the  finite  element  analysis 
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Figure  2.4:  ResiduaJ  stresses  after  step  1;  Material:  SiC/Ti-15-3 
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Figure  2.5;  Residual  shear  stress  distribution  alon 
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Figure  2.6:  Propagation  of  interface  failure  in  thin-slice  push-out  test 
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Interfacial  Shear  Stress,  MPa 
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Fiber  Displacement,  mm 

Figure  2.8:  Predicted  push-out  behavior  at  T  =  23^  C  and  T  =  400^  C;  r*  =  400.0MFa 


Temperature,  C 

Figure  2.9:  Variation  of  failure  initiation  load  with  temperature;  r*  =  400.0MPa 
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Relative  Distance 

Figure  2.10:  Variation  of  residual  radial  stresses  along  the  interface  with  temperature 
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Figure  2.11:  Calibration  curve  to  obtain  interfacial  shear  strength  from  experimental  data: 


Chapter  3 

Effect  of  Interface  Chemistry  on  the 
Fracture  Properties  of  Titanium 
Matrix  Composites 


3.1  Introduction 

Metal  matrix  composites  (MMC)  are  being  considered  as  the  material  systems  of  choice  for ' 
the  jet  engines  of  next  generation  aircraft,  owing  to  their  high  specific  stiffness  and  strength. 
In  MMCs  titanium  matrices  reinforced  by  silicon  carbide  fibers  have  shown  considerable 
potential  for  high  temperature  aerospace  applications.  A  critical  issue  to  be  addressed  in  the 
use  of  Titanium  matrix  composites  is  the  high  reactivity  between  the  fiber  and  the  matrix 
phases,  which  leads  to  the  formation  of  a  reaction  zone  at  the  fiber/matrix  interface.  The 
cnemistry  and  geometry  of  this  zone  can  influence  fiber- matrix  interfacial  properties,  and 
hence  the  overall  thermomechanical  response  of  these  Titanium  based  MMCs. 

SiC  fibers  are  of  considerable  interest  as  reinforcements  for  metal  and  ceramic  matrix  com¬ 
posites  because  of  their  high  strength  and  stiffness  even  at  elevated  temperatures.  In  MMCs 
a  reaction  layer  often  forms  between  the  outermost  regions  of  the  fiber  and  the  adjacent 
matrix.  This  reaction  layer  constitutes  the  interfacial  region.  The  microstructure  of  the 
outer  region  of  the  fiber  and  the  fiber /matrix  interphase  play  a  vital  role  in  determining  the 
properties  of  the  composite.  The  microstructure  of  the  SCS-6  fiber  (developed  by  Textron 
Specialty  Materials  Inc.)  has  been  reported  in  the  literature  [1].  The  SCS-6  fiber  consists 
of  several  layers  deposited  using  CVD  technique,  as  schematically  shown  in  Figure  3.1.  The 
fiber  comprises  of  three  broad  regions  namely  a  carbon  core,  many  silicon  carbide  layers 
and  an  outer  coating  rich  in  carbon.  The  outer  coating  which  is  the  region  of  interest,  is  a 
carbonaceous  layer  comprising  carbon  matrix  with  SiC  particles  embedded  in  it.  This  outer 
(carbon)  coating  itself  is  composed  of  three  distinct  layers  (layersl,  2  and  3)  with  varying  size 
and  density  of  SiC  particles.  Layer  2  is  very  thin  compared  to  the  other  two  and  also  has  the 
lowest  density  of  silicon  carbide  particles  in  it  [1].  Push-out  experiments  on  SCS-6/RBSN 
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CMCs  [1]  have  shown  that  the  failure  initiation  occurs  predominantly  in  this  region.  Similar 
experiments  by  Roman  and  Jero  [2]  on  as-processed  SCS-6/Ti-6Al-4V  MMC  also  revealed 
that  the  interface  failure  has  occurred  in  the  outer  carbon  coating  and  not  in  the  reaction 
zone. 

Titanium  is  highly  reactive  at  elevated  temperatures  reacts  with  the  carbon  and  silicon  in 
the  coating  to  form  the  reaction  zone  (Figure  3.2).  These  constituents  (carbon,  silicon  and 
titanium)  inevitable  in  the  Titanium  Matrix  Composite  systems  ;  hence,  a  clear  understand¬ 
ing  of  the  reaction  between  these  elements  at  elevated  temperatures  and  the  effect  of  the 
reaction  products  on  the  mechanical  properties  is  essential.  Studies  of  reaction  zone  forma¬ 
tion  in  TMCs  have  been  carried  out  by  several  researchers  [3,  4,  5].  None  of  these  studies 
focus  on  establishing  a  relationship  between  the  mechanic^  properties  and  the  chemical 
composition  of  the  interfacial  region  as  the  interface  evolves.  This  is  the  objective  of  the 
current  study.  This  tjrpe  of  study  is  critical  in  not  only  understanding  the  evolution  of  the 
interfacial  reaction  zone,  but  also  in  providing  means  to  optimize  the  interface  chemistry  for 
performance. 

Fiber  Push-out  Test  has  recently  emerged  as  a  popular  test  for  characterizing  the  mechanical 
properties  of  the  interface  due  mainly  to  simpler  specimen  preparation  and  equipment.  This 
test  has  been  used  to  measure  the  interfacial  properties  of  MMCs  at  room  [6,  1,  8]  and 
elevated  temperatures  [9,  10].  Chandra  and  coworkers  [3,  7,  5]  analyzed  the  thin-slice  piish- 
out  test  and  simulated  the  initiation  of  interface  failure  and  its  propagation  using  non-linear 
finite  element  methods.  This  work  was  able  to  explain  many  of  the  experimental  observations 
[9]  of  push-out  behavior  at  room  and  elevated  temperatures.  They  also  proposed  a  method 
for  evaluating  the  interfacial  shear  strength  [14]  and  fracture  toughness  [7]  from  the  push-out 
data. 

The  present  work  involves  elevated  temperature  exposure  experiments  on  MMCs  to  study  the 
effect  of  exposure  temperature  and  time  on  the  reaction  zone,  and  -he  effect  of  evolution  on 
mechanical  properties.  SCS-6/Timetal  21s  (Ti-15Mo-3Al-2.7Nb-0.‘>5Si)  composite  system  is 
selected  as  the  model  material  system  and  the  relationship  between  reaction  zone  evolution 
and  interface  mechanical  properties  is  investigated.  Push-out  tests  in  combination  with 
the  computational  fracture  mechanics  approach  developed  by  the  authors  [5,  7]  are  used  to 
evaluate  the  fracture  properties. 


3.2  Evaluation  of  Interfacial  Fracture  Toughness  of 

MMCs 


In  this  paper,  a  numerical  method  based  on  the  equivalent  domain  integral  technique  devel¬ 
oped  by  the  authors  [5,  7]  is  used  to  simulate  the  fracture  process  during  a  thin  slice  push-out 
test  and  extract  the  interfacial  fracture  toughness  Gic  of  SCS-6/Timetal21s  composites  from 
the  experimental  data. 
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The  thin-slice  push  out  test  involves  pushing  a  fiber  out  of  a  composite  specimen  by  the 
application  of  compressive  load  through  a  flat  indenter.  Modeling  such  a  test  involves  three 
major  stages.  Step  1  involves  modeling  the  composite  consolidation  at  high  temperatures 
and  then  cooling  down  to  room  temperature  which  induces  residual  stresses  in  the  composite. 
Step  2  is  the  specimen  preparation  stage  consisting  of  slicing  a  thin  push-out  specimen  from 
the  bulk  composite  which  results  in  stress  redistribution  at  the  two  ends  of  a  fiber.  Step  3  is 
modeling  the  loading  stage  where  compressive  load  is  applied  to  the  fiber  until  it  is  pushed  out 
of  the  specimen.  A  detailed  description  of  the  above  sequence  of  steps  is  provided  elsewhere 
[3].  The  test  is  analyzed  using  an  axisymmetric  cylindrical  model.  A  three  phase  model 
consisting  of  fiber,  matrix  and  a  composite  phase  with  eflFective  properties  has  been  used. 
The  interface  is  modeled  as  a  contact  surface  and  the  bonding  at  the  interface  is  simulated 
using  high  stiffness  springs.  The  interface  failure  process  is  modeled  using  a  strain  energy 
based  failure  criteria  given  by: 

Gi  ^  Gic 

where  Gi  is  the  strain  energy  release  rate  for  the  interface  crack  and  Gic  is  the  critical  value  of 
the  strain  energy  release  rate.  The  interface  crack  propagates  further  when  its  strain  energy 
release  rate  exceeds  the  critical  value.  The  strain  energy  release  rate  of  the  propagating 
interface  crack  is  computed  using  the  Equivalent  Domain  Integral  (EDI)  [8,  9,  4]  method. 
The  EDI  method  was  implemented  in  the  commercial  finite  element  analysis  software  MARC 
[6]  through  user  subroutines.  A  typical  mesh  used  is  shown  in  Figure  5.1.  A  small  interface 
crack  of  the  order  of  0.5r/  (where  ry  is  the  fiber  radius)  is  introduced  at  both  the  ends 
of  the  specimen  during  step  2  of  the  simulation  for  the  purpose  of  computation  of  energy 
release  rates. 


3.3  Results  and  Discussion 


The  effect  of  elevated  temperature  exposure  on  the  interfacial  fracture  behavior  of  MMCs  is 
studied  by  subjecting  SCS-6/Timetal  21s  specimens  to  selected  temperatures  (450°  C,  700° 
C  and  927°  C)  in  vacuum  for  varying  periods  of  time.  The  time  periods  used  are  25,  70  and 
120  hours.  Most  of  the  specimens  for  the  exposure  study  were  cut  from  a  composite  panel 
consolidated  by  Textron  Specialty  Materials,  Inc.  Some  of  the  specimens  were  fabricated  in- 
house  using  foil-fiber-foil  technique.  Scanning  Electron  Microscope  (SEM)  is  used  to  measure 
the  thickness  of  the  interface.  Interfacial  chemistry  is  monitored  using  Energy  Dispersive 
X-ray  Spectroscopy  (EDX).  Push-out  tests  are  conducted  on  the  heat  treated  specimens  and 
the  results  are  compared  with  those  of  as-processed  ones.  Numerical  simulations  are  used  to 
extract  interfacial  fracture  toughness  values  from  the  experimental  data. 


3.3.1  Change  in  Reaction  Zone  Size  and  Composition 

The  change  in  reaction  zone  thickness  is  measured  as  a  function  of  time  for  two  temperatures 
927°  C  and  700°  C  and  compared  with  the  reaction  zone  size  for  as  processed  conditions. 
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The  presence  of  a  reaction  zone  is  identified  by  the  change  in  the  morphology  in  the  SEM 
micrographs  (Figure  3.2).  The  coating  size  of  unconsolidated  fibers  is  measured  and  taken  as 
the  reference.  The  variation  in  the  reaction  zone  and  the  coating  thicknesses  for  specimens 
exposed  to  927°  C  for  25,  70  and  120  hours  were  measured  and  plotted  in  Figure  3.4  along 
with  the  as-processed  sample.  It  is  seen  from  this  figure  that  the  initial  growth  of  the 
reaction  zone  for  the  25  hours  case  is  by  partial  consumption  of  the  coating  and  partial 
transformation  of  the  matrix  region.  With  extended  exposure  ,upto  120  hours  the  reaction 
zone  growth  occurs  more  through  transformation  of  the  matrix  than  through  consumption  of 
the  coating.  The  exposure  studies  conducted  at  700°  C  showed  that  the  reaction  zone  growth 
with  time  and  also  the  change  in  the  thickness  of  the  outer  carbon  coating  are  significantly 
reduced.  An  important  observation  from  the  exposure  studies  is  that  the  coating  is  not  fully 
consumed  even  after  being  subjected  to  high  temperatures  for  extended  periods  of  time. 

Elemental  composition  and  diffusion  are  studied  with  the  help  of  EDX.  Line  scans  are  per¬ 
formed  across  the  reaction  zone,  showing  the  variation  in  the  chemical  composition.  The 
diffusion  of  alloying  elements  across  the  reaction  zone  is  studied  as  a  function  of  time  for  927° 
C  and  700°  C.  The  results  for  the  exposure  study  at  927°  C  for  120  hours  only  are  discussed 
in  detail  in  this  study  since  the  maximum  change  is  observed  in  this  case.  Comparisons  are 
made  with  the  as-processed  samples.  Since  EDX  system  used  in  this  study  is  not  capable 
of  detecting  lighter  elements  such  as  Carbon,  the  Carbon  percentage  is  deduced  from  the 
relative  counts  as  the  rest  of  the  chemical  composition.  The  concentration  of  Titanium  was 
found  to  be  the  maximum  in  the  reaction  zone,  followed  by  silicon.  The  concentration  of 
the  other  alloying  elements  such  as  Molybdenum,  Aluminum  and  Niobium  are  observed  to 
be  minimal. 

The  variation  of  Titanium  for  the  as-received  specimen  is  shown  in  Figure  3.5.  It  is  seen 
from  this  figure  that  Titanium  content  for  half  of  the  reaction  zone  is  of  the  same  order 
as  that  in  the  matrix.  The  Titanium  content  reduces  sharply  in  the  reaction  zone  near  the 
coating.  Figure  3.6  shows  the  line  scan  of  Titanium  from  a  sample  which  was  exposed  to  o2  7° 
C  for  120  hours.  It  can  be  seen  that  the  diffusion  of  Titanium  towards  the  carbon  coa+  ng 
is  more  profound  after  the  thermal  exposure  when  compared  to  the  as-processed  composite 
specimen.  Figure  3.7  shows  the  variation  of  silicon  across  the  reaction  zone.  From  this  figure, 
it  can  be  seen  that  the  highest  concentration  of  silicon  occurs  at  the  coating  reaction  zone 
interface.  This  silicon  concentration  decreases  across  the  reaction  zone  thickness  towards  the 
matrix.  The  behavior  of  silicon  after  the  thermal  exposure  is  observed  to  be  similar  to  that 
of  the  as-processed  case  as  can  be  seen  from  Figure  3.8. 

It  can  be  inferred  from  the  chemical  composition  studies  that  the  carbon  from  the  fiber 
coating  probably  diffuses  the  maximum  into  the  matrix  indicating  Titanium  Carbide  to  be 
the  compound  which  is  present  in  most  of  the  interfacial  region.  This  has  also  been  confirmed 
by  other  researchers  [3].  The  results  also  suggest  that  Titanium  Silicides  (ThSis)  may  be 
present  mainly  near  the  coating  reaction  zone  interface.  It  has  to  noted  that  TiC  in  the 
reaction  zone  is  formed  in  two  possible  ways;  by  the  diffusion  and  subsequent  reaction  of 
the  free  carbon  (from  the  outer  carbon  coating)  with  the  Titanium,  and  due  to  the  direct 
reaction  of  Ti  with  SiC  from  the  coating. 
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3.3.2  Evaluation  of  Interfacial  Fracture  Toughness  from  Thin- 
slice  Push-out  Tests 


Fiber  Push-out  test  was  carried  out  on  as-processed  and  heat  treated  specimens.  The  tests 
were  carried  out  at  NASA  Lewis  Research  center.  The  average  thickness  of  the  as-processed 
specimens  is  about  0.53  mm.  A  typical  load  displacement  curve  for  a  push  out  test  on  as- 
processed  specimen  is  shown  in  Figure  3.9.  Fro»u  this  figure  it  is  seen  that  after  attaining 
the  peak  load  there  is  a  sharp  load  drop  which  corresponds  to  the  fiber  being  completely 
debonded  and  the  onset  of  frictional  sliding.  Full  debonding  is  also  characterized  by  a  sharp 
acoustic  signal  also  shown  in  the  figure.  The  average  peak  load  obtained  for  the  as-processed 
specimens  is  of  the  order  of  10  N  which  corresponds  to  an  average  peak  debond  shear  stress 
of  about  38  MPa.  Figure  3.10  shows  the  push-out  behavior  after  extended  exposure  at  927° 
C.  It  can  be  seen  that  the  peak  load  in  this  has  significantly  increased.  It  is  also  observed 
from  this  figure  that  there  is  no  distinct  sharp  load  drop  associated  with  complete  debonding 
as  is  seen  for  the  as-processed  case.  There  is  also  no  sharp  acoustic  signal  which  marks  the 
event  of  complete  debonding  as  observed  in  the  earlier  case.  The  push-out  test  results  at 
700°  C  and  at  450°  C  are  observed  to  be  similar  to  that  the  as-processed  case  characterized 
by  abrupt  load  drops  after  complete  debonding  accompanied  by  sharp  acoustic  signals. 

A  composite  of  the  test  results  for  the  as- processed  specimens  and  those  subjected  to  450°, 
700°  and  927°  C  for  different  exposure  time  periods  is  shown  in  Figure  3.11.  Each  point  is 
the  average  of  12  push-out  test  results.  The  scatter  in  the  data  is  also  shown  in  the  figure. 
As  mentioned  earlier  the  average  peak  debond  shear  stresses  are  significantly  higher  for  the 
specimens  subjected  to  927°  C  heat  treatments  than  those  for  the  as-processed  specimens 
and  the  specimens  heat  treated  at  450°  and  700°  C.  The  results  for  the  927°  C  heat  treated 
specimens  also  show  a  large  degree  of  scatter.  It  can  also  be  seen  from  this  figure  that 
the  exposure  time  did  not  have  a  significant  effect  on  the  test  results,  especially  at  lower 
temperatures.  The  results  for  the  tests  conducted  for  the  as  received  specimen  and  those 
subjected  to  the  lower  temperature  ranges  (<  7  j0°  C)  suggest  that  there  is  not  much  change 
in  interfacial  properties  at  these  temperatures.  Thus  the  most  significant  change  in  interfacial 
behavior  is  observed  for  the  specimens  subjected  to  927°  C  heat  treatments. 

The  fracture  mechanics  based  approach  developed  previously  [5,  7]  is  used  here  to  extract 
interfacial  fracture  toughness  values  from  the  push-out  test  results.  Numerical  simulations 
of  push-out  tests  representing  conditions  of  specimens,  exposed  to  927°,  700°  C  and  450° 
foi-  120  hours  and  as- processed  specimens  are  carried  out.  The  fracture  toughness  values 
obtained  from  the  numerical  simulations  are  shown  in  Table  1.  It  is  observed  that  the  fracture 
toughness  values  for  specimens  heat  treated  at  700°  C  and  450°  C  are  not  very  different  from 
that  of  the  as-processed  specimen  and  are  in  the  range  of  55  to  60  J jTr? .  However  when  the 
specimens  were  exposed  at  927°  C,  the  mode  II  fracture  toughness  is  observed  to  actually 
increase  significantly  to  about  70  J/m^. 
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3.3.3  Observation  of  the  Fracture  Surface  of  Push-out  Specimens 


In  order  to  obtain  a  better  idea  about  the  mechanics  of  the  failure  process  the  fracture 
surfaces  of  the  failed  push-out  specimens  were  examined  with  the  help  of  an  SEM.  It  is 
observed  that  in  all  the  cases,  debonding  initiated  from  the  bottom  of  the  push-out  specimen, 
i.e.,  the  side  opposite  to  that  of  the  point  of  application  of  the  compressive  load.  This  is 
consistent  with  theoretical  predictions  [3]  and  experiments  carried  out  in  other  as-processed 
Titanium  Matrix  Composite  systems  where  residual  stresses  are  present  [6,  8].  The  bottom 
surface  was  examined  for  fracture  characteristics  as  crack  initiation  takes  place  there.  Typical 
fracture  surfaces  for  700°  and  927°  C  are  shown  in  Figures  3.12  and  3.13  respectively. 

The  fracture  surfaces  for  the  700°  C  specimens  reveal  that  cracks  initiate  in  different  regions 
of  the  coating,  occurring  more  frequently  in  the  middle  and  in  the  region  near  the  coating 
reaction  zone  interface.  The  nature  of  the  fracture  surfaces  at  700°  C  indicate  a  sharp  abrupt 
fracture  process  thus  consistent  with  the  abrupt  load  drops  seen  in  the  load  displacement 
curves  and  the  sharp  acoustic  emission  results.  The  element  diffusion  results  and  the  negligi¬ 
ble  reaction  growth  at  700°  C  suggests  that  the  middle  region  of  the  coating  is  the  weak  link 
in  these  specimens.  The  fracture  surfaces  for  the  as-processed  specimens  and  those  subjected 
to  450°  C  heat  treatments  are  observed  to  be  similar  to  that  of  the  700  C  case. 

The  fracture  surfaces  of  the  927°  C  heat  treated  specimens  are  significantly  different  from 
that  of  the  700°  heat  treated  specimens.  It  is  observed  from  scanning  electron  micrographs 
that  debonding  initiates  more  in  the  interface  region  than  in  the  coating.  The  fracture 
surfaces  are  rough  with  the  presence  of  large  amount  of  debris.  Cracks  are  seen  to  form 
in  the  middle  of  the  reaction  zone  or  near  the  reaction  zone  matrix  interface.  Comparison 
of  the  fracture  surfaces  of  the  push-out  specimens  exposed  to  different  time  periods  did  not 
reveal  significant  differences  in  the  fracture  behavior.  The  large  reaction  zone  powth  and  the 
extensive  diffusion  of  titanium  into  the  coating  indicates  the  possible  formation  of  titanium 
carbides  and  silicides  near  the  coating.  Also  since  there  is  a  considerable  consumption  of 
the  coating,  the  weaker  region  (layer  2)  in  the  coating  is  replaced  by  the  strongly  bonded 
high  concentration  titanium  region.  The  fracture  process  in  the  927°  C  specimens  thus  is 
characterized  by  a  slow  fragmentation  of  the  reaction  zone  and  high  peak  debond  and  sliding 
frictional  stresses. 


3.4  Summary  and  Conclusions 

The  evolution  of  interface  in  SCS-6/Timetal  21s  composite  system  is  studied  by  means  of 
heat  treatment  experiments,  with  the  objective  of  establishing  a  link  between  the  evolution 
of  chemistry  and  geometry  of  the  reaction  zone  and  the  interfacial  properties. 

•  Microstructural  studies  reveal  that  reaction  growth  is  due  to  transformation  of  both 
coating  and  the  matrix.  The  growth  of  the  reaction  zone  into  the  matrix  is  much  more 
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than  the  growth  towards  the  fiber,  as  carbon  is  able  to  diffuse  much  farther  into  the 
matrix  than  the  diffusion  of  Titanium  into  the  fiber. 

•  Reaction  zone  sizes  increase  rapidly  with  time  at  927°  C.  At  lower  temperatures  such 
as  700°  C  and  450°  C  no  significant  increase  in  reaction  zone  size  is  observed  for  this 
system. 

•  Push-out  tests  conducted  on  heat  treated  specimens  reveal  that  the  effect  of  thermal 
exposure  on  interfacial  properties  is  predominantly  due  to  exposure  temperature.  Time 
of  exposure  is  found  to  have  only  a  secondary  effect. 

•  Fracture  toughness  values  remain  more  or  less  constant  for  normal  use  temperatures 
in  the  range  of  450°  to  700°  C.  The  values  however,  increase  when  the  exposure  tem¬ 
perature  is  raised  to  927°  C. 

•  Fracture  surfaces  of  the  specimens  exposed  at  927°  C  show  that  debonding  took  place 
in  the  reaction  zone  whereas  in  as- processed  specimens  and  those  heat  treated  at  the 
lower  temperature  ranges  debonding  initiated  in  the  coating.  This  is  because  the  weak 
region  in  the  coating  is  replaced  by  the  stronger  interfacial  region. 
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Material  System: 

SCS-6/Timetal  21s 

fiber  volume  fraction=  0.35 

average  specimen  thickness=0.53  mm 

Peak  Load 
(Experimental) 

(N) 

^ic 

(Simulation) 

(J/m2) 

As-processed 

10.12 

53.5 

Heat  treated  at  450oC  for  25  hours 

12.01 

55.0 

Heat  treated  at  700oC  for  120  hours 

16.37 

58.0 

Heat  treated  at  921°C  for  120  hours 

29.28 

70.0 

Table  3.1:  Interfacial  fracture  toughness 
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Silicon  Carbide  Layers 


Figure  3.2:  Interfacial  region  in  a  SCS-6/Timetal  21s  composite 
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nF.TAILS  OF  THE  INTERFACE 
1 


•  2' 

Springs  connecting 
the  interface  nodes 


1-1’,  2-2' 

duplicate  nodes 

(with  the  same  co-ordinates) 

at  the  interface 

1-2,  l'-2’ 

Contact  edges  which  cannot 
penetrate 

n,  t 

normal  and  tangential 
directions 

kn.kt 

Stiffness  of  springs  in  normal 
and  tangential  direction^ 

Figure  3.3;  Finite  element  model  and  details  of  the  interface 
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As  processed  specimen  ^ 


Coating 

RZ 

Matrix 

3  |i,m 

0123456  78 

Distance  from  beginning  of  coating  (microns) 


Specimen  subjected  to  927  C  for  25  hours 


Coating 

RZ 

Matrix 

2.47  jxm 

2  |J.m 

0  1  2 

1  f  ^  —\ 

3  4 

1  1-  1  •  1  • 

5  6  7  8 

Distance  from  beginning  of  coating  (microns) 
Specimen  subjected  to  927  C  for  70  hours 


Coating 
2.29  p.m 

RZ 

2.2  p.m 

Matrix 

T — ' — [—■ — H 
0  1  2 

•  1  ‘  1 

3  4 

1  1  j  t  , 

5  6  7  8 

Distance  from  beginning  of  coating  (microns) 
Specimen  subjected  to  927  C  for  120  hours 


Coating 
2.14  pm 

RZ 

3.1  pm 

Matrix 

T — ' — r— ^ — r 
0  1  2 

1  .  ,  1 
3  4 

5 

*  1  ■  1  •  1 

6  7  8 

Distance  from  beginning  of  coating  (microns) 

Figure  3.4:  Variation  in  reaction  zone  size  with  exposure  time  for  SCS-6/Timetal  21s  com¬ 
posite  heat  treated  at  927°  C 


Figure  3.5:  Titanium  x-ray  scan  of  an  as-processed  SCS-6/Timetal  21s  specimen  (a)  line 
scan  (b)  dot  map 


Figure  3.6:  Titanium  x-ray  scan  of  a  SCS-6/Timetal  21s  specimen  exposed  to  927  C  for  120 
hours  (a)  line  scan  (b)  dot  map 
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Figure  3.7:  Silicon  x-ray  scan  of  an  as-processed  SCS-6/Timetal  21s  specimen  (a)  line  scan 
(b)  dot  map 


Figure  3.8:  Silicon  x-ray  scan  of  a  SCS-6/Timetal  21s  specimen  exposed  to  927° 
hours  (a)  line  scan  (b)  dot  map 


C  for  120 


Figure  3.9:  Typical  load  displacement  curve  for  a  push-out  test  conducted  on  an  as-processed 
specimen 


Figure  3.10:  Typical  load  displacement  curve  for  a  push-out  test  conducted  on  a  specimen 
subjected  to  92"^  C  for  25  hours 
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Figure  3.11*  Composite  of  the  push-out  test  results 
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Figure  3.12:  Bottom  fracture  surface  for  a  700°  C  heat  treated  push-out  specimen 
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Figure  3.13:  Bottom  fracture  surface  for  a  927°  C  heat  treated  push-out  specimen 
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Chapter  4 

Effect  of  Residual  Stresses  on  the 
Interfacial  Fracture  Behavior  of 
MMCs 


4.1  Introduction 


Metal  matrix  composites  (MMC)  are  being  considered  as  potential  material  systems  of  choice 
for  the  advanced  propulsion  systems  of  twenty-first  century  aircraft,  owing  to  their  high  spe¬ 
cific  stiffness  and  strength.  One  of  the  critical  issues  in  the  successful  application  of  these 
composites  is  the  behavior  of  the  fibre-matrix  interface.  The  role  of  interface  is  critical  to 
the  performance  of  the  composites.  The  strength,  stiffness  and  fracture  behavior  of  the 
composites  depend  or;  the  interfacial  conditions.  In  MMCs  experimental  investigations^ 
give  clear  evidence  premature  interfacial  debonding,  emphasizing  the  need  to  study  the 
interfacial  behavior  of  this  class  of  composites.  The  push-out  test  is  emerging  as  an  impor¬ 
tant  experimental  tool  for  characterizing  the  interface  behavior  of  MMCs  mainly  after  its 
widespread  use  in  the  study  of  interfaces  in  ceramic  matrix  composites  (CMC).  The  pres¬ 
ence  of  high  levels  of  residual  stresses  in  MMCs  alters  the  mechanics  of  the  interfacial  failure 
process,  making  the  interpretation  of  test  results  difficult.  Moreover  the  relatively  high  bond 
strength  of  MMCs  and  residual  clamping  stresses  at  the  interface  limit  the  specimen  thick¬ 
ness  for  push-out  tests.  The  experiments  under  these  conditions  indicate  that  interface 
failure  initiates  from  the  bottom  of  the  specimen  and  not  from  the  loading  end  as  in  the  case 
of  CMCs.  The  cause  of  bottom  initiation  during  thin  slice  push-out  tests  of  MMCs  has  been 
dealt  in  detail®’®  and  is  found  to  be  predominantly  due  to  residual  stresses.  This  is  different 
from  two  way  debonding  reported  in  the  analyses  of  fiber  pull-out  tests^’®  on  polymer  and 
ceramic  matrix  composites,  where  fiber  volume  fraction  and  Young’s  modulus  ratio  are  the 
dominant  factors. 

In  recent  years  there  have  been  several  efforts  to  model  the  push-out  test  analytically  with 
the  objective  of  understanding  the  interfacial  behavior®"^®.  The  most  comprehensive  study 
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thus  far  which  is  applicable  for  MMCs  is  by  Liang  and  Hutchinson^®;  but  this  work  still  does 
not  address  the  case  of  failure  initiation  from  the  bottom  surface.  Chandra  et  al®’®’^^  analyzed 
the  thin-slice  push-out  test  where  the  initiation  of  interface  failure  and  its  propagation  were 
simulated  using  non-linear  finite  element  methods.  This  work  was  able  to  explain  many  of 
the  experimental  observations^®  of  push-out  behavior  at  room  and  elevated  temperatures. 
The  failure  criteria  used  in  their  approach  is  based  on  local  stresses  and  hence  is  affected  by 
stress  concentration  effects  which  are  a  function  of  mesh  density  in  the  finite  element  analysis. 
This  error  was  minimized  to  a  certain  extent  by  adopting  proper  averaging  schemes  near  the 
crack  tip.  A  more  consistent  quantity  for  monitoring  the  fracture  process  is  the  strain 
energy  release  rate  (Gi)  of  the  interface  crack.  Gic,  the  critical  value  of  Gi  is  a  measure  of 
the  fracture  toughness  of  the  interface. 

In  this  paper,  a  numerical  method  based  on  the  equivalent  domain  integral  technique  has 
been  developed  to  compute  the  strain  energy  release  rates  Gi  of  the  interfaces  in  composites. 
This  method  has  been  applied  to  evaluate  the  Gi  of  as-processed  titanium  based  MMCs  and 
compared  with  some  CMC  systems.  The  interface  failure  process  in  a  thin  slice  push-out  test 
is  simulated  using  the  developed  methodology.  It  is  shown  that  processing  induced  residual 
stresses  significantly  affect  the  initiation  and  propagation  of  interface  cracks  in  MMCs. 


4.2  Evaluation  of  Energy  Release  Rates 


Several  methods  are  available  to  calculate  the  energy  release  rate  G  for  an  existing  crack 
such  as  virtual  crack  closure  technique  (VCCT),  contour  integral  and  the  equivalent  domain 
integral  method  (EDI).  The  contour  integral  method  is  relatively  independent  of  mesh  den¬ 
sity,  however  it  requires  a  focussed  mesh  i.e.,  a  ring  of  singular  elements  around  the  crack 
tip.  This  hampers  crack  propagation  studies.  The  EDI  method  does  not  require  the  use  of 
sii  ^ular  elements  for  computing  energy  release  rates  and  is  ideally  suited  for  studying  the 
pi  vjpagation  of  cracks.  J-integral  is  a  closed  contour  integral  of  the  strain  energy  density  and 
the  work  done  by  tractions  around  the  crack  tip.  It  is  a  path- independent  parameter  and  is 
equivalent  to  the  rate  of  change  of  total  potential  energy  with  reference  to  the  crack  length. 
For  the  purposes  of  numerical  implementation  the  J-integrals  are  converted  to  equivalent 
area  or  domain  integrals^®“^\  hence  the  name  equivalent  domain  integral  (EDI).  The  con¬ 
version  of  line  integrals  to  domain  integrals  is  very  advantageous  because  all  the  quantities 
necessary  for  computation  of  the  domain  integrals  are  readily  available  in  a  finite  element 
analysis.  The  EDI  method  has  been  shown^^  to  give  accurate  results  for  the  J-integral  for 
mode  I,  mode  II  and  mixed  mode  problems.  For  an  arbitrary  closed  contour  F  around  the 
crack  tip  (Figure  4.1),  the  J-integral  is  defined  in  the  absence  of  any  body  forces  as: 


f  r  dui 

-'•*  =  Jr  r 


dr 


•(4.1) 
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(4.2) 


where  k=l,2  and  W  is  the  total  strain  energy  density  defined  as: 

W  =  f  '  CTijdeij 

Jo 

The  integrals  and  are  two  path-independent  integrals  that  compute  the  tota,l  aniount 
of  the  energy  flux  leaving  the  contour  r  in  the  two  directions  Xi  and  0:2  respectively.  is 
known  as  the  J-integral.  Jx^  is  called  the  product  integral.  A  closed  contour  integral  along 
DEFCBAD  (Figure  4.1)  can  be  defined  bv  adding  and  subtracting  the  line  integrals  on  the 

crack  faces  FC  and  AD  and  can  be  written  as: 

J  =_  /■  QSc[r+  f  QSdT+  I  QSdr+  fQdr  +  J  Qdr  (4.3) 

J  DEFCBAD^  JFC  JAD  JCO  JOA 


Jx  —  f  QSdT  -t-  {Jxk)i  ine  (4.4) 

*  Jdefcbad 


where 


dui 

dxk  J 


(4.5) 


For  implementation  in  a  finite  element  code  the  above  contour  integral  on  DEFCBAD  for  J 
can  be  converted  into  a  domain  or  area  integral  as. 


{^Jxi^  domain 


dui  OS' 


f  \dW 


Hence  the  total  integral  Jx^  can  be  written  as 

t/x(t  ~  lain  "F  {Jxi^)line  (^•'^) 

g  _  5'(xi,a:2)  is  an  arbitrary  but  continuous  function  and  has  a  value  of  O  on  the  outer 
contour  and  a  value  of  1  on  the  inner  contour.  In  order  to  be  able  to  separate  the  indi¬ 
vidual  modes  (J/  and  Jji)  from  the  domain  integral,  a  procedure  called  the  decomposition 
method^^"^^  was  adopted  in  which  the  stress  and  displacement  fields  are  separated  into  sym¬ 
metric  and  anti-symmetric  parts.  Consider  two  points  P(a:i,X2)  and  P  (xi,  — X2)  that 
the  immediate  neighbourhood  of  the  crack  tip  and  are  symmetric  about  the  crack  line.  Then 
the  symmetric  and  anti-symmetric  parts  of  the  displacements  can  be  expressed  m  terms  of 
the  displacements  at  points  P  and  P’  as 

f  1  =  1  /  I  (4-8) 

lti2Js  2\u2p-U2P'i 

f  ^  i  (  \  (4.9) 

W2  J  2  I  U2P  +  U2P'  ) 
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Similarly,  the  symmetric  and  anti-symmetric  components  of  the  stresses  can  be  expressed  in 
terms  of  the  stresses  at  points  P  and  P’  as 

"  O’!!  1  f  ' 

^  <T22  I  ^  1  J  <^22P  +  (^22P'  I  (4,10) 

(T33  2  0’33P  +  1^33^' 

,  Cri2  J  5  I  <7i2P  -  cri2P'  . 

'  cm  ^  r  <^11P  ~  ^11^’'  ' 

^  <722  1  _  i  J  <^22P  -  '722P'  I  (4,11) 

(733  I  2  CTssp  —  Cr33P' 

.  C^12  J  AS  ^  ^12P  +  <7i2P'  , 

These  symmetric  and  anti-symmetric  displacements  and  stresses  can  be  used  to  evaluate  the 
four  integrals  Js.,,  Jsx,,  Jasx,,  Jasx,  using  eqn.(7).  The  integrals  Jsx,  and  Jasx,  (domam 
and  line  components  individually)  will  be  identically  zero  because  of  the  symmetric  and 
anti-symmetric  nature  of  the  stress  and  displacement  fields.  The  individual  modes  J;  and 
Jij  can  then  be  written  as 

Jj  =  Jsxi 

Jll  =  JaSx\ 

Jtotal  =  Jsxi  +  JaSxi  (4-12) 


The  line  integrals  Jsx,  and  Jasx,  using  the  symmetric  and  antisymmetric  components  of 
stresses  and  displacements  can  be  written  as 


fC  duo  „  ,  „  ^^2 

Us.,).™  =  2X-  .>22^. 


{Jas  xi)line 
{JASx2)line 


fC  ^  dui 


(4.13) 


As  mentioned  earlier  (Jsxjime  and  (J^sx^nne  vanish  due  the  symmetric  and  anti-symmetric 
nature  of  the  stress  and  displacement  fields.  (Jsijiine  and  (JASxJime  are  non-zero  only  m 
the  presence  of  crack  face  tractions.  Simplified  expressions  for  the  above  line  integrals  are 
available  in  Ref.  21.  The  total  J-integral  J*,  can  be  expressed  as 

Jxi  =  {Jsxi)domain  +  iJsxi)line  +  {J ASxt)  domain  +  {JASxi)line  (^-14) 

The  formulation  described  above  is  implemented  through  user  subroutines  in  the  finite  el¬ 
ement  code  MARC^"*.  The  program  computations  were  validated  with  typical  problems 
including  ones  with  cracks  along  bimaterial  interfaces  subjected  to  remote  and  crack  face  nor¬ 
mal  and  shear  tractions.  Many  of  the  results  were  also  checked  with  those  m  the  literature  . 
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About  5  paths  were  used  to  check  for  path  independency.  A  typical  path  used  is  shown  in 
Figure  4.2.  The  maximum  variation  of  the  total  J-integral  (eqn.  14)  was  found  to  be  less  than 
0.5%  between  first  (innermost)  and  the  last  (outermost)  paths  for  all  the  cases.  Since  the 
decomposition  method  allows  separation  of  individual  modes,  mode  I  and  II  contributions 
were  also  examined.  It  was  found  that  the  individual  modes  were  path  dependent  (variation 
of  4%  between  first  and  last  paths)  even  though  the  total  integral  was  path-independent. 
This  is  consistent  with  observations  in  Ref.  21.  The  individual  me  ies  even  though  being 
path  dependent  can  still  provide  rough  estimates  of  the  mode  I  and  TI  contributions. 


4.3  Numerical  Modeling  of  the  Interface  Using  Finite 
Element  Methods 

The  fiber-matrix  interface  is  modeled  using  a  contact-friction  formulation  and  the  bonding 
at  the  interface  is  simulated  with  high  stiffness  springs^.  The  fiber  and  the  matrix  phases 
are  modeled  using  quadrilateral  elements.  Duplicate  nodes  are  created  at  the  interface  on 
fiber  and  matrbe  sides.  The  fiber-matrbe  bonding  is  achieved  by  connecting  these  duplicate 
nodes  by  means  of  springs  with  a  stiffness  greater  than  the  elastic  modulus  of  the  fiber  by 
a  factor  of  100.  The  mutual  penetration  of  the  matrbe  and  the  fiber  elements  is  avoided 
by  defining  the  fiber  and  matrbe  phases  as  deformable  contact  surfaces.  A  typical  finite 
element  discretization  used  along  with  details  of  the  interface  element  are  shown  in  Figure 
4.3.  Refinement  was  carried  out  until  the  variations  in  displacement  fields  between  successive 
refinements  were  less  than  0.01%.  The  presence  of  a  smaU  pre-existing  crack  is  essential  for 
the  computation  of  energy  release  rates  using  the  EDI  method.  Thus  an  interface  crack  of 
size  less  than  0.5  rj  (where  rf  is  the  fiber  radius)  is  assumed  to  be  present  at  each  end  of 
the  push-out  specimen.  Numerical  investigations  have  shown  that  the  presence  of  cracks  of 
this  magnitude  does  not  affect  Gi  or  the  peak  load.  The  element  S’^e  around  the  crack  front 
was  a/20,  where  a  is  the  length  of  the  pre-existing  crack.  The  r“;_,i‘jns  used  for  cornputing 
the  J-integral  have  a  mesh  which  is  symmetric  about  the  crack  plane.  This  is  a  requirement 
for  the  decomposition  method.  The  stresses  at  the  interface  at  each  increment  are  monitored 
in  the  user  subroutine.  Initially  a  relative  motion  flag  is  set  to  0.  When  the  interface  failure 
criterion  is  satisfied  then  the  flag  is  changed  to  1  and  the  stiffness  of  the  springs  is  dropped  to 
zero,  indicating  the  onset  of  frictional  sliding.  Coulomb’s  law  is  used  to  model  the  frictional 

behavior. 

An  important  consideration  when  computing  energy  release  rates  in  the  presence  of  Coulomb 
type  of  sliding  friction  is  the  path  independency  of  the  J-integral.  Deng  ’  using  an  asymp¬ 
totic  approach  reported  that  the  strength  of  the  crack  tip  field  singularity  will  be  weaker 
than  for  a  moving  interface  crack  with  frictional  surface  contact,  as  a  result  of  which 
the  conventional  strain  energy  release  rate  G  will  be  zero  at  the  crack  tip.  This  is  a  however 
a  mathematical  phenomenon.  According  to  Deng^®’^^  G  will  have  finite  values  away  from 
the  crack  tip.  The  variation  of  G  with  distance  from  the  crack  tip  was  studied  in  a  push-out 
test  simulation  after  substantial  debonding  had  occured.  For  the  push-out  problem  with 
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the  geometry  as  described  previously  and  with  the  coefficient  of  friction  used,  even  with  a 
very  fine  mesh  near  the  crack  tip,  we  could  not  numerically  observe  any  appreciable  decrease 
in  G,  save  the  numerical  error  (less  than  0.5%)  which  was  present  even  in  the  validation 
studies.  A  possible  explanation  for  the  failure  to  observe  the  decrease  in  G  may  be  due  to 
the  fact  that  in  the  case  of  thin-slice  push-out  test  where  the  cracks  initiate  at  the  support 
end,  frictional  effects  are  not  found  to  be  appreciable.  This  is  mainly  due  to  the  relaxation 
of  residual  stresses  in  the  debonded  region  and  the  absence  of  additional  radial  compres¬ 
sive  stresses  caused  by  Poisson’s  expansion  of  the  fiber  (as  in  the  case  of  loading  end  crack 
initiation  in  CMCs).  According  to  Deng^®’^^  a  strain  energy  release  rate  type  of  fracture 
criterion  and  any  contour  around  the  crack  tip  can  be  used  to  represent  G  values  as  long 
as  one  is  consistent  with  the  path  distance  from  the  crack  tip.  In  this  paper  an  integral 
path  with  radius  Tc  =  0.075r/  is  consistently  used  in  all  the  cases.  In  the  case  of  bimaterial 
interfaces  with  Coulomb  friction  since  the  classical  strain  energy  release  rate  vanishes  at  the 
crack  tip  the  evaluated  value  of  the  J-integral  cannot  be  interpreted  as  an  intrinsic  fracture 
energy  quantity  in  precise  terms.  However,  this  value  can  be  used  for  comparison  purposes 
between  different  composite  systems.  In  all  the  cases  studied  here  (except  section  4.3)  the 
crack  length  is  extremely  small  and  the  crack  is  stationary.  Negligible  sliding  and  friction 
occurs  near  the  crack  tip;  hence  the  J-integral  values  computed  for  these  cases  was  path 
independent.  The  line  integral  contributions  are  found  to  be  insignificant  for  these  cases  and 
the  value  of  the  J-integral  is  from  the  domain  integral  contributions. 

Elastic  constitutive  behavior  is  assumed  for  the  fiber  phase.  Matrix  is  assumed  to  be  a 
rate  independent  elastic-plastic  material.  The  region  in  the  vicinity  of  the  interface  crack  is 
assumed  to  be  elastic  as  experimental  observations^®  have  shown  that  crack  tip  plasticity  is 
not  present  for  interfax^e  cracks  in  the  metal  and  ceramic  matrix  composite  systems  considered 
in  this  study.  The  temperature  dependency  of  the  elastic  and  inelastic  properties^’^®“®^  of 
the  constituent  phases  are  included  in  the  analysis.  The  material  properties  used  in  this 
study  are  given  in  Figure  4.4  and  Table  1. 


4.3.1  Simulation  of  the  Single  Fiber  Push-out  test 

The  single  fiber  push-out  test  is  analyzed  using  an  axisymmetric  cylindrical  model.  A  three 
phase  model  consisting  of  fiber,  matrix  and  a  composite  phase  with  homogenized  properties 
has  been  used.  The  push-out  test  is  modeled  following  a  sequence  of  steps  described  in 
detail  elsewhere®.  The  appropriate  boundary  conditions  for  the  steps  are  shown  in  Figure 
4.5.  Residual  thermal  stresses  are  induced  during  the  composite  consolidation  process  due 
to  the  GTE  mismatch  between  the  fiber  and  the  matrix.  A  brief  outline  of  the  steps  involved 
is  given  below: 

1.  Cooling  of  the  bulk  composite  from  the  reference  or  stress  free  temperature. 

2.  Preparing  the  push-out  specimen  by  cutting  a  thin  slice  from  the  bulk  composite.  This 
process  results  in  shear  stresses  at  the  interface  due  to  differential  axial  residual  strains 
between  the  fiber  and  matrix. 
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3.  The  push-out  load  is  applied  by  prescribing  incremental  displacements  to  the  top  face 
of  the  fiber. 


The  state  of  residual  stress  at  the  interface  can  be  modified  by  pre-straining  the  composite 
to  levels  sufficiently  high  so  as  to  introduce  plastic  strains  in  the  matrix.  In  this  case  an 
additional  step  is  introduced  in  between  1  and  2  with  the  same  boundary  conditions  as  that 
of  step  1.  The  composite  is  strained  in  the  longitudinal  (or  fiber)  direction  by  prescribing 
displacements  to  the  top  face  of  the  composite. 


4.4  Results  and  Discussion: 


The  push-out  test  is  simulated  following  the  sequence  of  steps  listed  in  the  previous  section. 
SCS-6/Ti-15-3,  SCS-6/Ti-6-4  and  SCS-6/Timetal-2lS  are  the  typical  MMC  systems  consid¬ 
ered  in  this  study.  The  thickness  of  the  push-out  specimens  in  all  the  cases  is  taken  to  be  0.5 
mm  and  the  fibre  volume  fraction  was  chosen  to  be  35%.  Both  mode  I  and  mode  II  energy 
release  rates  are  measured  for  the  push  out  problem.  Preliminary  studies  of  the  push-out 
problem  showed  that  the  mode  II  energy  release  rate  contribution  to  the  total  J-integral  was 
10“*  times  that  of  mode  I.  Thus  the  mode  I  contribution  is  negligible  and  the  problem  of 
interface  debonding  due  to  push-out  load  can  be  treated  essentially  as  a  mode  II  problem. 


4.4.1  Gi  after  Specimen  Slicing 


As  mentioned  earlier  the  specimen  preparation  stage  which  is  modeled  as  per  step  2  in  the 
simulation  sequence,  results  in  the  redistribution  of  the  residual  stresses  near  the  composite 
ends.  This  causes  stress  concentration  at  the  pre-existing  cracks  at  both  the  ends  of  tk' 
specimen.  The  calculations  of  Gir  (the  value  of  Gi  after  slicing  and  before  the  application 
of  indenter  load)  from  the  EDI  method  for  different  material  systems  are  shown  in  Table  2. 
These  are  values  obtained  from  the  total  J-integral.  It  can  be  seen  that  Gir  of  the  order  of 
40  J /tit?  are  present  at  the  specimen  ends.  The  variation  of  interfacial  shear  stress  along  the 
specimen  thickness  (Figure  4.6)  indicates  that  the  stress  state  at  both  ends  of  the  specimen 
is  symmetric  in  nature  and  so  are  the  other  components  of  stress.  This  leads  to  symmetric 
values  for  the  energy  release  rates  at  both  the  ends  of  the  specimen.  This  also  indicates 
the  tendency  of  the  fiber  to  protrude  out  of  the  matrix  at  either  end.  Gto  which  is  the 
critical  value  required  for  debonding  is  taken  to  be  higher  than  Gir  in  these  simulations,  as 
microscopic  studies  of  MMC  specimens  do  not  reveal  significant  interfacial  debonding  before 
the  application  of  the  mechanical  load.  Hence  the  values  of  Gir  listed  in  Table  2  serve  as 
a  lower  bound  for  Gic  of  these  MMC  systems.  To  illustrate  this,  step  2  is  simulated  by 
assuming  Gic  value  to  be  lower  than  40  J/m.^  in  the  case  of  SCS-6/Ti-15-3  composite.  In 
this  case  debonding  occurs  after  the  end  of  the  processing  stage  (even  before  application 
of  load  by  the  punch)  and  proceeds  with  unstable  crack  propagation  from  either  end  for 
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most  of  the  specimen  thickness.  It  should  however  be  pointed  out  that  surface  roughness  in 
interfaces  may  play  a  role  in  arresting  this  unstable  crack  propagation. 

The  processing  simulations  are  also  performed  on  some  CMCs  like  SCS-6/CAS  I  and  SCS- 
6/Borosilicate  to  compare  energy  release  rate  values  with  those  obtained  for  the  MMCs. 
These  energy  values  are  also  listed  in  Table  2.  As  expected  the  Gir  values  for  the  CMCs 
were  at  least  one  order  less  than  those  for  the  MMCs.  This  is  to  be  expected  since  the 
residual  stresses  generated  in  the  CMCs  are  much  smaller  compared  to  that  in  MMCs. 
To  examine  the  contribution  to  the  energy  release  rates  from  residua)  stresses,  numerical 
simulations  were  performed  in  which  the  MMC  specimens  after  processing  were  plastically 
strained  before  slicing.  The  results  for  0.5%  and  1.3%  total  strain  are  shown  in  Table  2. 
The  effect  of  inducing  plastic  strain  in  the  matrix  is  to  relax  the  residual  stresses.  This  is 
manifested  as  a  reduction  in  the  Gir  values  obtained  after  slicing.  We  notice  from  Table  2 
that  with  increase  in  the  amount  of  plastic  strain  the  energy  release  rates  reduce  to  levels 
which  are  comparable  with  those  obtained  for  CMCs. 


The  relative  contribution  of  the  fiber  axial  residual  stress  to  the  overall  residual  stress  state 
is  also  examined.  The  strain  energy  in  the  fiber  can  be  estimated  from  the  axial  residual 
stress  (neglecting  other  stress  components)  in  the  fiber  using  the  following  relation: 


(4.15) 


where  E  is  the  Young’s  modulus  of  the  fiber  and  r  is  its  radius.  The  strain  energy  release 
rate  contribution  from  the  fiber  can  be  computed  by  differentiating  the  above  relation  with 
respect  to  the  crack  surface  area  A.  This  is  given  by 


dU  0-2 
dA  “ 


(4.16) 


Eqn.  4.16  gives  the  contribution  of  axial  residual  stress  towards  the  strain  energy  release 
rate  Gir-  It  can  be  seen  that  the  energy  computed  from  this  equation  depends  only  on  the 
stress  level  and  the  Young’s  modulus  of  the  fiber.  To  examine  the  validity  of  this  equation, 
a  parametric  study  is  conducted  for  different  levels  of  elastic  mismatch  given  by  Em/ Ep 
which  is  related  to  the  Dundur’s  parameter  The  strain  energy  release  rates  computed 
from  the  EDI  method  are  shown  in  Table  3  along  with  the  corresponding  residual  stress 
values.  The  error  in  the  estimation  of  G<r  due  to  residual  stress  from  Eqn.  4.16  and  the 
EDI  method  is  plotted  in  Figure  4.7.  It  can  be  seen  that  the  %  error  decreases  as  a  0 
(the  case  of  Err^atrix  =  Efiber)-  TMs  can  be  attributed  to  the  crack  driving  force  being  of 
the  same  order  in  the  fiber  and  in  the  matrix  as  a  ^  0  .  This  analysis  also  implies  that 
axial  residual  stress  is  a  major  contributor  to  Gir-  For  the  composite  system  SCS-6/Ti-15-3 
[Em / Ef=Q.2Z)  under  consideration  in  this  study,  the  error  in  the  Gir  calculation  is  about 
25%  (Figure  4.7). 
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4.4.2  Debond  Initiation 


Previous  studies®  have  indicated  that  in  the  presence  of  processing  induced  residual  stresses, 
debonding  initiated  at  the  support  end  for  the  thicknesses  of  the  order  of  0.5  mm  or  lower  con¬ 
sidered  in  this  work.  This  observation  is  supported  by  experimental  push-out  test  results^^. 
In  this  paper  the  energy  release  rates  for  both  the  ends  at  the  time  of  crack  initiation  are 
examined.  The  applied  displacement  u  has  been  normalized  with  respect  to  Umax  in  the 
plots  where  Umax  is  the  maximum  applied  displacement  at  which  debond  initiation  occurs. 
The  results  for  the  loading  and  support  ends  are  shown  in  Figure  4.8.  From  this  figure  it  is 
observed  that  the  energy  release  rate  for  the  support  end  keeps  on  increasing  (from  the  value 
at  the  end  of  the  slicing  stage)  till  the  assigned  critical  value  of  the  energy  release  rate  Gic  of 
50  J/m^  is  reached  at  which  point  crack  initiation  occurs.  This  is  in  contrast  to  the  energy 
release  rate  curve  for  the  loading  side  which  shows  negligible  change  from  its  value  at  the  end 
of  the  slicing  stage.  This  phenomenon  is  a  direct  consequence  of  the  presence  of  processing 
induced  residual  stresses.  As  mentioned  in  the  previous  section,  even  before  the  application 
of  load,  the  fibers  have  a  tendency  to  protrude  from  the  ends  due  to  the  redistribution  of 
residual  stresses  after  slicing.  When  the  load  is  applied  the  axial  residual  stresses  at  the 
support  end  are  in  the  same  direction  as  the  applied  load  whereas  at  the  loading  end  they 
oppose  the  applied  load.  This  causes  the  stress  intensity  at  the  support  end  to  increase  and 
thus  debonding  initiates  from  this  end.  In  this  case  the  strain  energy  at  the  support  end 
shoots  up  and  reaches  the  critical  value  even  before  any  appreciable  change  is  seen  at  the 
loading  end. 

The  application  of  push-out  load  is  simulated  after  altering  the  residual  stress  state  by 
introducing  plastic  strains  in  the  matrix  by  subjecting  the  specimen  to  tensile  loads  (refer 
section  4.1).  The  variation  of  these  stresses  with  equivalent  plastic  strain  is  shown  in  Figure 
4.9.  Results  for  three  specimens  which  are  subjected  to  three  different  applied  strains  of  (a) 
0.9%  (b)  1.3%  and  (c)  1.5%  are  presented  in  this  section.  The  strain  energy  release  rates 
for  crack  initiation  for  these  cases  are  taken  to  be  different  (21  J/m^  for  cases  (a)  and  (b) 
and  7  J/rri^  for  case  (c)  )  from  the  zero  plastic  strain  case  discussed  earlier.  This  is  done  for 
keeping  the  loads  about  the  same  for  crack  initiation.  The  emphasis  here  is  not  to  extract 
the  critical  energy  release  rates,  but  to  study  the  variation  of  the  strain  energy  release  rates 
with  increasing  loeid.  The  variation  of  strain  energy  release  rates  at  the  loading  and  support 
ends  are  discussed  in  the  subsequent  paragraphs. 

In  the  first  case  the  specimen  is  subjected  to  0.9%  strain,  sliced  and  then  compressive 
load  is  applied  on  the  fiber.  The  variation  of  energy  release  rates  at  the  two  ends  of  the 
specimen  after  slicing  and  upto  initiation  of  debonding  are  shown  in  Figure  4.10.  The 
behavior  exhibited  in  this  case  by  the  energy  release  rates  is  similar  to  the  case  when  no 
tensile  strain  is  applied.  The  support  end  Gi  continues  to  increase  until  debonding  is 
initiated  at  this  end.  However  Gi  for  the  loading  end  decreases  in  magnitude  from  the  value 
at  the  end  of  slicing  and  reaches  very  small  values  and  maintains  these  values  for  the  rest  of 
the  initiation  stage.  The  decrease  in  Gi  at  the  loading  end  is  due  to  the  drop  in  the  stress 
intensity  at  this  end  as  the  influence  of  the  axial  residual  stresses  is  reduced. 
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In  the  second  case  the  specimen  is  subjected  to  1.3%  tensile  strain  before  slicing.  The  Gir 
value  after  slicing  is  only  of  the  order  of  4  J/m^  compared  to  the  as-received  value  of  the 
order  of  40  J/m^.  The  results  for  this  case  are  shown  in  Figure  4.11.  The  variation  of  the 
energy  release  rates  show  a  markedly  different  behavior  from  the  previous  cases.  Towards 
the  beginning,  the  energy  release  rates  for  both  ends  increase  more  or  less  similarly.  With 
increase  in  displacement  the  loading  end  energy  release  rate  attains  values  far  higher  than 
that  of  the  support  end.  Debonding  finally  initiates  at  the  loading  end.  It  has  to  be  noted 
that  in  this  case  the  ratio  of  the  energy  release  rates  at  the  two  ends  is  much  lower  compared 
to  the  earlier  cases.  This  case  is  very  near  the  transition  of  the  debond  initiation  switching 
from  the  support  to  the  loading  end. 

In  the  third  case  the  specimen  is  strained  to  1.5%  strain  uniaxially.  The  axial  residual 
compressive  stresses  in  the  fiber  are  completely  removed.  The  energy  release  rates  for  both 
ends  of  the  specimen  are  shown  in  Figure  4.12.  In  this  case  debonding  is  observed  to  initiate 
from  the  loading  side  similar  to  the  behavior  exhibited  by  CMCs^°’^®  where  residual  stresses 
are  negligible.  The  energy  release  rate  for  the  loading  side  gradually  increases  in  magnitude 
till  debonding  occurs.  The  support  end  energy  release  rates  remains  constant  and  do  not 
increase.  This  effect  is  observed  experimentally  by  Majumdar  et  al  . 

To  evaluate  the  response  of  CMCs,  the  push  out  test  is  simulated  for  the  SCS-6/Borosilicate 
system  for  the  same  thickness.  The  residual  stresses  induced  in  this  case  due  to  processing 
are  much  smaller  compared  to  the  stresses  induced  in  MMCs.  The  energy  release  rates  upto 
debond  initiation  are  shown  in  Figure  4.13.  With  application  of  load  by  the  indenter  the 
loading  end  energy  release  rates  continue  to  increase  until  debonding  initiates  from  this  end. 
The  energy  release  rates  for  the  support  end  show  negligible  change  from  the  value  at  the 
end  of  the  slicing  stage  which  leads  us  to  conclude  that  debonding  is  not  likely  to  initiate 
from  this  end.  Fiber  will  get  pushed  out  when  the  loading  end  crack  reaches  the  support 
end.  These  results  illustrates  the  considerable  role  played  by  processing  induced  residual 
stresses  on  interface  debonding  initiation. 


4.4.3  Application  of  the  EDI  Method  to  Interfacial  Crack  Prop¬ 
agation 

The  EDI  method  is  advantageous  for  crack  propagation  studies  as  non-singular  elements  can 
be  used  in  the  region  around  the  crack  tip.  This  eliminates  the  need  to  remesh  with  every 
stage  of  crack  advance.  The  crack  is  advanced  when  the  energy  release  rate  for  the  crack  tip 
Gi  exceeded  the  prescribed  Gic  of  50  J/m^.  The  behavior  exhibited  by  the  energy  release 
rates  during  the  course  of  crack  propagation  is  examined  for  the  SCS-6/Ti-15-3  system. 
After  the  processing  stage,  the  specimen  is  sliced  and  load  is  applied  with  a  punch  till  the 
fiber  is  completely  debonded.  The  sequence  of  debonding  is  shown  in  Figure  4.14.  In  this 
figure  the  crack  tip  locations  for  both  sides  are  plotted  as  a  function  of  fiber  displacement 
u  which  has  been  normalized  with  respect  to  maximum  fiber  displacement  Umax-  u-jumax 
=  0  indicates  the  end  of  the  slicing  stage  and  start  of  the  application  of  load  by  the  punch. 


52 


u/umax  —  1  indicates  the  stage  where  the  fiber  gets  pushed  out.  The  location  of  the  crack  tip 
X  is  normalized  with  respect  to  the  specimen  thickness  L.  The  curve  which  represents  the 
support  side  crack  tip  location  rises  steadUy  till  it  intersects  the  curve  which  represents  the 
loading  end  crack.  At  this  point  the  fiber  gets  pushed  out.  The  curve  which  represents  the 
loading  end  crack  maintains  zero  slope  for  all  of  the  duration  of  the  push-out  test  indicating 
that  this  crack  does  not  propagate. 

The  energy  release  rates  for  both  ends  are  measured  during  the  course  of  debonding  and  are 
shown  in  Figure  4.15.  Prom  this  figure  we  see  that  the  energy  release  rates  for  the  support  end 
rises  until  debonding  occurs.  From  the  point  of  debond  initiation  till  the  fiber  gets  pushed 
out,  the  energy  release  rate  remains  essentially  constant.  When  complete  debonding  takes 
place  and  the  fiber  gets  pushed  out,  the  energy  levels  drop  to  a  negligible  value.  The  energy 
release  rate  curve  for  the  loading  side  exhibits  a  behavior  different  from  that  of  the  support. 
In  this  case  the  energy  release  rate  shows  negligible  change  from  the  value  obtained  at  the 
end  of  the  slicing  stage  until  the  support  side  crack  tip  crosses  the  centre  of  the  specimen. 
Then  the  energy  release  rates  continue  to  decrease  until  very  near  the  stage  at  which  the 
fiber  gets  pushed  out.  This  is  because  the  effect  of  the  compressive  residual  stresses  (acting 
opposite  to  the  direction  of  the  applied  load)  decreases  which  leads  to  a  reduction  in  stress 
intensity  at  the  loading  end  crack  tip.  When  the  fiber  gets  pushed  out  the  support  end 
energy  release  rate  drops  to  a  very  low  value  as  the  resistance  to  further  sliding  is  only  due 
to  friction. 


4.5  Summary  and  Conclusions 

A  fracture  mechanics  approach  based  on  the  J-integral,  computed  using  the  Equivalent 
Domain  Integral  (EDI)  method,  is  used  to  examine  the  interfacial  debonding  process  in 
MMCs  and  CMCs  during  a  fiber  push-out  test.  Strain  energy  release  rates  for  the  fiber- 
matrix  interface  due  to  processing  induced  residual  stresses  are  determined  by  the  above 
method.  These  rates  serve  as  a  lower  bound  for  the  Mode  II  fracture  toughness  of  the 
composite  systems  examined.  The  phenomenon  of  bottom  (support  end)  debonding  observed 
in  MMC  systems  during  the  push-out  test  is  explained  by  the  study  of  the  variation  of  energy 
release  rates  for  the  loading  and  support  end  cracks  upto  debond  initiation.  The  effect  of 
residual  stresses  on  crack  initiation  is  studied  by  altering  the  residual  stress  state  by  plastic 
straining.  It  is  found  that  crack  initiation  switches  from  the  support  to  the  loading  end  with 
reduction  in  residual  stresses  with  increase  in  plastic  strain.  For  the  specimen  geometry  and 
material  system  studied,  this  switching  occurs  when  the  Gir  values  drop  to  values  of  the  order 
of  jrr?  from  the  as-received  value  of  around  jvn} .  It  is  shown  that  the  EDI  method  can 
be  used  to  simulate  interface  crack  propagation  without  the  need  for  remeshing  with,  crack 
advance.  Evaluation  of  strain  energy  release  rates  through  the  use  of  J-integral  for  cracks 
along  bimaterial  interfaces  in  the  presence  of  friction  is  still  a  “mathematical  problem”  and 
care  needs  to  be  exercised  in  the  interpretation  of  the  numerical  predictions. 
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Table:  Material  properties  for  ceramic  matrices 


Matrices 

E 

(GPa) 

G 

(GPa) 

a 

(X  10-6 /OC) 

Borosilicate 

68 

28 

3.2 

CAS  I 

88 

36 

5.0 

Material  properties  obtained  from:  Ref.  31 


Table  4.1:  Material  properties  for  ceramic  matrices 


Table:  Computed  strain  energy  release  rates  for  different  composite  systems 


Material  System 

(fiber  volume  fraction=  0.35) 

AT 

(OC) 

Fiber  Residual  Stress 
Axial  Radial 

(MPa)  (MPa) 

G;  (Residual) 
(J/m2) 

SCS-6/Ti-15-3  ❖ 

790 

-750 

-223 

37.6 

SCS-6m-6-4  * 

875 

-858 

-272 

42.3 

SCS-6/Timetal  21s  t 

875 

-985 

-304 

47.6 

SCS-6/CASI  -r 

1230 

-220 

-100 

7.86 

SCS-6/Borosilicate 

600 

20 

8 

0.19 

SCS-6/Ti-15-3  (after  straining  to  0.5%  strain) 

-582 

-248 

14.5 

SCS-6/Ti-15-3  (after  straining  to  1.3%  strain) 

-423 

-94 

3.79 

Material  properties  obtained  from: 

♦>  Ref  30 
*  Ref  1 
t  Ref  32 
4-  Ref  31 

Table  4.2:  Material  properties  for  ceramic  matrices 
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Table  :  Parametric  study  to  examine  effect  of  elastic  mismatch 


Material  System  * 

(fiber  volume  fraction  =  0.35) 

Dundur'^ 

Parameter 

a 

Fiber  Residual  Stress 
Axial  Radial 

(MPa)  (MPa) 

Gj  (Residual) 
(J/m2) 

Eni  =  0.23Ef 

0.619 

-750 

-223 

37.6 

Em  =  0.3  Ef 

0.489 

-906 

-387 

62.3 

Em  =  0.5  Ef 

0.320 

-1120 

-458 

73.9 

Em  =  0.75  Ef 

0.128 

-1350 

-497 

84.5 

Em=  Ef 

-0.015 

-1450 

-536 

94 

*  Thermal  properties  used  is  the  same  as  that  of  SCS-6/Ti-15-3  composite 
Ef  -  Young's  Modulus  of  Fiber 

-  Young's  Modulus  of  Matrix 


Table  4.3;  Material  properties  for  ceramic  matrices 


58 


(a)  Coordinate  system  and  contour  G 
around  the  crack  tip 


(b)  Domain  bounded  by  contours 
Fq  andF] 


Figure  4.1:  Contours  used  in  the  Equivalent  Domain  Integral  Method 
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Figure  4.2:  Typical  path  used  in  the  energy  release  rate  computations 
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DETAILS  OF  REGION  A 


1-1',  2-2’ 

duplicate  nodes 

(with  the  same  co-ordinates) 

at  the  interface 

1-2,  l'-2' 

Contact  edges  which  can  not 
penetrate 

n,  t 

normal  and  tangential 
directions 

Stiffness  of  springs  in  normal 
and  tangential  directions 

Figure  4.3:  Finite  element  model  and  details  of  the  interface 
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Fiber:  Silicon  Carbide  (SCS-6);  v  =  0.25; diameter  =  140  M:n 
_  400-1 - 


B  B  B 


200  400  600  800 

o 

Temperature  C 

Matrix; 

Ti-15V-3AI 
V  =  0.36 


L.  hhh  h 

hH  ^  B 

H  B  B  BB 


0  200  400  600^  800  1000 

Temperature  C 


Timetal  21s 
V  =  0.34 


Ti-6AI-4V 


V  =  0.3 


SRI 


■  Timetal  21s  ‘35 

•  Ti-6-4  5 

▲  Ti-15-3 

— I - 1 - 1 - 1 - 

200  400  600  800  1000 

■  Timetal  21s 

•  Ti-6-4  1 

A  Ti-15-3 


J  H  j  B 
Hh 


H  H  ? 

— I - \ - 1 - 1 1 

200  400  600  800  1000 


0  200  400  600  800  '  1000 


Figure  4.4:  Variation  of  material  properties  with  temperature 


iFiberm 


(a)  (b)  (c) 


Figure  4.5:  Boundary  conditions  used  in  the  Finite  Element  Analysis 
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Interfacial  Shear  Stress,  MPa 


MATRIX 

FIBER 


Figure  4.6:  Shear  stress  distributions  along  the  interface  for  specimen  thickness  of  0.5  mm 


%  Error  in  energy  release  rates  computed 
from  fiber  stress 


Figure  4.7:  Percentage  error  in  energy  release  rates  calculated  using  the  axial  fiber  residual 
stress  plotted  as  a  function  of  the  ratio  of  the  matrix  modulus  to  the  fiber  modulus 


Figure  4.8:  Variation  of  the  energy  release  rates  from  the  end  of  the  processing  stage  till 
debond  initiation  as  a  result  of  indenter  load  application 
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Fiber  Residual  Stress  (absolute,  MPa) 


Axial 

Radial 


Figure  4.9:  Variation  of  axial  and  radiaJ  components  of  residual  stress  in  the  fiber  with 
equivalent  plastic  strain  induced  in  the  matrix  as  a  result  of  straining 
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Figure  4.10:  Variation  of  energy  release  rates  till  debond  initiation  upon  application  of  load 
by  the  indenter  after  straining  to  0.9  %  total  strain 
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Figure  4.11:  Variation  of  energy  release  rates  till  debond  initiation  upon  application  of  load 
by  the  indenter  after  straining  to  1.3  %  total  strain 
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Figure  4.12:  Variation  of  energy  release  rates  till  debond  initiation  upon  application  of  load 
by  the  indenter  after  straining  to  1.5  %  total  strain 
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Figure  4.13:  Variation  of  energy  release  rates  till  debond  initiation  upon  application  of  load 
by  the  indenter  for  ceramic  matrix  composites 
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Material: 
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Propagating 

Crack  Crack  Tip 


Figure  4.14:  Crack  tip  locations  for  the  loading  and  support  ends  from  the  end  of  processing 
stage  till  the  fiber  gets  pushed  out 
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Figure  4.15;  Variation  of  energy  release  rates  from  the  end  of  processing  stage  till  the  fiber 
gets  pushed  out 


Chapter  5 


Evaluation  of  Fracture  Toughness  of 
MMC  Interfaces  Using  Thin-slice 
Push-out  Tests 


5.1  Introduction 

Metal  matrix  composites  (MMC)  are  being  considered  as  potential  material  systems  of  choice 
for  the  advanced  propulsion  systems  of  twenty-first  century  aircraft,  owing  to  their  high  spe¬ 
cific  stiffness  and  strength.  The  role  of  fiber-matrix  interface  in  MMCs  is  critical  to  the 
performance  of  these  composites.  Push-out  test  is  emerging  as  an  important  experimental 
tool  for  characterizing  the  interface  behavior  of  MMCs  mainly  after  its  widespread  use  in  the 
study  of  interfaces  in  ceramic  matrix  composites  (CMC).  The  presence  of  high  levels  of  resid¬ 
ual  stresses  in  MMCs  alters  the  mechanics  of  the  failure  process,  making  the  interpretation 
of  test  results  and  evaluation  of  interfacial  properties  very  difficult.  Moreover  the  relatively 
high  bond  strength  of  MMCs  and  residual  clamping  stresses  at  the  interface  limit  the  speci¬ 
men  thickness  for  push-out  experiments.  In  the  thin  slice  push-out  tests  it  is  observed  that 
[1]  interface  failure  initiates  from  the  bottom  of  the  specimen  and  not  from  the  loading  end  as 
in  the  case  of  CMCs  [2].  Chandra  and  co-workers  [2,3,4]  analyzed  the  thin-slice  push-out  test 
and  simulated  the  initiation  of  interface  failure  and  its  propagation  using  non-linear  finite 
element  methods.  This  work  was  able  to  explain  many  of  the  experimental  observations  [5] 
of  push-out  behavior  at  room  and  elevated  temperatures.  They  also  proposed  a  method  for 
evaluating  the  interfacial  shear  strength  from  the  push-out  data.  However  fracture  tough¬ 
ness  (Gic)  is  increasingly  being  considered  as  a  more  appropriate  quantitative  measure  of 
the  bond  at  the  interface.  The  interfacial  fracture  energies  of  metal  matrbc  composites  are 
not  currently  available  in  the  literature  because  of  the  difficulty  in  evaluating  these  from  the 
experimental  results. 

In  this  paper,  a  numerical  method  based  on  the  equivalent  domain  integral  technique  de¬ 
veloped  previously  is  used  to  simulate  the  fracture  process  during  a  thin  slice  push-out  test 
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and  determine  the  critical  interfacial  fractme  toughness  Gic-  Experimental  correlations  are 
made  with  push-out  test  results  of  ceramic  matrix  composites  to  validate  the  Gic  predictions 
using  the  EDI  method.  This  approach  is  then  used  to  evaluate  the  Gic  of  some  as  processed 
SCS-6/Titanium  matrix  composites  from  the  experimental  data. 


5.2  Numerical  Simulation  of  the  Push-out  Test  Using 
Finite  Element  Method 


The  thin  slice  push  out  test  involves  pushing  a  fiber  out  of  a  composite  specimen  by  the 
application  of  compressive  load  through  a  flat  indenter.  Modeling  such  a  test  involves  three 
major  stages.  Step  1  involves  modeling  the  composite  consolidation  at  high  temperatures 
and  then  cooling  down  to  room  temperature  which  induces  residual  stresses  in  the  composite. 
Step  2  is  the  specimen  preparation  stage  consisting  of  slicing  a  thin  push-out  specimen  fi'om 
the  bulk  composite  which  results  in  stress  redistribution  at  the  two  ends  of  a  fiber.  Step  3  is 
modeling  the  loading  stage  where  compressive  load  is  applied  to  the  fiber  until  it  is  pushed 
out  of  the  specimen.  A  detailed  description  of  the  above  sequence  of  steps  is  provided 
elsewhere  [2].  The  test  is  analyzed  using  an  axisymmetric  cylindrical  model.  A  three  phase 
model  consisting  of  fiber,  matrix  and  a  composite  phase  with  homogeneous  properties  has 
been  used.  The  interface  is  modeled  as  a  contact  surface  and  the  bonding  at  the  interface  is 
simulated  using  high  stiffness  springs.  The  interface  failure  process  is  modeled  using  a  strain 
energy  based  failure  criteria  given  by: 

Gi  >  Gic  (5-1) 

where  Gi  is  the  strain  energy  release  rate  for  the  interface  crack  and  Gic  is  the  critical  value  of 
the  strain  energy  release  rate.  The  interface  crack  propagates  further  when  its  stra,in  energy 
release  rate  exceeds  the  critical  value.  The  total  strain  energy  release  rate  of  the  propagating 
interface  crack  is  computed  using  the  Ekjuivalent  Domain  Integral  (EDI) [6,7]  method.  The 
EDI  method  was  implemented  in  the  commercial  finite  element  analysis  softwme  MARC  [8] 
through  user  subroutines.  A  schematic  of  the  finite  element  model  is  shown  in  Figure  5.1. 
A  small  interface  crack  of  the  order  of  0.5r/  (where  rf  is  the  fiber  radius)  is  introduced  at 
both  the  ends  of  the  specimen  during  step  2  of  the  simulation  for  the  purpose  of  computation 
of  energy  release  rates.  A  path  of  radius  Tc  =  0.075r/  from  the  crack  tip  is  consistently  used 
in  all  the  cases  to  compute  Gi  as  theoretical  studies  [9]  have  shown  that  the  strain  energy 
release  rate  is  mathematically  zero  at  the  crack  tip  for  a  crack  at  a  bimaterial  interface  in 
the  presence  of  Coulomb  friction. 


5.3  Results  and  Discussion 

The  numerical  method  described  in  the  previous  section  is  applied  to  different  composite 
systems.  The  material  properties  and  the  reference  temperatures  are  taken  from  available 
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literature  [10,11,12,13].  The  levels  of  the  axial  residual  stress  in  the  fiber  predicted  by  the 
numerical  simulation  for  SCS-6/Ti-15-3  (fiber  volume  fraction  =  0.35)  shown  in  Table  1  is 
of  the  order  of  750  MPa.  This  compares  reasonably  well  with  the  values  of  749±47  MPa 
measured  experimentally  by  Pickard  et  al.[14]  using  the  dissolution  technique.  The  values 
obtained  for  other  composite  systems  are  also  shown  in  Table  1. 

Slicing  of  a  push-out  specimen  from  a  bulk  composite  (step  2)  results  in  the  redistribution 
of  the  residual  stresses  near  the  composite  ends.  This  causes  stress  concentration  at  the 
pre-existing  cracks  at  both  the  ends  of  the  specimen.  The  calculations  of  Gir  (the  value  of 
Gi  after  slicing  and  before  the  application  of  indenter  load;  purely  due  to  residual  stresses; 
from  the  EDI  method  for  different  material  systems  are  also  shown  in  Table  1.  It  can  be  seen 
that  a  Gir  of  37.6  J/m^  is  present  at  the  specimen  ends  for  SCS-6/Ti-15-3.  The  stress  state 
at  both  ends  of  the  specimen  is  symmetric  in  nature.  This  leads  to  symmetric  values  for  the 
energy  release  rates  at  both  the  ends  of  the  specimen.  This  also  indicates  the  tendency  of 
the  fiber  to  protrude  out  of  the  matrbi;  at  either  end.  Gic,  which  is  the  critical  value  required 
for  debonding  is  taken  to  be  higher  than  Gir  in  these  simulations,  as  microscopic  studies  of 
MMC  specimens  do  not  reveal  significant  interfacial  debonding  before  the  application  of  the 
mechanical  load.  Hence  the  values  of  Gir  listed  in  Table  1  serve  as  a  lower  bound  for  Gic  of 
these  MMC  systems.  The  processing  simulations  are  also  performed  on  typical  CMCs  such 
as  SCS-6/CAS  I  and  SCS-6/Borosilicate  to  compare  energy  release  rate  values  with  those 
obtained  for  the  MMCs.  As  expected  the  Gir  values  for  the  CMCs  were  at  least  one  order 
less  than  those  for  the  MMCs.  This  is  to  be  expected  since  the  residual  stresses  generated 
in  the  CMCs  are  much  smaller  compared  to  that  in  MMCs. 

Previous  studies  [4]  have  indicated  that  in  the  presence  of  processing  induced  residual  stresses 
debonding  initiated  at  the  support  end  for  MMCs  for  the  thicknesses  of  the  order  of  0.5  mm. 
Profilograph  measurements  are  conducted  on  interrupted  push-out  test  specimens  before  and 
after  the  tests  and  are  shown  in  Figure  5.2.  The  loaded  fiber  (Fiber  2  in  Figure  2)  does  not 
show  any  change  in  height  after  the  tests  on  the  loading  side  whereas  the  measurements  on 
the  support  end  indicate  that  the  loaded  fiber  pops  significantly  out  of  the  matrix.  Positions 
of  the  fiber  at  both  ends  of  the  specimen  before  and  after  application  of  the  load  from 
the  numerical  simulation  is  shown  in  Figure  5.3.  Prom  this  figure  it  can  be  seen  that  the 
numerical  simulations  show  a  behavior  similar  to  that  observed  in  the  experiments.  The 
shear  stress  distribution  along  the  interface  obtained  from  the  simulation  is  shown  in  Figure 
5.4  when  the  debonded  length  has  reached  0.6t  where  t  is  the  thickness  of  the  specimen. 
The  crack  tip  is  characterized  by  the  high  shear  stress  levels.  From  this  figure  it  can  be 
clearly  observed  that  the  debonding  is  proceeding  from  the  support  end.  A  typical  load 
displacement  curve  obtained  from  the  simulation  and  variation  of  energy  release  rates  for 
the  loading  and  support  end  cracks  is  shown  in  Figure  5.5.  From  the  energy  variation  curve 
it  is  seen  that  the  energy  release  rates  for  the  support  end  increase  from  the  value  at  the  end 
of  the  slicing  stage  till  the  critical  value  when  debonding  initiates.  For  the  period  in  which 
steady  propagation  takes  place,  the  support  end  energy  release  rates  maintain  a  constant 
level  close  to  the  critical  value.  Unsteady  propagation  is  characterized  by  the  energy  release 
rates  shooting  up  and  then  falling  as  the  fiber  is  about  to  be  pushed  out.  The  energy  release 
rates  for  the  loading  end  maintains  a  value  which  is  much  smaller  than  the  critical  value  for 
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most  of  the  duration  of  debonding  and  rises  up  to  meet  the  support  end  value  just  as  the 
fiber  is  about  to  be  pushed  out. 

The  debonding  process  in  MMCs  is  significantly  affected  by  the  presence  of  processing  in¬ 
duced  residual  stresses  and  results  in  debonding  initiating  from  the  support  end.  The  state 
of  residual  stress  at  the  interface  can  be  modified  by  prestraining  the  composite  to  levels 
sufficiently  high  so  as  to  introduce  plastic  strains  in  the  matrix.  Push-out  experiments  con¬ 
ducted  by  Majumdar  et.  al  [15]  on  SCS-6/Ti-15-3  specimens  which  were  prestrained  to 
0.9%  strain  in  tension  prior  to  push-out  showed  that  debonding  initiation  switched  from  the 
support  end  to  the  loading  end.  This  test  is  also  simulated  numerically.  The  shear  stress 
distribution  along  the  interface  from  the  simulation  is  shown  in  Figure  5.6  when  the  speci¬ 
men  is  loaded  to  28  N  which  is  considerably  less  than  the  push-out  load.  It  can  be  seen  from 
this  figure  that  the  debonding  is  propagating  from  the  loading  end  which  is  similar  to  the 
behavior  observed  experimentally.  It  was  also  seen  experimentally  that  the  push-out  load 
increased  considerably  when  the  specimen  is  prestrained.  From  the  numerical  simulations 
it  was  observed  that  for  SCS-6/Ti-15-3  specimen  which  was  prestrained  debonding  did  not 
even  initiate  for  the  load  of  «  20  N  (Table  2)  which  corresponds  to  the  peak  load  for  as 
received  specimens. 

As  mentioned  previously  interfacial  fracture  properties  for  MMCs  are  not  currently  available 
in  the  literature  as  there  are  no  established  procedures  for  evaluating  them  from  test  data. 
These  are  essential  to  fully  characterize  interface  behavior.  The  numerical  simulations  are 
repeated  at  varying  values  of  Gi-  The  peak  loads  obtained  from  the  simulations  are  correlated 
with  experimental  data  [16,17,18]  to  arrive  at  the  critical  strain  energy  release  rate  of  the 
interface  Gic  which  is  a  material  property.  In  order  to  validate  the  numerical  predictions  the 
push-out  test  was  simulated  for  Silicon  Carbide/Reaction  bonded  Silicon  carbide  ceramic 
matrix  composite  as  the  residual  stresses  are  minimum  in  CMCs  and  the  fracture  energies 
are  already  reported  [19].  The  predicted  value  of  Gic  using  our  methodology  matched  closely 
with  the  value  of  2.1  J/m^  reported  for  the  experimental  peak  load  as  shown  in. Table  2. 
The  above  methodology  is  applied  to  different  MMC  systems  (Table  2).  As  expected  the  Gic 
values  of  the  MMCs  are  considerably  higher  than  the  CMC  systems  considered  in  this  study. 
It  can  also  be  seen  that  the  predicted  fracture  toughnesses  are  about  10  J/rn^  higher  than 
the  values  of  Gir  the  energy  obtained  after  slicing  for  the  MMCs.  This  indicates  that  most 
of  the  energy  contribution  to  the  failure  process  comes  from  the  residual  stresses.  Also  the 
interfacial  fracture  energies  of  the  three  systems  are  not  far  apart  being  in  the  range  of  50-60 
J/rri^.  These  values  are  on  the  lower  side  as  they  are  based  on  the  minimum  energies  required 
for  fiber  push-out  assuming  that  there  are  no  pre-existing  flaws  such  as  matrix  cracks  in  the 
vicinity  of  the  interface.  The  interface  model  used  in  this  work  does  not  take  into  account  the 
mechanical  locking  effect  due  to  asperities  at  the  interface  which  may  increase  the  apparent 
Gic  values.  Efforts  are  underway  to  incorporate  this  effect. 
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5.4  Summary  and  Conclusions 


A  fracture  mechanics  based  approach  is  used  to  compute  residual  and  interfacial  fracture 
energies  of  some  MMC  and  CMC  systems  during  push-out  tests.  The  phenomenon  of  bottom 
debonding  in  MMC  thin  slice  specimens  is  examined  using  experimental  results  and  strain 
energy  computations.  Residual  stress  predictions  are  correlated  with  experimentally  reported 
values.  The  main  contribution  to  the  interfacial  fracture  toughness  in  MMCs  is  found  to  be 
from  the  residual  stresses  and  thus  the  values  of  Gir  can  be  considered  as  a  lower  bound  for 
the  critical  strain  energy  release  rate  of  the  interface  Gic-  This  factor  is  also  demonstrated 
in  the  push-out  test  simulations  for  pre-strained  specimens,  where  in  the  absence  of  residual 
stresses,  a  far  higher  push-out  load  is  necessary  for  the  same  critical  strain  energy  release  rate; 
consequently  no  fiber  push-out  is  observed  at  lower  loads.  Interface  fracture  toughnesses  for 
certain  typical  MMC  and  CMC  systems  are  predicted  using  the  developed  methodology. 
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Figure  5.1:  finite  element  model 
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Profile  (A)  Profile  (A) 


fiber  2  to  37  N 


(c)  Bottom  face  of  as  received  specimen  before  loading  (d)  Bottom  face  of  as  received  specimen  after  loading 

fiber  2  to  37  N 

Figure  5.2:  profilometry-as  received  specimen 
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Fiber  Position  (microns) 


Figure  5.3:  profilometry-as  received  specimen 
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Figure  5.4:  shear  stress  for  residual  stress  case 
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Figure  3.  Position  of  center  of  fiber  for 
loading  and  support  ends 


Indenter  Displacement  (microns) 


Figiu-e  4.Shear  stress  distribution  along  the  interface 
when  residual  stresses  are  present 


Indenter  Displacement  (microns) 
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Figure  5.5:  profilometry-as  received  specimen 


Figure  5.6:  shear  stress  for  no  residual  stress  case 
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Chapter  6 

Effect  of  Fiber  Fracture  and 
Interfacial  Debonding  on  the 
Evolution  of  Damage  in  Metal 
Matrix  Composites 


6.1  Introduction 

The  failure  of  fiber-reinforced  metal  matrix  composites  is  a  complex  process  involving  var¬ 
ious  mechanisms  such  as  fiber  fracture,  interfacial  debonding  and  matrix  plasticity.  The 
interaction  of  these  mechanisms  makes  the  damage  in  composites  a  highly  non-linear  pro¬ 
cess.  While  progress  has  been  made  in  understanding  and  quantifying  damage  in  composite 
mater'  als,  there  is  no  estabhshed  procedure  base  on  how  to  assess  and  characterize  damage 
[!]• 

Voyiadjis  [2]  studied  the  initiation  and  evolution  of  damage  in  a  composite  lamina  using 
a  damage  tensor  to  describe  the  initiation,  growth  and  coalescence  of  voids,  cracks  and 
debonding.  The  damage  tensor  on  this  case  was  developed  based  on  laminate  experimental 
data  and  not  from  the  micromechanical  level.  Ladeveze  [3]  focused  on  modeling  damage 
to  include  micromechanical  effects  at  a  macroscopic  level.  Continuum  damage  mechanics 
principles  were  used  to  describe  matrix  micro-cracking  and  fiber-matrix  debonding.  The 
longitudinal,  transverse  and  shear  behavior  of  a  carbon-epoxy  lamina  was  arrived  at  experi¬ 
mentally  and  applied  to  laminates.  However,  fiber  fracture  and  debonding  were  not  modeled 
explicitly.  Allen  investigated  the  issue  of  delamination  of  composite  laminates  and  treated 
the  problem  using  a  stochastic  approach.  Robertson  and  Mall  [4]  investigated  the  problem 
of  micromechanical  damage  using  the  ‘Method  of  Cells’  approach  by  Aboudi  [5]  to  predict 
the  mechanical  response  of  metal  matrix  composite  (MMC)  laminates  under  fatigue  loading. 

Recent  work  by  Lacy  et  al.  [6]  have  emphasized  the  need  for  accounting  for  the  interaction 
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of  damage  entities.  They  showed  that  the  distribution  and  interaction  of  damage  entities 
at  the  micro-scale  (within  a  Representative  Volume  Element)  can  affect  the  evolution  of 
the  damage  (and  hence  the  final  failure)  significantly,  while  having  only  a  minor  effect  on 
the  effective  moduli.  Hence,  the  size,  distribution  and  interaction  of  the  defects  need  to  be 
explicitly  taken  into  account  in  the  damage  evolution  studies. 


6.2  Damage  Mechanics 


A  general  formulation  for  the  thermomechanical  response  of  a  composite  body  undergoing 
damage  can  be  derived  from  [7]  the  Helmholtz  free  energy  potential  tp: 

(6.1) 
(6.2) 


Cij  =  P 


dtp 
defj 


where 


'ipp  -  plastic  part  of  the  potential  which  is  not  dependent  on  the  damage 
p  -  Density 

Cljki  -  Elastic  tensor  modified  by  damage 

C*jki  =  {lijrs  -  Dijrs)  Cyski  where  AjTs  is  the  damage  tensor 

Cijki  -  Elastic  tensor  without  damage 


in  equation  6.2  is  the  elastic  strain  which  is  defined  as  -  eh,  wt-ere  eh  is  the  total 
strain  and  eh  is  the  inelastic  strain  due  to  plasticity  and  creep.  For  the  jpecial  case  of  2-D 
transversely  isotropic  material  which  is  the  case  of  the  composite  lamina,  or  the  stiffness 
matrix  ^ijkl  which  is  its  inverse  reduces  to  the  following: 


S^ki- 


Jl2 

q* 

^22 
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0 
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-Dll,  -D22,  Di2,  D>2i  and  -Dee  are  the  damage  parameters.  As  S12  =  S21,  D12  and  £>21  are 
related.  This  reduces  the  number  of  independent  damage  parameters  for  the  2-D  transversely 
isotropic  case  to  four.  These  parameters  can  be  determined  from  longitudinal  and  transverse 
tension  tests  as  well  as  shear  tests. 

Currently,  experimental  methods  are  being  used  to  investigate  the  role  of  local  failure  mech¬ 
anisms  on  the  macroscopic  behavior  of  composites.  Local  damage  in  a  fiber-reinforced  com¬ 
posite  can  be  affected  by  several  geometric  and  material  parameters,  e.g.  fiber  volume 
fraction,  fiber  arrangement,  n^echanical  and  statistical  properties  of  the  fiber,  matrix  con¬ 
stitutive  behavior,  interfacial  strength  and  any  local  perturbations  of  the  above  parameters. 
This  makes  experimental  quantification  of  damage  under  different  conditions  (geometrical, 
material  and  loading)  difficult.  In  this  study  computational  techniques  are  used  to  simulate 
the  micromechanical  tests  for  extracting  the  damage  parameters.  This  enables  a  parametric 
study  of  the  influence  of  raicromechanical  parameters  such  as  fiber,  matrix  and  mterfacial 
properties,  fiber  arrangement  and  volume  fraction  on  the  damage  evolution. 


6.3  Monte  Carlo  simulation 

As  mentioned  earlier,  fracture  of  fiber-reinforced  metal  matrix  composites  (MMCs)  is  a  com¬ 
plex  process  involving  several  mechanisms  that  cause  deterioration  of  composite  properties 
such  as  fiber  fracture,  interfacial  debonding  and  matrix  plasticity.  Among  the  above  mech¬ 
anisms,  matrix  plasticity  and  interfacial  debonding  are  essentially  deterministic,  but  the 
fracture  of  brittle  fibers  is  stochastic  in  nature.  Fiber  failure  occurs  over  a  range  of  stress 
values  determined  by  the  random  distribution  of  microscopic  flaws  in  the  fiber  (Curtin  [2,  3]). 

A  simulation  of  damage  in  a  composite  laminate  can  account  for  the  stochastic  variations  of 
the  strengths  of  the  fiber  using  a  Monte  Carlo  simulation  [10].  This  method  has  been  used  to 
analyze  a  wide  range  of  physical  processes  [11].  In  general,  it  involves  a  computer  simulation 
using  random  numbers  to  solve  stochastic  problems  [12].  Such  simulations  may  be  regarded 
as  ‘numerical  experiments’.  Voleti  et  al.  [13]  applied  Monte  Carlo  simulation  procedure  to 
model  the  fiber  failure  process  in  single  fiber  composite  test.  The  stochastic  variation  of  the 
fiber  strength  along  its  length  is  modeled  using  a  ‘two-parameter  Weibull  model  [14],  which 
is  described  in  detail  in  the  next  subsection. 


6.3.1  Statistical  Model  for  Fiber  Strength  and  its  Distribution 

A  statistical  method  commonly  used  to  determine  the  strength  of  a  brittle  fiber  is  the  Weibull 
model  [14]-  In  this  model,  it  is  assumed  that  the  fiber  material  is  isotropic  and  statistically 
homogeneous  [15].  This  method  has  been  used  to  model  fiber- fragmentation  test  or  single 
fiber  composite  test,  which  has  been  used  in  the  study  of  fiber-matrix  interface  properties 
[13,  11,  12]  of  polymer  matrix  composites. 
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Coleman  [18]  examined  the  strength  of  long  fibers  taken  from  a  common  source  (e.g.  a  spool 
of  fiber).  He  showed  that  the  cumulative  strength  distribution  function,  has  the  form 

of  a  Weibull  distribution  (a/  is  the  failure  stress  of  a  fiber).  The  cumulative  probability  of 
failure  for  a  fiber  of  length  ratio  L  is  given  by: 

Pf{af)  =  1  -  (6-4) 

where  P/  is  the  probability  of  failure  of  a  fiber  at  a  stress  level  equal  to  or  less  than  aj, 
(70  is  the  Weibull  scale  parameter  for  the  unit  fiber  length  ratioy-'m’  is  the  Weibull  shape 
parameter  and  ‘L’  is  the  length  ratio  with  respect  to  a  reference  length  (the  fiber  length  at 
which  the  cq  and  ‘m’  are  determined).  The  Weibull  scale  and  shape  parameters  for  a  fiber 
material  are  determined  from  several  experiments  [1]. 


The  above  model  is  adapted  to  perform  the  Monte  Carlo  simulation  of  a  single-fiber-composite 
test  using  the  following  approach.  In  this  method,  the  length  of  the  given  fiber  is  divided 
into  a  fixed  number  of  equally-sized  ‘links’  [20].  The  strength  of  a  link  (<t/)  of  any  length  {1) 
is  given  by  the  following  equation: 


<Ti  =  (To 


1/m 


(6.5) 


where,  Lo  is  the  gage  length  at  which  the  Weibull  scale  parameter  cto  is  determined  exper¬ 
imentally.  In  a  fiber,  the  strength  of  the  links  has  a  random  distribution.  This  strength 
distribution  is  assigned  using  the  equation  below: 


<Tr(i)  =  <Ti 


In 


( 


1 

1  -  r{i) 


) 


1/m 


(6.6) 


where  ar{i)  is  the  strength  of  the  link,  i  is  the  number  of  the  fiber  link  and  r(i)  is  a 
random  number  between  zero  and  one.  ‘m’  is  the  Weibull  shape  factor.  This  equation 
assigns  a  strength  value  to  each  link  in  the  fiber,  based  on  the  value  of  the  random  number 
used  for  that  link. 


6.3.2  Fiber-matrix  interface 


Interface  (or  interphase)  is  the  zone  of  transition  separating  the  matrix  from  the  fiber  with 
thermal,  physical  and  mechanical  properties  that  can  be  different  from  that  of  fiber  and 
matrix.  The  role  of  interface  is  very  critical  to  the  performance  of  the  composites.  Recent 
experimental  investigations  [21,  1,  23]  give  clear  evidence  of  premature  interfacial  debonding 
in  metal  matrix  composites.  This  emphasizes  the  need  to  account  for  the  interfacial  behavior 
of  this  class  of  composites.  Chandra  and  Ananth  [7]  modeled  the  interface  using  spring  layer 
model.  In  this  model,  the  interface  is  considered  as  a  region  of  negligible  thickness  possessing 
the  required  stiffness  and  strength.  Debonding  is  postulated  to  occur  under  the  combined 
action  of  normal  tensile  stress  which  could  cause  mode  I  failure  and  shear  stress  which  could 
cause  mode  II  failure  at  the  interface.  A  quadratic  stress  based  failure  theory 
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is  applied  [25,  26]  where  4  is  the  interfacial  strength  in  the  normal  direction  resisting  crack 
opening,  and  r-f  is  the  shear  strength  of  the  interface.  In  the  above  equation  the  effect  of 
a(  is  included  only  when  it  is  positive.  Subsequent  to  the  debonding  process,  interfacial 
frictional  sliding  occurs  in  the  mating  debonded  surfaces  according  to  Coulomb’s  law. 


6.4  Computational  Models 


Nonlinear  finite  element  analysis  is  used  to  capture  local  damage  in  composite  lamina.  The 
commercial  Finite  Element  Analysis  software  MARC  [6]  is  used  for  this  purpose.  External 
subroutines  are  linked  to  the  main  program  to  model  phenomena  such  as  fiber  failure  and 
interfax^ial  debonding. 


6.4.1  Numerical  Implementation  of  the  Weibull  Model 

The  two-parameter  Weibull  model  described  in  section  on  statistical  model  for  fiber  strength 
is  used  for  representing  the  experimentally  observed  strength  distribution  of  the  silicon  car¬ 
bide  fiber.  The  fiber  used  in  the  study  is  considered  as  being  made  up  of  a  number  of  links 
or  divisions.  Lo  in  equation  6.5  is  taken  as  25  mm  and  I  =  0.3  mm.  The  values  for  the 
Weibull  parameters  used  in  this  case  are  shown  in  Table  6.1. 


6.4.2  Numerical  Implementation  of  the  Monte  Carlo  Method  for 
modeling  fiber  breakage 

The  Monte  Carlo  simulation  procedure  to  model  th  j  process  of  fiber  fracture  is  shown  in 
the  flow  chart  in  Figure  6.1.  The  finite  element  model  (described  in  detail  in  the  following 
section)  comprises  of  five  fibers  embedded  in  a  matrix.  Each  fiber  is  divided  into  an  arbitrary 
number  of  links.  Initially,  a  specific  strength  value,  based  on  equation  6.6,  is  assigned  to 
each  of  the  links  in  all  the  fibers.  Next,  the  finite  element  model  is  subjected  to  a  series 
of  increasing  displacements.  For  each  load  step,  the  stresses  in  the  fiber  and  matrix  are 
calculated.  The  MARC  subroutine  ELEVAR,  linked  to  this  program,  extracts  the  fiber 
stresses  from  the  solution  and  these  values  are  compared  with  the  strength  of  each  fiber 
link  through  the  use  of  another  user-subroutine,  HYPELA.  If  the  fiber  stress  exceeds  the 
assigned  link  strength,  the  link  breaks.  This  is  simulated  by  dropping  the  elastic  modulus  of 
the  fiber  link  (element)  to  a  very  low  value.  After  each  load  increment,  dummy  increments 
are  applied  to  capture  new  breaks  which  may  occur  due  to  a  re-distribution  of  stresses.  The 
incremental  loading  process  is  continued  till  a  steep  drop  in  the  stress-strain  curve  occurs  (as 
will  be  shown  in  the  results)  indicating  the  failure  of  the  lamina.  At  this  point,  the  lamina 
has  little  resistance  to  applied  load. 


84 


6.4.3  Niim©rical  Modeling  of  interfacial  debonding 

The  spring-layer  approach  described  ifi  the  section  on  fiber-matrix  interface  is  used  to  study 
interfacial  debonding  in  0°  and  90°  models.  Duplicate  nodes  are  created  at  the  interface  on 
fiber  and  matrix  sides.  The  fiber-matrix  bonding  is  achieved  by  connecting  these  duplicate 
nodes  by  means  of  high  stiffness  springs  with  a  stiffness  greater  than  the  elastic  modulus  of 
the  fiber  by  a  factor  of  100.  The  penetration  of  matrix  elements  into  the  fiber  elements  is 
avoided  by  defining  the  fiber  and  matrix  phases  as  deformable  contact  surfaces.  A  typical 
finite  element  discretization  used  along  with  details  of  the  interface  element  axe  shown  in  Fig. 
6.2.  The  stresses  at  the  interface  are  monitored  in  this  subroutine.  Initially  a  relative  motion 
flag  is  set  to  0.  When  the  interface  failure  criterion  is  satisfied  then  the  flag  is  changed  to  1 
and  the  stiffness  of  the  springs  is  dropped  to  zero  resulting  in  the  separation  of  the  interface 
nodes.  The  frictional  stresses  are  neglected  in  this  study. 

In  the  case  of  45°  lamina  model  the  spring-layer  model  is  not  used  because  once  the  fiber 
debonds  completely  there  will  be  a  total  separation  between  the  fiber  and  matrbc,  thus 
splitting  the  composite  lamina.  This  does  not  correctly  capture  the  actual  case  where  the 
matrix  continues  to  take  the  load  even  after  debonding.  An  alternative  approach  known  as 
‘matrix  layer’  approach  is  used  to  simulate  the  stiffness  reduction  due  to  interface  m  the  45 
lamina.  In  this  model,  fiber  elements  in  the  debonded  region  are  assigned  matrbc  properties. 
The  failure  criterion  used  is  the  same  as  in  the  spring  layer  approach. 


6.4.4  Finite  Element  Models  and  Boundary  Conditions 


Longitudinal  Model:  .  n.  *  i 

A  schematic  of  the  0°  lamina  is  shown  in  Figure  6.3.  The  lamina  contains  5  fibers  (the  central 

.:Der  has  two  neighbors  on  either  side)  with  the  inter-fiber  spacing  determined  by  the  volume 
faction  (16%  or  35%).  The  length  of  the  lamina  is  18  mm  which  is  found  sufficiently  longer 
than  the  load  transfer  length  of  this  fiber-matrbc  system.  A  4-node,  plane-stress  element  was 
used  for  the  model.  The  fiber  is  elastic,  and  the  random  strength  distribution  of  the  links 
is  determined  by  the  Weibull  equation.  A  length  of  about  0.1  L  (where  L  is  the  length  of 
the  composite)  is  made  infinitely  strong  on  both  the  ends  to  minimize  the  end-effects.  The 
boundary  conditions  and  loads  applied  axe  shown  in  Figure  6.3.  The  nodes  along  the  top 
and  the  lower  edges  of  the  lamina  axe  constrained  as  shown,  to  simulate  the  representative 
volume  element  or  a  repeating  unit  cell. 

The  Monte  Carlo  simulation  model  described  above  can  not  predict  the  local  stresses  around 
in  the  presence  of  a  fiber  break  accurately.  To  examine  the  effect  of  stress  redistribution  due 
to  a  broken  fiber,  a  finite  element  model  representing  three  fibers  with  the  broken  fiber  in 
the  middle  is  used.  Due  to  the  symmetry  only  one  intact  fiber  is  considered  adjacent  to  the 
broken  fiber  as  shown  in  Fig.  6.4.  The  boundary  conditions  and  applied  loads  axe  also  shown 
in  the  same  figure.  This  model  helps  in  understanding  the  cause  of  preferential  fractures  (if 
any)  in  the  Monte  Carlo  simulation  model  explained  in  the  previous  paragraph. 
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Transverse  Model 

A  schematic  of  the  transverse  (90°)  loading  model  is  shown  in  Figure  6.5.  Since  fiber  frac¬ 
ture  does  not  occur  at  transverse  strains  of  about  3%  or  less,  a  simple  finite  element  model 
capable  of  simulating  interfacial  debonding  and  matrix  plasticity  was  used.  The  boundary 
conditions  are  seen  in  the  schematic.  The  boundary  conditions  are  also  indicated  in  Fig¬ 
ure  6.5.  Debonding  is  modeled  using  the  spring-layer  model  through  the  user-subroutine 
USPRNG. 

Shear  Model; 

The  shear  behavior  of  the  lamina  is  studied  using  a  45°  ply  model.  A  schematic  of  the  model 
is  shown  in  Figure  6.6.  The  lamina  has  several  fibers  oriented  at  an  angle  of  45®  to  the  load 
axis  which  is  the  most  accurate  angle  to  orient  a  specimen  for  shear  modulus  measurement 
[28].  The  inter-fiber  spacing  is  260  microns  corresponding  to  a  volume  fraction  of  35%. 
The  boundary  conditions  are  also  indicated  in  Figure  6.6.  As  discussed  earlier,  debonding 
is  simulated  using  the  matrix  layer  model.  Debonding  and  fiber  fracture  are  incorporated 
through  user  subroutines. 


6.5  Results  and  Discussion 


SCS-6/Ti-15-3  composite  is  studied  as  the  model  system.  SCS-6  is  a  silicon  carbide  fiber 
(140  microns  diameter)  made  by  Textron  Specialty  Materials.  Matrix  is  elastic-plastic  with 
linear  work  hardening.  The  fiber  and  matrix  properties  are  given  in  Table  6.1  (Mital  et.  al. 
[29],  Gambone  [15]). 


6.5.1  Longitudinal  Behavior 

Figure  6.7  (a)  shows  the  stress-strain  curves  resulting  from  a  simulation  conducted  for  a 
35  %  volume  fraction  lamina  (curve  A)  and  a  16  %  volume  fraction  lamina  (curve  B).  The 
elastic  portion  of  curve  A  ends  at  point  T,  after  which  the  change  in  slope  indicates  the 
onset  of  matrix  plasticity.  The  first  fiber  break  occurs  at  a  stress  value  very  close  to  the  peak 
value  in  the  figure.  Subsequent  breaks  occur  rapidly  within  a  very  small  increase  in  applied 
strain,  causing  a  steep  drop  in  the  stiffness  of  the  lamina.  This  indicates  a  catastrophic 
failure  of  the  lamina.  The  failure  of  the  lamina  immediately  after  the  first  break  is  confirmed 
from  experimental  observations  by  Majumdar  [13]  and  Lerch  and  Saltsman  [32].  Figure  6.7 
(b)  shows  the  schematic  of  the  location  of  fiber  breaks  that  occurred  in  the  35  %  volume 
fraction  case.  The  breaks  are  marked  by  a  circle  and  the  numbers  indicate  the  sequence  of 
breaks,  signifying  the  evolution  of  damage.  For  the  35  %  volume  fraction  lamina,  detailed 
finite  element  analyses  by  Voleti  et  al.  [20]  have  shown  that  the  stress  elevation  ini  the 
fiber  adjacent  to  a  break  occurs  in  the  same  plane,  increasing  the  probability  of  a  co-planar 
fracture.  This  co-planar  stress  elevation  is  magnified  when  plasticity  occurs. 

Figure  6.8  shows  the  fracture  pattern  for  the  16  %  volume  fraction  lamina.  The  fracture 
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pattern  is  at  an  angle  to  the  y-axis  (termed  as  ‘inclined  fracture’  in  this  work).  It  may 
be  mentioned  that  the  stochastic  distribution  of  fiber  flaws  for  both  the  volume  fractions 
are  identical.  Finite  element  analysis  of  the  fracture  plane  for  this  volume  fraction  [20]  has 
shown  that  the  peak  stress  due  to  fiber  fracture  occurs  at  an  angle  to  the  fracture  plane  due 
to  the  propagation  of  the  plastic  front  increasing  the  probability  of  a  non-coplanar  cracking. 

In  both  the  cases,  a  cluster  of  preferential  breaks  (co-planar  and  inclined  fracture)  caused 
the  failure  of  the  lamina.  These  results  agree  with  the  single  ply  tensile  experiments  [13]  on 
Titanium  matrix  ccinposites  having  SiC  as  the  reinforcement.  Metallurgical  examination  of 
the  tested  specimens  showed  the  presence  of  local  intense  slip  bands  (Figure  7.16  a)  in  the 
matrix  near  a  fiber  break.  The  fiber  spacing  in  this  case  is  about  0.7  mm  which  is  equivalent 
to  16.6%  fiber  volume  fraction.  Ultrasonic  Shear  Back  Reflection  analysis  (Figure  7.16  b) 
revealed  fiber  fractures  were  more  correlated  occurring  in  a  narrow  band.  Fiber  breaks  were 
predominantly  found  to  occur  at  locations  where  the  slip  bands  contact  the  neighboring 
fiber.  This  validates  the  numerical  predictions. 

Figure  6.10  shows  the  results  from  a  simulation  where  the  fiber  fractures  were  forced  to 
occur  in  a  pre-determined  pattern;  co-planar,  inclined  and  random.  These  load-displacement 
curves  reflect  the  qualitative  behavior  of  the  lamina.  The  stiffness  response  of  the  lamina 
is  distinctly  different  when  the  fractures  are  randomly  located  in  the  lamina  as  compared 
to  the  case  where  the  breaks  occur  preferentially.  The  co-planar  failure  shows  the  weakest 
lamina  response.  The  lamina  failure  involving  ‘inclined’  fracture  showed  marginally  better 
lamina  response  than  the  co-planar  case;  but  in  the  actual  simulation  little  difference  was 
observed,  and  catastrophic  failure  resulted  soon  after  the  first  break.  These  simulations  do 
not  include  the  effect  of  residual  stresses.  It  was  found  from  other  studies  that  inclusion 
of  residual  stresses  caused  the  initiation  of  matrix  plasticity  at  a  lower  strain  level  while 
delaying  the  strain  at  which  the  first  fiber  break  occurred.  The  stress-strain  behavior  was 
otherwise  similar. 

Effect  of  Interface 

To  study  the  role  ot  interface  in  the  propagation  of  fiber  fracture,  a  3-fiber  model  (Figure 
6.4)  of  a  composite  lamina  is  examined.  The  middle  fiber  in  this  system  was  fractured 
and  the  stress  elevation  due  to  the  fracture  on  the  neighboring  intact  fiber  was  examined. 
Debonding  reduces  the  stress  concentration  near  the  broken  end  of  the  fiber  and  causes  a 
reduction  in  the  maximum  plastic  strain.  However,  debonding  in  the  broken  fiber  increases 
the  load  in  the  adjacent  matrix  creating  diffused  plastic  zone  as  shown  in  Figure  7.14.  The 
volume  fraction  used  in  this  case  is  35%.  The  net  effect  is  a  reduction  in  the  peak  stress 
in  the  intact  fiber  but  a  greater  length  of  the  intact  fiber  experiences  higher  stresses  (figure 
7.15).  This  reduces  the  probability  of  coplanar  (0°)  fracture.  In  this  study  frictional  effects 
in  the  debonded  interface  are  neglected. 
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6.5.2  Transverse  Behavior 


The  transverse  modulus  of  the  lamina  model  is  shown  in  Table  6.2,  and  compared  with  values 
derived  from  theories  (Rule  of  mixtures  and  Tsai-Halpin  equations).  The  table  reflects  the 
undamaged  elastic  modulus.  There  is  a  good  match  between  the  Tsai-Halpin  and  finite 
element  results  for  the  case  of  strong  bond,  thus  validating  the  numerical  model.  The 
damaged  lamina  response  was  computed  for  two  different  values  of  interfacial  bond  strength, 
200  MPa  and  1500  MPa,  representing  weak  and  strong  bonds,  respectively.  The  stress-strain 
curves  for  weak  and  strong  interfacial  bonds  are  shown  in  Figure  6.13. 

In  the  case  of  higher  interfacial  bond  strength,  matrix  plasticity  initiates  much  earlier  to 
debonding  as  seen  by  the  first  non-linearity  in  the  curve  (Point  T).  Debonding  occurs 
subsequently  causing  a  major  drop  in  stress  values  (Point  A  to  B).  The  weak  bond  case 
shows  the  characteristic  3-stage  behavior  seen  in  the  experiments  [33].  The  first  stage  (upto 
Point  A')  is  completely  elastic.  The  second  stage  is  dominated  by  debonding.  Debondmg 
initiates  at  at  this  point  and  continues  until  the  total  debonding  occurs  at  Point  B  .  Local 
plasticity  in  the  matrix  is  seen  at  about  0.33%  and  steadily  increases  with  load.  The  third 
stage  (after  Point  B')  is  completely  dominated  by  the  plastic  deformation  of  the  matrix. 


6.5.3  Shear  Behavior 

The  shear  modulus  of  the  lamina  evaluated  from  the  model  is  shown  in  Table  6.2,  and 
compared  with  the  values  derived  from  two  different  theories.  Rule  of  mixtures  and  Tsai- 
Halpin  model.  There  is  a  good  match  between  the  two  values.  Figure  6.14  shows  the  results 
from  the  Monte-Carlo  simulation.  Cases  were  run  for  two  different  interfacial  bond  strengths 
of  200  MPa  and  1500  MPa.  These  bond  strengths  make  the  fiber  debond  before  and  after 
matrix  becomes  plastic,  corresponding  to  weak  and  strong  interfacial  bond  strengths.  Some 
of  the  observations  are: 

•  No  fiber  fracture  is  observed  in  the  simulations  indicating  that  debonding  and  matrix 
plasticity  are  the  only  mechanisms  invoked. 

•  In  the  case  of  weak  bond  strength,  debonding  initiates  at  a  longitudinal  strain  of  0.24 
%  and  matrix  plasticity  initiates  subsequently. 

•  For  higher  bond  strength,  matrix  plasticity  initiates  at  0.57  %  (figure  6.14)  and  debond¬ 
ing  initiates  subsequently. 

•  In  both  the  cases  failure  of  the  lamina  does  not  occur  even  at  2%  strain  agreeing  with 
the  experimental  observations  [32]. 
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6.5.4  Evaluation  of  the  damage  tensor 

Each  lamina  along  its  material  (fiber)  direction  is  assumed  to  be  transversely-isotropic  and 
requires  five  independent  constants  (four  in  the  plane  of  isotropy)  to  completely  describe  its 
elastic  behavior.  The  stiffness  tensor  for  two-dimensional  response  with  damage  is  given  in 
equation  6.3.  In  this  study  the  evolution  of  the  Poisson’s  ratio  D\2  is  not  considered.  Du, 
D22,  Dee  are  derived  from  longitudinal,  transverse  and  shear  tests  described  earlier.  The 
variation  of  Dn,  D22,  Dee  with  applied  load,  models  the  evolution  of  damage  in  each  coor¬ 
dinate  direction.  The  stress  state  at  any  applied  strain  can  be  determined  uniquely  using 
equation  6.2  by  knowing  the  evolution  of  the  Damage  parameters  (or  the  deterioration  of 
the  elastic  moduli)  and  the  inelastic  strain.  This  data  can  be  obtained  from  the  numerical 
predictions.  The  reduction  in  the  elastic  moduli  is  determined  by  numerically  simulating 
unloading  and  reloading  experiments  at  different  applied  strain  levels.  The  damage  param¬ 
eters  in  the  transverse  direction  (D22)  is  shown  in  in  Table  6.3  as  a  function  of  the  applied 
strain. 


6.6  Summary  and  Conclusions 


A  micromechanics  based  modeling  approach  is  developed  to  study  the  local  failure  mech 
anisms  such  as  fiber  fracture  and  interfacial  debonding  in  continuous  fiber  metal-matrix 
composites.  The  fiber,  matrbc  and  the  interface  are  modeled  explicitly.  Monte  Carlo  simu¬ 
lation  procedure  was  used  to  model  the  stochastic  nature  of  the  fiber  failure.  The  fracture 
patterns  predicted  from  the  simulation  of  0^-lamina  response  matched  well  with  experimen¬ 
tal  results.  The  effect  of  interfacial  debonding  on  the  stress  elevation  in  the  presence  of  a 
fiber  fracture  is  investigated.  It  is  found  that  debonding  results  in  diffused  plastic  region 
in  the  matrix  near  the  fiber  break  thereby  rhifting  the  location  of  the  peak  stress.  The 
evolution  of  damage  under  different  interfacv  i  conditions  is  studied  using  the  models.  A 
damage-dependent  constitutive  tensor  is  developed  for  the  lamina  from  the  numerical  simu¬ 
lations.  This  approach  will  help  in  studying  the  influence  of  micromechanical  parameters  on 
the  evolution  of  damage  in  composites  not  only  under  static  loads  but  also  when  subjected 
to  cyclic  loading. 
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Material  Properties 

Weibull  Properties 

Fiber 

E  =  428  GPa,  v  =  0.3,  CTE  =  3.56  E-06  /^C 

Matrix 

E  =  85  GPa,  V  =0.32,  Yield  =  760  MPa,  CTE  =  7.812  E-06  /^C 

Work  Hardening  =  3.2  E3  MPa 

Table  6.1:  Material  properties  of  the  fiber  and  matrix 


Theoretical  Model 

FEM  Model 

0  degree  lamina 
(longit.  modulus) 

205  GPa  (Rule  of.  Mix) 

205  GPa 

90  degree  lamina 
(transverse  modulus) 

149.0  GPa  Tsai-Halpin 

118.1  GPa  Rule  of  Mix. 

155  GPa 

45  degree  lamina 
(shear  modulus) 

44.81GPa  Rule  of  Mix. 

52.0  GPa  Tsai-Halpin 

45.35  GPa 

Table  6.2:  Elastic  modulli  calculated  from  FEM  and  theoretical  solutions  . 
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Finite  Element 
Mesh  of  Lamina 


Figure  6.1:  Flow-chart  of  the  Monte  Carlo  simulation 


95 


1 


i-r,  2-2' 

duplicate  nodes 

(with  the  same  co-ordinates) 

at  the  interface 

1-2,  r-2' 

Contact  edges  which  cannot 
penetrate 

n,  t 

normal  and  tangential 
directions 

kn.kt 

Stiffness  of  springs  in  normal 
and  tangential  directions 

Figure  6.2:  Spring-layer  model  for  fiber-matrix  interface 
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Figure  6.3:  Schematic  of  the  0°  lamina  and  boundary  conditions 
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Figure  6.4:  Schematic  of  the  3-fiber  model  and  boundary  conditions 
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Figure  6.5:  Schematic  and  boundary  conditions  of  the  finite  element  model  to  study  trans¬ 
verse  behavior 
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Figure  6.6:  Schematic  diagram  of  the  loading  and  boundary  conditions  of  the  45°  ply  model 
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(a) 


(b) 

Figure  6.7:  (a)  Stress-strain  curves  for  0°  degree  laminae,  (b)  fracture  pattern  for  a  35  % 
volume  fraction  lamina 
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Figure  6.10:  Load-displacement  curves  for  pre-determined  damage  patterns 
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Figure  6.11:  Propagation  of  plastic  front  in  the  case  of  broken  fiber  with  debonding 
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Ratio  of  stress  elevation  (o/ct  ) 


Figure  6.12;  Variation  of  axial  stress  in  the  intact  fiber  upon  debonding  of  the  broken  fiber 


Figure  6.13:  Stress-strain  curve  for  a  weakly-bonded,  transversely-loaded  composite  • 
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Chapter  7 


Effect  of  Fiber  Fracture  and  Matrix 
Yielding  on  Load  Sharing  in 
Continuous  Fiber  Metal  Matrix 
Composites 

7.1  Introduction 

The  tensile  response  of  a  composite  lamina  is  governed  by  the  strength  of  the  fiber  and  the 
load  transfer  properties  of  the  matrix.  In  continuous  fiber  composites,  individual  fibers  can 
fail  well  before  the  ultimate  failure  of  the  composite.  Any  fiber  break  causes  a  re-distribution 
of  load  [1],  to  other  fibers  causing  stress  elevation  r>u  the  adjacent  fibers.  This  results  in  a 
non-uniform  stress  distribution  along  the  length  of  those  fibers.  The  strength  of  the  fiber 
also  varies  randomly  across  its  length  due  to  the  stochastic  failure  behavior  of  the  brittle 
fibers.  Under  these  conditions,  the  fiber  break  occurs  not  necessarily  at  the  weakest  link, 
but  where  the  stress  exceeds  the  local  strength  of  the  fiber. 

The  effect  of  an  isolated  fiber  break  on  the  occurrence  of  subsequent  breaks,  and  their 
effect  on  the  stiffness  of  the  composite  are  important  issues.  Different  rules  for  load  sharing 
have  been  proposed  in  this  context,  such  as  equal  (or  global)  load  sharing  [2]  and  local  load 
sharing  [3].  The  prevalence  of  one  over  the  other  is  determined  by  various  factors  including 
the  fiber  and  matrix  properties  and  fiber  volume  fraction.  In  the  case  of  local  load  sharing, 
the  stress  elevation  in  the  neighboring  fibers  caused  by  load  re-distribution  could  be  sufficient 
to  create  preferential  fiber  fractures  or  coordinated  fractures  (such  as  coplanar  or  angled), 
leading  to  premature  failure  of  the  compsoite.  Thus,  a  study  of  the  load  re-distribution  is 
critical  to  the  evaluation  of  the  tensile  properties  of  a  lamina. 

The  problem  of  load  transfer  in  composites  has  been  studied  using  shear-lag  theory  by 
Rosen  [2]  and  others  [4,  5,  6].  Most  of  these  analyses  are  developed  for  Polymer  Matrix 
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Composites  (PMC)  where  the  assumptions  of  shear-lag  theory  are  valid.  Metal  Matrix 
Composite  (MMC)  which  is  the  focus  of  the  present  study  significantly  differs  from  the 
PMC’s  in  that  the  behavior  of  fiber,  matrix  and  interface  is  quite  different.  There  are 
issues  such  as  interfacial  debonding  [7],  processing  induced  residual  stresses  [8]  and  the 
inelastic  deformation  of  the  matrix  [1]  that  are  specific  to  MMCs.  Nicholas  et  al.  [10] 
studied  the  problem  of  strain  accumulation  from  fiber  breaks  in  MMCs  during  in-phase 
thermomechanical  fatigue.  Iremonger  et.  al.[ll]  performed  simplified  computational  analysis 
of  the  problem.  He  et.  al.  [12]  used  linear  elastic  model  to  study  the  stress  elevation  in  the 
next,  and  the  next-to-next  fibers  in  the  presence  of  a  broken  fiber.  Probablistic  arguments 
were  then  used  to  determine  whether  global  load  sharing  or  local  load  sharing  would  occur 
in  those  composites.  But  the  effect  of  local  plasticity  and  debond  length  on  the  stress  re¬ 
distribution  were  not  addressed  in  that  study.  Majumdar  et.  al.  [13]  performed  0°  tensile 
tests  on  a  SiC/Ti-6-4  lamina  to  study  the  effect  of  fiber  breakage  and  matrix  plasticity  on  the 
fracture  behavior  of  the  lamina.  Their  experimental  results  showed  the  presence  of  extensive 
slip  band  formation  due  to  fiber  breaks.  The  impingement  of  those  slip  bands  on  adjacent 
fibers  caused  the  latter  to  fail  at  the  location  of  impingement.  The  presence  of  slip  bands 
indicate  that  plasticity  might  have  a  significant  effect  on  local  load  sharing  in  the  case  of 
MMCs.  The  objective  of  the  present  study  is  to  examine  the  role  of  matrix  plasticity  and 
interfacial  debonding  on  the  global/local  load  sharing  characteristics  of  MMCs. 


7.2  Computational  Model 


A  3-fiber  unit  cell  is  selected  to  study  the  stress  elevation  in  the  adjacent  fiber  in  the  pres¬ 
ence  of  a  broken  fiber  in  the  middle.  Due  to  symmetry  only  one  half  of  the  unit  cell  is 
considered  as  shown  in  Figure  7.1.  8-noded  plane  stress  quadrilateral  elements  are  used 
with  boundary  conditions  as  shown  in  Figure  7.1.  Plane  stress  assumption  represents  the 
closest  2-D  approximation  of  the  thin  section  with  layers  of  fibers  and  matrix  (as  shown  in 
the  figure).  MARC  [6],  a  non-linear  finite  element  code,  is  used  to  perform  the  numerical 
simulations.  Duplicate  nodes  are  created  along  the  fiber-matrix  interface  and  are  connected 
using  high  stiffness  springs  to  create  a  bonded  interface.  Debonding  is  simulated  by  reducing 
the  spring  constant  to  a  very  low  value.  Contact-friction  algorithm  inherent  in  the  software 
is  used  to  model  the  post-debonding  friction  and  also  to  prevent  the  penetration  of  nodes  at 
the  contact  surface.  This  model  has  been  sucessfully  used  in  a  number  of  micromechanical 
studies  of  interfacial  debonding  [7].  Nodal  displacements  are  applied  to  the  top  edge  of  the 
model  to  simulate  the  application  of  uniform  uniaxial  strain. 


7.3  Results  and  Discussion 


The  three-fiber  unit  cell  model  is  used  to  study  the  effect  of  stress  redistribution  due  to  a 
fiber  break  on  the  adjacent  fiber.  SiC/Ti-15-3  MMC  system  is  considered  in  this  study.  The 
fiber  and  matrix  properties  are  shown  in  Table  1  [15].  The  fiber  diameter  is  140  microns. 
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Fiber  spacing  is  computed  based  on  the  load  carrying  capacity,  thus  the  area  ratio  of  the  fiber 
to  total  (fiber+matrix)  represents  the  fiber  volume  fraction.  If  s  is  the  distance  betweenj;he 
edges  of  the  two  fibers  and  df  is  the  diameter  of  the  fiber,  i.e,  the  volume  fraction  Vf  df+s' 
Fiber-matrix  interface  is  perfectly  bonded  prior  to  the  propagation  of  fracture  along  the 
interface.  Two  types  of  constitutive  models  are  considered  for  the  matrix  material:  (i)  linear 
elastic  (ii)  elastic-perfectly  plastic.  First  the  results  from  purely  elastic  stress  analysis  are 
discussed  followed  by  a  section  devoted  to  elastic-plastic  behavior. 


7.3.1  Effect  of  Fiber  Fracture:  Elastic  Analysis 

Fiber  fracture  causes  redistribution  of  stresses  in  the  vicinity  of  the  break  resulting  in  a 
localized  zone  of  high  stresses  in  the  neighboring  fibers.  This  is  primarily  due  to  the  loss  in 
the  load-carrying  ability  near  the  broken  end  of  the  fiber.  This  is  evident  from  the  results 
from  Figure  7.2  which  show  a  drop  in  the  axial  stress  along  the  length  of  the  broken  fiber. 
It  can  be  seen  from  this  figure  that  the  axial  stress  drops  from  the  far  field  value  to  almost 
zero  at  the  broken  end.  The  neighboring  fibers  have  to  take  ths  extra  load  due  to  the  fiber 
break.  Results  from  the  three-fiber  model  also  indicate  that  the  peak  stress  along  the  fiber 
axis  (Figure  7.3)  in  the  intact  fiber  shows  a  11%  increase  when  compared  to  the  far-field 
fiber  axial  stress  ao  which  is  the  same  in  the  broken  as  well  as  the  intact  fiber.  In  order  to 
determine  the  effect  of  the  broken  fiber  on  the  second  neighbor,  the  analysis  was  conducted 
using  a  separate  model  with  two  neighboring  fibers.  It  was  found  that  the  stress  elevation  on  . 
the  second  neighbor  was  less  than  1%  for  a  fiber  volume  fraction  of  35%.  These  effects  are 
expected  to  decrease  further  for  lower  volume  fractions.  Hence,  the  3-fiber  (or  one  neighbor) 
model  is  chosen  in  all  the  cases  considered  in  this  study.  It  has  to  be  noted  that  the  stress 
elevation  from  the  above  analyses  are  much  lower  than  that  predicted  by  the  shear-lag  theory 
[1],  which  neglects  the  axial  stresses  in  the  matrix.  With  decrease  in  the  stiffness  of  the 
matrix,  the  peak  stress  in  the  intact  fiber  increases  to  reach  the  levels  predicted  by  the 
shear-lag  analysis,  consistent  with  the  assumptions  in  the  theory. 

The  stress  elevation  in  the  neighboring  fibers  is  also  caused  by  the  stress  concentration 
effects  at  the  crack  tip  and  due  to  bending  of  the  neighboring  fiber.  Fiber  fracture  causes 
stress  concentration  at  the  edge  of  the  fractured  surface.  The  variation  of  the  axial  stresses 
along  the  fracture  plane  is  shown  in  Figure  7.4.  It  can  be  seen  that  the  stress  is  maximum 
near  the  crack  tip  and  drops  rapidly  with  increasing  distance  away  from  the  broken  end. 
To  capture  the  effects  of  the  stress  singularity  adequately  the  crack  tip  was  modeled  using 
singular  elements  (collapsed  8-noded  quadrilaterals  with  quarter  point  elements).  Compari¬ 
son  of  the  results  indicated  that  though  the  stresses  at  the  broken  end  of  the  fiber  increases 
sharply  with  the  use  of  singular  elements,  the  stresses  in  the  neighboring  fiber  show  negligible 
change  (for  fiber  volume  fraction  of  35%).  Based  on  this  analysis  it  can  be  concluded  that 
for  the  fiber  spacing  under  consideration  special  singular  elements  are  not  needed  to  study 
the  stress  concentration  effects  on  adjacent  fibers,  and  hence  not  used  subsequently. 

Fiber  bending  is  caused  by  the  shear  strain  gradients  due  to  the  relative  displacement 
between  the  broken  and  the  intact  fiber  in  the  loading  direction.  Figure  7.5  shows  the 
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variation  of  the  x-displacements  of  the  intact  fiber  along  its  axis.  As  can  be  seen,  the  x- 
displacement  is  maximum  along  the  fracture  plane  resulting  in  a  maximum  bending  stress 
at  location  B.  This  also  accounts  for  the  non-uniformity  of  the  axial  stresses  in  this  fiber 
along  the  x-axis  (Figure  7.6)  with  the  peak  on  the  left  edge  where  the  bending  stresses  are 
tensile.  It  is  clear  from  the  above  that  the  peak  stress  in  the  intact  fiber  is  the  net  effect  of 
the  above  local  mechanisms.  In  this  analysis  a  debonding  of  the  order  of  0.1  r/  is  always 
assumed  to  be  present  at  the  interface  around  the  broken  fiber,  to  relieve  the  high  levels  of 
shear  stresses  present  at  the  crack  tip. 


7.3.2  Effect  of  Matrix  Plasticity 

Inelastic  behavior  of  the  matrix  affects  the  stress  redistribution  due  to  a  fiber  break  in  many 
ways.  Local  plastic  deformation  near  the  fiber  break  creates  a  crack  shielding  effect  that 
minimizes  the  effect  of  stress  concentration  and  also  avoids  matrix  cracking.  Though  fiber 
fracture  may  occur  arising  firom  a  stochastic  distibution  of  strength,  in  general,  some  fibers 
may  have  very  high  strength  leading  to  matrbc  yielding  prior  to  any  fiber  breakage.  However, 
when  a  fiber  break  occurs  in  the  early  loading  stages,  local  plasticity  occurs  much  before 
bulk  matrbc  yielding.  With  increase  in  the  applied  strain  the  plastic  front  grows  at  an  angle 
towards  the  adjacent  fiber.  Figure  7.7  shows  the  propagation  of  plastic  strains  for  the  case 
of  35%  volume  fraction.  With  the  initiation  of  plasticity  the  ratio  of  peak  stress  in  the 
intact  fiber  to  the  far  field  stress  (the  stress  elevation  ratio,  q)  increases  and  deviates  from 
the  elastic  case  where  it  has  a  constant  low  value  independent  of  the  applied  strain  (Figure 
7.8).  This  is  anticipated  since  the  plastic  regions  cannot  take  any  additional  load,  increasing 
the  load  on  the  fiber.  It  can  also  be  seen  from  Figure  7.7  and  7.8  that  as  the  plastic  frorit 
approaches  the  inteict  fiber  (applied  strains  of  >  0.6%)there  is  a  significant  increase  in  q  . 
At  strain  levels  of  1%  global  plasticity  sets  in  and  ‘q’  does  not  increase  any  further.  It  should 
be  noted  that  the  sti  ^s  elevation  in  this  case  is  much  higher  (maximum  value  of  the  order 
of  50%)  when  compc-xcd  to  11%  (see  Figure  7.8)  in  the  elastic  case. 

Figure  7.9  shows  the  variation  of  the  plastic  deformation  of  the  matrix  along  the  length 
of  the  specimen  after  the  plastic  front  has  reached  the  intact  fiber.  It  can  be  seen  that  for 
this  volume  fraction,  though  the  plastic  front  propagates  at  an  angle,  the  equivalent  plastic 
strain  near  the  interface  of  the  intact  fiber  along  the  fracture  plane  (A-B  in  the  figure)  is 
almost  75%  of  the  peak  value  which  occurs  at  a  distance  above  the  lower  end  (B).  This  is 
very  much  different  from  the  case  of  16%  volume  fraction  to  be  dicussed  in  the  next  section. 
This  is  because  at  this  high  volume  fraction  the  fibers  are  close  enough  that  by  the  time  the 
plastic  front  reach  the  neighboring  intact  fiber,  plasticity  has  initiated  in  the  whole  region 
between  the  two  fibers  up  to  the  fracture  plane  plane  (Figure  7.7).  Under  this  conditmn  ‘q’ 
is  found  to  be  maximum  at  B  making  it  the  location  which  is  the  most  probable  to  fail  next 
(Figure  7.10). 

Figure  7.11  shows  the  peak  stress  along  the  intact  fiber  for  a  fiber  volume  fraction  of 
16%.  Figure  shows  that  the  peak  stress  for  this  case  occurs  at  an  angle  of  about  45°  from 
the  fratnre  plane.  This  can  be  explained  from  the  propagation  of  the  plastic  front  shown  in 
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Figure  7.12.  As  in  the  case  of  35%  volume  fraction,  the  plastic  front  grows  at  an  angle  to  the 
fracture  plane  and  reaches  the  intact  fiber.  It  can  be  seen  that  at  this  stage  the  matrix  in 
the  vicinity  of  the  intact  fiber  along  the  fracture  plane  shows  negligible  plastic  strain  which 
is  very  much  different  from  the  35%  case.  This  is  also  clear  from  the  variation  of  plastic 
strain  shown  in  Figure  7.13.  Moreover  at  this  low  volume  fraction  levels  the  effect  of  stress 
concentration  due  to  the  fiber  break  on  the  adjacent  fiber  is  also  considerably  reduced.  This 
makes  matrix  plasticity  as  the  predominant  factor  in  causing  the  stress  elevation. 


7.3.3  Effect  of  Interfacial  Debonding 

Interfacial  debonding  followed  by  a  fiber  break  modifies  the  local  stress  state.  Debonding 
reduces  the  stress  concentration  near  the  broken  end  of  the  fiber  and  causes  a  reduction  in 
the  maximum  plastic  strain.  However,  debonding  in  the  broken  fiber  increases  the  load  in 
the  adjacent  matrbc  creating  a  diffused  plastic  zone  as  shown  in  Figure  7.14.  A  debonded 
length  of  0.3  Iq  is  used  in  this  case.  The  net  effect  is  a  reduction  in  the  peak  stress  in  the 
intact  fiber  while  a  greater  length  of  the  intact  fiber  is  subjected  to  higher  stresses  (Figure 
7.15).  This  reduces  the  probability  of  coplanar  (0°)  fracture.  In  this  study  frictional  effects 
in  the  debonded  zone  are  neglected. 


7.3.4  Experimental  Comparisons 

Majumdar  et.  al.  [13,  16]  conducted  single  ply  tensile  experiments  on  titanium  matrix  com¬ 
posites  having  SiC  as  the  reinforcement.  Metallurgical  examination  of  the  tested  specimens 
showed  the  presence  of  local  intense  slip  bands  (Figure  7.16  a)  in  the  matrbc  near  a  fiber 
break.  The  fiber  spacing  in  this  case  is  about  0.7  mm  which  is  equivalent  to  16.6%  fiber  vol¬ 
ume  fraction.  Ultrasonic  Shear  Back  Reflection  analysis  (Figure  7.16  b)  revealed  correlated 
fiber  fractures  occured  in  a  narrow  band.  Fiber  breaks  were  predominantly  found  to  occur 
at  locations  where  the  slip  bands  contact  the  neighboring  fiber.  This  validates  the  numerical 
predictions.  This  also  clearly  indicates  the  prevalence  of  local  load  sharing  compared  to 
global  load  sharing  behavior  (where  the  fiber  fractures  are  expected  to  occur  randomly). 


7.3.5  Summary  and  Conclusions 


The  stress  elevation  on  the  adjacent  fiber  in  the  presence  of  a  broken  fiber  in  a  MMC 
lamina  subjected  to  longitudinal  tension  was  examined  using  micromechanical  models.  The 
role  of  various  local  mechanisms  such  as  the  crack  tip  stress  concentration  effect,  interfacial 
debonding  due  to  the  fiber  break  were  studied  for  the  case  of  a  ductile  matrix.  The  stress 
elevation  was  found  to  be  due  to  the  net  effect  of  all  these  mechanisms.  At  lower  volume 
fractions  the  stress  concentration  in  the  vicinity  of  the  crack  tip  (due  to  the  fiber  break) 
was  found  to  have  minimum  effect  on  the  stress  elevation  in  the  adjacent  fiber.  However  it 
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initiates  local  matrix  plasticity  which  propagates  towards  the  adjacent  fiber  much  earlier  to 
the  occurrence  of  global  matrbc  plasticity.  This  results  in  a  significant  increase  in  the  stress 
elevation  in  that  region.  At  higher  volume  fractions  the  probability  of  coplanar  firacture 
was  found  to  increase.  Interfacial  debonding  was  able  to  decrease  the  magnitude  of  stress 
elevation  and  also  decreased  the  probability  of  coplanar  breaks  which  cause  catastrophic 
failure  of  the  composite.  The  results  from  this  study  strongly  favor  the  existence  of  local 
load  sharing  in  MMCs  at  higher  volume  frac^^' ons. 
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Table  7.1:  Material  properties  of  fiber  and  matrix 
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Figure  7.1:  Finite  element  model  with  boundary  conditions 
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Figure  7.4:  axial  elastic  x 
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Figure  7.9:  variation  pi  str  fiber  length  vf  35 
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Figure  7.11:  variation  stress  fiber  length  vf  16 
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Figure  7.12:  propagation  vf  16 
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Figure  7.14:  plasticity  with  debonding 
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Chapter  8 

Effect  of  Interfacial  Properties  on  the 
Fiber  Fragmentation  Process  in 
Polymer  Matrix  Composites 


8.1  Introduction 


The  Fiber  Fragmentation  Test  (also  called  Single  Fiber  Composite  Test,  SFC)  has  become 
a  commonly  used  method  of  characterizing  interfaces  in  composites  and  evaluating  quality 
of  bond  strength.  First  developed  by  Kelly  and  Tyson  [1]  to  study  copper-tungsten  com¬ 
posites,  this  technique  has  been  applied  to  Polymer  Matrbc  Composites  (PMC)  for  better 
understanding  the  mechanics  and  statistics  of  fiber  fragmentation  process  [2,  3].  The  testing 
procedure  relies  on  multiple  fragmentation  of  a  single  fiber  embedded  in  a  matrbc  under  the 
application  of  an  axial  strain.  The  axial  strain  is  progressb-dy  increased  until  a  ‘saturation’ 
point  is  reached,  when  there  are  no  more  fiber  breaks.  The  test  yields  a  number  of  fiber 
fragments  which  is  related  to  the  fiber-matrix  interfacial  shear  strength.  Interpretation  of 
the  fragmentation  data  is  generally  based  on  the  concept  of  critical  fiber  length,  Ic,  and  a 
suitable  load  transfer  model  such  as  the  shear  lag  theory.  The  averap  shear  strength  r  of 
the  interface  is  then  estimated  from  the  equation  r  =  ^  where  as  is  the  tensile  strength 
of  the  fiber  at  the  critical  len^h  and  d  is  the  fiber  diameter.  Ic  is  derived  from  the  mean 
fragment  length  Im,  at  saturation. 

A  large  literature  base  available  in  the  application  of  SFC  to  Polymer  Matrix  Compos¬ 
ites  [4,  5,  6,  7,  8,  9,  10].  Though  it  has  been  accepted  that  SFC  gives  useful  information 
about  the  inter  facial  bond  quality  and  the  fiber  strength,  there  is  no  clear  consensus  on 
how  to  interpret  the  data  and  quantify  the  interfacial  properties  [llj.  The  average  shear 
strength  given  by  Kelly-Tyson  equation  based  on  the  mean  fragment  length  at  saturation 
from  the  experiment  is  the  most  common  result  obtained  from  the  test;  but,  this  is  purely  a 
qualitative  measure  and  does  not  consider  the  history  of  the  fragmentation  process  and  ne¬ 
glects  factors  such  as  interfacial  debonding  which  can  significantly  influence  the  test  results. 
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To  understand  the  fragmentation  process  better,  numerical  modeling  of  the  test  has  been 
adopted  recently  [4,  12,  7].  It  has  been  also  recognized  that  for  better  understanding  of  the 
failure  process  during  SFC,  it  is  necessary  to  obtain  the  test  data  not  only  at  saturation,  but 
also  at  various  levels  of  applied  strain  [12]. 

Interfacial  bond  strength  and  frictional  stress  play  a  critical  role  in  the  load  transfer  in 
a  single-fiber  composite  test;  hence,  understanding  the  effect  of  these  parameter^  is  crucial  to 
predicting  the  test  response.  Thus  far,  there  has  been  little  effort  at  quantifying  the  effect  of 
these  parameters  on  the  fragmentation  characteristics.  The  influence  of  interracial  properties 
(bond  strength  and  frictional  stress)  on  the  progression  of  the  fragmentation  process  is 
investigated  in  this  work  using  numerical  simulation.  A  parametric  study  is  conducted  using 
the  model  to  quantify  the  effect  of  interfacial  bond  strength  and  frictional  sliding  stress  on 
the  mean  fragment  length  and  debonding  progression.  Some  of  the  numerical  issues  relating 
to  the  simulation  of  the  SFC  test  are  also  addressed. 


8.2  Interface  Model 


In  this  study  the  load  transfer  at  the  fiber-matrix  interface  is  modeled  using  shear-lag  theory 
[13].  The  axial  and  shear  stress  equations  firom  the  shear-lag  model  are  [14,  15]: 


cTf  =  eEf{l  — 


coshjPi^  -  x)) 
cosh{l3^) 


(8.1) 


r  =  0.beEfrf(3 


sinh{l3{^  —  x)) 
cosh{P\) 


(8.2) 


where: 


a  —  /  2.QGm  m/2 


(8.3) 


and  0  <  X  <  /.  In  the  above  equations,  aj  is  the  axial  stress  in  the  fiber,  e  is  the  applied 
strain  and  T  is  the  fragment  length.  Ef  is  the  fiber  modulus  and  Gm  is  the  matrix-shear 
modulus.  Tf  is  the  fiber  radius  and  ‘R’  is  the  matrix  width  (taken  about  50  r/).  The  iilterface 
is  assumed  to  be  perfectly  bonded  in  deriving  these  equations.  The  axial  and  shear  stress 
profiles  at  the  interface  are  shown  in  Figure  8.1.  It  can  be  seen  from  Figure  8.1(a)  that  at 
both  the  ends  of  the  fragment  there  is  a  region  in  which  the  axial  stress  rises  exponentially 
reaching  the  far  field  value  at  the  end  of  this  region.  This  region  where  the  stresses  are  lower 
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than  the  far  field  value  is  known  as  the  ineflFective  length.  The  shear  stresses  drop  to  zero 
at  the  end  of  the  ineffective  length  as  shown  in  Figure  8.1(b).  in  the  above  equation  is  a 
material  constant  and  influences  the  steepness  of  the  exponential  region  of  the  axial  stress 
curve.  As  can  be  seen  from  equation  8.3,  /3  is  a  function  of  the  shear  modulus  of  the  matrix 
and  represents  the  effectiveness  of  the  load  transfer  between  the  fiber  and  the  matrix.  A 
higher  0  leads  to  a  rapid  increase  in  the  axial  stress,  thereby  reducing  the  ineffective  length. 

Interfacial  debonding  is  simulated  based  on  maximum  shear  strength  criterion.  Debond¬ 
ing  results  in  drop  in  the  stress  levels  increasing  the  ineffective  length.  In  the  debonded 
region  axial  force  equilibrium  arguments  yield  the  following  stress  distributions.  A  constant 
frictional  shear  stress  condition  is  assumed  in  this  region. 

af  =  2rfrx/rf  (8.4) 

T  =  Tfr  (8.5) 

where  0  <  x  <  Id-  ld  =  ml/2  is  the  debonded  length  where  ‘m’  is  the  debonding  ratio.  The 
length  of  debonding  is  determined  based  on  the  axial  force  equilibrium  and  by  imposing  the 
continuity  of  the  axial  stresses  at  the  end  of  the  debonded  zone.  In  the  bonded  region  shear 
lag  equations  are  applicable  and  are  modified  to  account  for  the  debonded  length  [15]. 


af  =  eEf  -  [eEf  -  2TfrSm] 


cosh{0{^  —  x)) 
cosh{0rfs{l  —  m)) 


(8.6) 


and. 


T  =  O.50rf[€Ef  -  2TfrSm] 


sinh{0{^  —  x)) 
cosh{0rfs{\  —  m)) 


(8.7) 


where  s  =  5  is  the  aspect  ratio  (where  d  is  the  diameter  of  the  fiber).  The  debond  ratio  ‘m’ 
is  determined  by  iteratively  solving  the  following  equation  [15]: 


m  — 


eEf  —  2^coth{0rfs{l  —  m)) 
2.0TfrS 


(8.8) 


The  axial  stress  and  shear  stress  profiles  in  the  presence  of  debonding  are  shown  in 
Figure  8.2.  As  can  be  seen  from  equation  8.8  the  debonding  ratio  ‘m  is  a  function  of  the 
bond  strength  Tdeb  and  the  frictional  stress  r/r-  The  significance  of  these  two  important 
interfacial  parameters  will  be  discussed  in  detail  in  section  8.5.  The  frictional  stress  r/r  =  tip 
where  p,  is  the  coefficient  of  friction  and  ‘p’  is  the  radial  stress.  The  radial  stress  gets  the 
contribution  from  three  possible  sources;  thermal  residual  stress,  differential  Poisson  effect 
and  the  fiber  surface  asperities. 


126 


8.3  Statistical  Model  for  Fiber  Strength  Distribution 


The  two-parameter  Weibull  probability  theory  is  used  to  extrapolate  the  fiber  strength  from 
the  test  gage  length  to  the  desired  unit  length  used  in  the  simulation.  According  to  this 
theory  the  probability  that  a  fiber  of  length  ‘L’  breaks  at  a  stress  lower  than  a  is  given  by: 


P{a)  =  1  -  exp{-^[—f)  (8-9) 

where  ‘q’  is  the  Weibull  modulus  (shape  factor),  Lo  is  a  reference  length  which  can  be 
arbitrarily  chosen,  and  Go  is  the  scale  parameter  for  unit  length  ratio  (;^  =  1). 

The  variation  of  fiber  strength  along  the  length  of  the  fiber  due  to  the  distribution  of 
flaws  is  taken  into  account  by  introducing  the  Weibull  Model  in  the  frame  work  of  Monte 
Carlo  simulation.  The  intial  length  of  the  fiber  is  divided  in  to  a  number  of  small  links. 
Random  values  of  link  strengths  are  derived  by  inverting  equation  8.9.  Lo  is  taken  as  unity. 
Thus 

a(i)  = 

where  R{i)  is  a  random  number  between  zero  and  unity  which  represents  the  probability  of 
failure  of  i  th  link  and  ‘i’  is  the  number  of  the  fiber  link. 


8.4  Implementation  of  the  model 

The  numerical  procedure  is  outlined  in  the  flow-chart  in  Figure  8.6.  The  fiber  is  divided  into 
a  numbe'r  of  small  links  (2500  links  for  all  the  cases  run  here  and  was  found  to  be  adequate). 
The  average  strength  of  one  link  is  evaluated  by  extrapolating  from  the  value  at  the  test  gage 
length  using  equation  ??.  The  statistical  distribution  of  strength  values  along  the  length  of 
the  fiber  is  given  by  equation  8.10.  The  geometric,  material  and  statistical  parameters  of  the 
fiber  and  matrix  are  input  into  the  program.  These  include  the  Young’s  modulii,  Weibull 
scale  parameter  and  Weibull  modulus  (shape  faxitor),  Poisson  ratio  and  radii  of  the  fiber  and 
matrix;  and  the  initial  length  of  the  fiber. 

A  strain  is  initially  applied  to  the  un-broken  fiber  and  the  resulting  stress  in  each  link 
is  compared  with  its  respective  strength.  If  the  stress  induced  exceeds  the  link  strength,  a 
break  occurs  in  this  link,  resulting  in  the  formation  of  two  fragments.  At  this  point,  a  check 
is  made  to  see  if  any  debonding  takes  place  in  the  entire  fiber.  If  there  is  any  debondihg,  the 
debond  length  is  calculated  through  the  iterative  solution  of  equation  8.8.  Then,  the  axial 
stresses  are  calculated  for  the  bonded  and  debonded  regions  of  the  fragment.  The  strain  is 
now  incremented  and  the  procedure  is  repeated.  This  iterative  procedure  is  repeated  until 
no  further  breaks  occur  with  increase  in  strain.  A  state  of  ‘saturation  is  reached  at  this 
strain  level. 
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8.5  Results  and  Discussion 


The  numerical  simulation  of  the  fiber  fragmentation  test  is  conducted  for  a  carbon-epoxy 
system.  The  material  properties  of  the  fiber  and  matrix  and  the  Weibull  parameters  of  the 
fiber  [14]  are  given  in  Table  1.  A  specimen  of  gage  length  10  mm  is  used  in  all  the  cases. 
To  study  the  effect  of  interfacial  conditions  on  the  fragmentation  test  results,  two  extreme 
conditions  are  considered:  a  perfectly  bonded  interface  (rjeb  =  oo)  and  a  weak  interface 
with  a  low  Tdeb  and  Tfr-  Perfect  bonding  though  a  hypothetical  situation,  is  considered  here 
mainly  to  bring  out  the  effect  of  debonding  on  fragmentation  process  by  comparing  the  two 

cases. 


8.5.1  Stress  Distributions 

Figure  8.6  shows  the  axial  stress  profile  in  the  fiber  as  a  function  of  increasing  applied  strain 
for  the  case  of  strongly  bonded  interface.  It  can  be  seen  that  the  number  of  breaks  continue 
to  increase  with  increase  in  applied  strain  till  the  fragment  lengths  are  smaller  than  the  load 
transfer  length  required  for  the  axial  strain  (or  stress)  to  reach  the  far  field  value  (Figure  8.6 
b);  further  breaks  are  then  considerably  reduced.  When  the  bonding  at  the  interface  is  not 
perfect,  debonding  occurs  in  the  regions  of  maximum  shear  stress  which  are  at  the  end  of 
the  fragments.  A  bond  strength  of  60  MPa  is  used.  The  debonded  length  is  calculated  using 
the  equations  in  section  8.2.  Figure  8.6  shows  the  axial  stress  profiles  from  the  simulation  for 
the  weakly  bonded  case.  Debonding  reduces  the  effective  length  available  for  load  transfer, 
leading  to  a  drop  in  the  fiber  stresses.  Debonding  progresses  with  increasing  applied  strain; 
this  reduces  the  probability  of  further  breaks  resulting  in  saturation  at  a  lower  applied  strain 
when  compared  with  the  perfectly  bonded  case.  It  can  be  seen  from  Figure  8.6(b)  that  a 
significant  portion  of  the  specimen  length  has  debonded  at  this  stage. 


8.5.2  Mean  Fragment  Length  and  Saturation 

Numerical  simulation  shows  that  the  variation  of  cumulative  number  of  breaks  (or  fiber  frag¬ 
ments)  with  increasing  applied  strain  is  significantly  different  for  weak  and  strong  interfaces. 
For  the  case  of  perfect  bonding,  the  rate  of  increase  of  fiber  breaks  with  increasing  applied 
strain  is  almost  uniform  as  can  be  seen  from  Figure  8.6.  It  can  also  be  seen  that  saturation 
of  breaks  is  not  seen  even  at  10%  strain.  Figure  8.6  shows  the  value  of  mean  fragment  length 
Ijn  (=  total  length/  cum.  number  of  breaks)  as  a  function  of  applied  strain  for  both  strong 
and  weak  interfaces.  In  the  absence  of  debonding,  the  Im  value  decreases  with  increasing 
strain  (Figure  8.6(a)  )  due  to  continuing  process  of  fragmentation.  For  the  sake  of  claritythe 
variation  only  after  4.5%  strain  is  shown  in  the  Figure;  more  than  25  breaks  have  occurred 
at  this  stage  as  can  be  seen  from  Figure  8.6,  causing  a  drop  in  Im  from  10,000  /xm  to  400 
/xm.  The  slope  of  the  curve  in  Figure  8.6(a)  is  very  steep  in  the  initial  regions  where  the 
applied  strain  is  less  than  10%.  The  slope  of  the  curve  drops  after  this  point  as  the  fragment 
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lengths  become  smaller  than  the  load  transfer  length.  Saturation  is  not  reached  at  even  40 
%.  (Strains  above  5%  are  unrealistic  from  an  experimental  point  of  view,  but  the  plot  was 
taken  up  to  50%  to  show  the  trend  in  the  curve).  By  contrast,  saturation  is  seen  at  only  4% 
when  Tdeb  was  taken  at  60  MPa  (Figure  8.6(b)  and  8.6).  As  can  be  seen  from  Figures  8.6 
and  8.6  the  number  of  fiber  breaks  at  any  strain  is  much  lower  for  the  weakly  bonded  case. 
This  is  due  to  the  extensive  debonding  which  limits  the  number  of  breaks  and  also  results 
in  mud:  lower  strain  for  saturation. 

The  evolution  total  debonding  of  the  fiber  as  a  function  of  applied  strain  is  shown 
in  Figure  8.6.  A  comparison  of  saturation  strain  in  Figure  8.6  with  the  corresponding  % 
debonded  length  in  Figure  8.6  indicates  that  more  than  80%  of  the  fiber  length  has  debonded 
at  saturation.  Also,  comparing  the  debonding  with  the  fiber  breaks  (Figure  8.6),  it  is  seen 
that  the  debonding  occurs  most  when  breaks  are  occurring,  and  reduces  once  the  breaks 
stop.  This  is  because  the  frictional  forces  in  the  debonded  region  will  become  significant 
after  extensive  debonding  and  resist  further  debonding  in  the  broken  fragments.  Any  new 
break  will  have  minimum  resistance  from  friction  and  hence  can  debond  more  than  an  already 
debonded  fragment;  but,  as  noted  earlier,  debonding  considerably  reduces  the  probability  of 
further  breaks  in  the  broken  fragments. 


8.5.3  Effect  of  variation  in  interfacial  properties 

Bond  Strength 

A  parametric  study  is  conducted  to  examine  the  effect  of  bond  strength  on  the  mean  fragment 
length  at  saturation.  Figure  8,6  shows  the  variation  of  Im  with  bond  strength.  It  can  be 
seen  that  higher  bond  strengths  lead  to  a  drop  in  the  mean  fragment  length  (increase  in  the 
numb''x  of  breaks)  for  a  given  applied  strain,  as  the  initiation  of  debonding  is  postponed. 
This  k-  also  evident  from  the  increase  in  the  peak  value  in  the  histogram  showing  the  number 
of  breaks  versus  fragment  size  (Figure  8.6).  It  can  be  seen  that  with  increase  in  bond 
strength  the  spread  in  the  distribution  of  fragment  size  significantly  decreases  with  the  peak 
value  moving  to  the  left  (smaller  fragment  size).  The  histogram  is  plotted  from  simulations 
performed  with  more  than  40  different  seeds  to  generate  a  reasonable  volume  of  sampling 
data.  It  was  found  that  the  effect  of  sampling  size  is  much  more  significant  on  the  fragment 
length  distribution  than  on  the  mean  fragment  length. 


Frictioncd  Stress 

Increased  frictional  forces  cause  a  reduction  in  the  debonded  length.  This  results,  in  an 
increase  in  the  number  of  breaks  which  is  clear  from  the  reduction  in  Im  in  Figure  8.6.  A 
comparison  of  the  slopes  in  Figures  8.6  and  8.6  reveals  that  the  sensitivity  of  Im  to  change  in 
the  frictional  properties  is  much  higher  when  compared  to  a  corresponding  drop  in  the  bond 
strength,  especially  at  lower  values  of  Tfj.]  but  as  the  value  of  Tjr  approaches  Tdeb^  further 
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decrease  in  Im  is  negligible. 


8.5.4  Comparison  with  Experimental  Data 

Lacroix  et  al.  [12]  performed  single  fiber  fragmentation  tests  on  carbon-epoxy  composites  and 
measured  tlie  number  of  breaks  as  a  function  of  applied  strain.  Carbon  fibers  with  different 
surface  treatments  were  used  in  the  study.  They  found  that  fiber  surface  treatments  modify 
the  interfacial  conditions  and  affect  the  experimental  results.  Figures  8.6  shows  a  comparison 
of  the  results  of  fiber  fragmentation  simulation  with  experimental  data.  Comparisons  are 
made  for  a  fiber  treatment  of  10%  standard  surface  treatment  (SST).  The  figure  shows 
the  scatter  in  the  experimental  data  for  five  tests.  The  numerical  simulation  in  Fi^re 
8.6  is  the  average  of  a  number  of  runs  with  different  seeds  for  random  number  generation. 
The  material  constants  used  in  the  simulation  are  shown  in  Table  2.  A  good  correlation 
between  the  simulation  and  experiment  is  obtained  for  the  mean  fragment  length  values 
with  Tdeb  =  ITO.OMPa  and  =  SO.OMPa.  It  can  be  seen  that  in  the  experiment  fiber 
breaks  occur  at  a  lower  value  of  applied  strain  compared  to  the  simulation  (curve  A).  Fiber 
statistics  dominate  the  occurrence  of  the  first  few  breaks  and  hence  a  difference  in  the  Weibull 
parameters  possibly  due  to  the  deterioration  of  the  fiber  strength  in  the  composite  may  be 
the  cause  of  this  shift.  To  confirm  this  the  simulation  is  repeated  after  modifying  the  Weibull 
scale  parameter  (curve  B).  As  expected,  a  lower  value  of  the  scale  parameter  (7300  instead 
of  8110)  shifted  the  curve  to  the  left,  closer  to  the  experimental  results. 


8.6  Summary  and  Conclusions 


A  numeric?!  simulation  of  the  single  fiber  fragmentation  test  is  performed,  taking  into  ac¬ 
count,  the  ’'i'lndom  variation  of  the  flaw  strength  of  the  fiber.  The  effect  of  debonding  and 
frictional  stress  on  the  mean  fragment  length  Im  3-t  saturation  is  studied.  Numerical  simula¬ 
tion  shows  that  the  variation  of  cumulative  number  of  breaks  and  debond  progression  with 
increasing  applied  strain  is  significantly  affected  by  the  interfacial  properties.  It  is  found 
that  in  the  case  of  weak  interfaces  saturation  of  fiber  fragments  occurs  at  much  lower  ap¬ 
plied  strain  values  and  is  usually  associated  with  extensive  debonding  (more  than  80%  of  the 
specimen  length  in  some  cases).  It  is  seen  that  an  increase  in  either  Tdeb  and  T/r  decreased 
the  Im  value.  The  results  from  the  numerical  simulation  are  correlated  with  the  experiments 
on  carbon-epoxy  composite  systems  and  a  reasonably  good  match  is  observed. 
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(a) 


Figure  8.1:  Stress  distribution  with  no  debonding,  (a)  axial  stress  in  the  fiber;  (b) interfacial 
shear  stress 


(a) 


Figure  8.2:  Stress  distribution  with  partial  debonding,  (a)  axial  stress  in  the  fiber;  (b) 
interfacial  shear  stress 
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Figure  8.3:  Flow  chart  for  the  simulation  of  the  fiber  fragmentation  test 
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Figure  8.4;  Axial  stress  in  fiber  with  strong  bond  (a)  at  intermediate  applied  strain,  (b)  high 
strain 
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Figure  8.5;  Axial  stress  profile  with  debonding  (a)  at  initial  loading  and  intermediate  stages, 
(b)  at  saturation 
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Mean  fragment  length  (microns)  Mean  fragment  length  (microns) 


Figure  8.7:  (a);Variation  of  average  fragment  length  with  applied  strain  for  a  strong  bond. 
(b):Variation  of  average  fragment  length  with  applied  strain  for  a  weak  bond. 
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Figure  8.11:  Fragment  length  distribution  with  increasing 
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Figure  8.15:  exp  comp  Im 
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Chapter  9 


Summary  and  Conclusion 

9.1  Introduction 


The  elevated  temperature  specific  stiffness  and  strength  properties  of  Metal  Matrix  Com¬ 
posites  (MMCs)  have  made  them  candidate  materials  to  meet  the  stringent  requirements  of 
the  future  U.S.  Air  Force  aircrafts.  The  successful  application  of  such  composites  not  only 
relies  on  the  use  of  right  fiber  and  matrix  materials,  but  strongly  depends  on  the  interface 
between  the  fiber  and  matrix.  The  processing  conditions  prevailing  during  the  consolidation 
of  these  composites  play  a  critical  role  in  the  evolution  of  the  fiber-matrix  interface. 

A  fundamental  understanding  of  the  relationship  between  the  processing  conditions, 
chemistry  of  the  reaction  zone  at  the  interface,  and  the  mechanical  properties  of  the  interface 
is  important  to  the  successful  use  of  MMCs. 


9.2  Overall  goal  of  the  research  work 

The  primary  objectives  of  the  research  work  are 

•  To  establish  a  link  between  the  interface  chemistry  and  its  mechanical  properties  of 
Titanium  Matrix  Composites. 

•  To  study  the  effect  of  interfacial  properties  and  inelastic  behavior  of  the  matrix  on  the 
evolution  of  damage  in  Metal  Matrix  Composites. 

The  research  effort  involves  composite  consolidation,  controlled  exposure  studies,  chemical, 
metallurgical  and  mechanical  characterization  and  numerical  modeling  of  selected  Titanium 
Metal  Matrix  Composites. 
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9.3  Basic  Research  Issues 


The  role  of  fiber-matrix  interface  is  crucial  to  the  performance  of  the  composites.  Interfacial 
debonding  is  one  of  the  common  damage  mechanism  in  most  fiber-reinforced  MMCs.  The 
fatigue,  fracture,  transverse  (90°  to  the  fiber  direction)  stress-strain  behavior,  are  found  to 
be  strongly  influenced  by  the  interfacial  bonding  conditions  [1].  The  bond  strength  of  an 
interface  is  the  net  effect  of  the  contributions  from  chemical  and  mechanical  components 
[21.  The  former  is  due  to  the  reaction  at  the  interface  during  the  composite  processing  and 
depends  on  the  fiber  and  matrix  materials  and  the  processing  conditions.  The  mechanical 
component  is  due  to  the  sliding  resistance  at  the  interface  and  depends  on  the  presence  of 
residual  stresses  and  interfacial  surface  irregularities.  It  is  necessary  to  understand  the  role 
of  this  mechanical  component  while  characterizing  the  interfacial  strength  and  toughness. 

The  interfacial  chemistry  evolves  with  time  exposure  to  elevated  temperature  and/or 
environment,  and  influences  both  the  size  and  phases  within  the  reaction  zone  of  the  interface. 
Titanium  alloys  and  intermetallics  which  are  the  matrix  materials  in  this  work  are  very  highly 
reactive  at  elevated  temperatures  and  produce  nonstoichometric  carbides  and  silicides  within 
the  interphases.  Chemical  compositional  and  phase  evolution  through  an  understanding  of 
the  diffusional  process  or  otherwise  is  necessary. 

Residual  stresses  develop  when  the  MMCs  and  IMCs  are  cooled  after  processing  from 
a  reference  temperature  at  which  they  are  stress  free,  to  room  temperature.  This  is  due 
mainly  to  the  difference  in  the  coefficients  of  thermal  expansion  (CTE)  of  fiber  and  matrix 
material.  These  stresses  reach  considerable  levels  because,  in  typical  MMCs  and  IMCs 
considered  in  this  study,  the  mismatch  in  the  CTE  is  quite  significant  (amatrix  >  ^OLjiher) 
and  the  temperature  differential  involved  in  the  cooling  process  is  high  (about  800°C  to 
900°C).  Residual  stresses  are  affected  by  the  evolution  of  interfacial  chemistry,  and  by  the 
inelastic  response  of  matrix  material. 

The  mechanical  response  of  the  interfaces  is  usually  measured  using  a  thin  snce  push  out 
test,  and  transverse  tensile  testing.  However,  extracting  quantitative  interfaciaJ  properties 
from  the  push-out  test  is  not  trivial.  A  comprehensive  process  model  (including  composite 
consolidation  and  specimen  preparation)  to  understand  the  interfaciaJ  failure  process  dur¬ 
ing  push-out  tests  is  needed  for  better  interpretation  of  the  test  results.  In  the  numerical 
simulation,  either  a  stress  based  failure  criterion  or  an  energy  based  criterion  can  be  used  to 
model  the  failure  process. 

In  summary,  basic  research  issues  in  this  effort  focus  on  the  understanding  the  chemical 
evolution  (diffusional)  process  at  the  interface,  accurately  computing  the  residual  stresses 
in  the  interface  zone  (accounting  for  matrix  inelasticity,  and  load  shedding  during  cpack 
propagation),  developing  a  mechanistic  model  of  an  interface  based  on  stress  or  energy  failure 
criterion,  and  studying  the  effect  of  interface  properties  on  the  performance  of  the  composites. 
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9.4 


Technical  Approach 


♦ 


The  technical  approach  involves  both  experimental  and  computational  components.  For 
our  purpose,  two  model  composite  systems  were  selected,  Ti-6Al-4V/SSC-6  and  Timetal 
21S/SCS-6.  While  the  former  was  consolidated  within  our  laboratory,  the  latter  was  ob¬ 
tained  from  other  sources.  The  effect  of  elevated  temperature  was  studied  by  heat  treating 
specimens  at  selected  temperatures  (450°  C,  700°  C  and  927°  C)  in  vacuum  for  varying  peri¬ 
ods  of  time.  The  time  periods  used  are  25,  70  and  120  hours.  Scanning  Electron  Microscope 
(SEM)  is  used  to  measure  the  thickness  of  the  interface.  Interfacial  chemistry  and  elemental 
composition  are  monitored  using  Energy  Dispersive  X-ray  Spectroscopy  (EDX). 

Push-out  tests  to  measure  interfacial  mechanical  properties  are  conducted  at  NASA 
Lewis  Research  Center  (thanks  to  Dr.  J.  Eldridge),  on  the  heat  treated  specimens  and  the 
results  are  compared  with  those  of  as-processed  ones.  Modeling  push-out  test  involves  three 
major  stages.  Step  1  involves  modeling  the  composite  consolidation  at  high  temperatures 
and  then  cooling  down  to  room  temperature  which  induces  residual  stresses  in  the  composite. 
Step  2  is  the  specimen  preparation  stage  consisting  of  slicing  a  thin  push-out  specimen  from 
the  bulk  composite  which  results  in  stress  redistribution  at  the  two  ends  of  a  fiber.  Step  3  is 
modeling  the  loading  stage  where  compressive  load  is  applied  to  the  fiber  until  it  is  pushed 
out  of  the  specimen.  A  detailed  description  of  the  above  sequence  of  steps  is  provided 
elsewhere  [3].  The  interface  is  modeled  as  a  contact  surface  and  the  bonding  at  the  interface 
is  simulated  using  high  stiffness  springs.  The  interface  failure  process  is  modeled  using  a 
strain  energy  based  failure  criteria  given  by: 


Gt  >  Gi, 

where  Gt  is  the  strain  energy  release  rate  for  the  interface  crack  and  Gic  is  the  critical  value  of 
the  strain  energy  release  rate.  The  interface  crack  propagates  further  when  its  strain  energy 
release  rate  exceeds  the  critical  value.  A  numerical  method  based  on  the  equivalent  domain 
integral  technique(EDI)[8,  9,  4]  has  been  developed  by  the  authors[5,  7]  to  evaluate  the  strain 
energy  release  rate  of  propagating  interface  cracks.  The  EDI  method  was  implemented  in  the 
commercial  finite  element  analysis  software  MARC [6]  through  user  subroutines.  It  is  then 
used  in  the  modeling  of  interface  failure  during  thin-slice  fiber  push-out  test.  Figure  2  shows 
the  application  of  numerical  simulation  to  evaluate  fracture  toughness  from  the  experimental 
data. 


9.5  Significant  Results 


The  evolution  of  interface  in  Titanium  Matrix  Composites  (e.g.  SCS-6/Timetal  21s)  com¬ 
posite  system  is  studied  by  means  of  heat  treatment  experiments,  with  the  objective  of 
establishing  a  link  between  the  evolution  of  chemistry  and  geometry  of  the  reaction  zone 
and  the  interfacial  properties. 
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•  Microstructural  studies  reveal  that  reaction  growth  is  due  to  transformation  of  both 
coating  and  the  matrix.  The  growth  of  the  reaction  zone  into  the  matrix  is  much  more 
than  the  growth  towards  the  fiber,  as  carbon  is  able  to  diffuse  much  farther  into  the 
matrix  than  the  diffusion  of  Titanium  into  the  fiber. 

•  Reaction  zone  sizes  increase  rapidly  with  time  at  927°  C.  At  lower  temperatures  such 
as  700°  C  and  450°  C  no  significant  increase  in  reaction  zone  size  is  observed  for  this 
system. 

•  Push-out  tests  conducted  on  heat  treated  specimens  reveal  that  the  effect  of  thermal 
exposure  on  interfacial  properties  is  predominantly  due  to  exposure  temperature.  Time 
of  exposure  is  found  to  have  only  a  secondary  effect. 

•  Load  displacement  curves  for  as-processed  SCS-6/Timetal  21s  specimens  and  speci¬ 
mens  heat  treated  at  700°  C  and  450°  C  have  sharp  load  drops  at  the  point  of  complete 
debonding  marked  by  distinct  acoustic  emission  signals.  The  results  for  specimens  heat 
treated  at  927°  C  show  that  debonding  occurs  as  a  slow  process  with  no  observable 
sharp  load  drops  and  acoustic  emission  signals. 

•  Fracture  toughness  values  remain  more  or  less  constant  for  normal  use  temperatures 
in  the  range  of  450°  to  700°  C.  The  values  however,  increase  when  the  exposure  tem¬ 
perature  is  raised  to  927°  C. 

•  Fracture  surfaces  of  the  specimens  exposed  at  927°  C  show  that  debonding  took  place 
in  the  reaction  zone  whereas  in  as-processed  specimens  and  those  heat  treated  at  the 
lower  temperature  ranges  debonding  initiated  in  the  coating.  This  is  because  the  weak 
region  in  the  coating  is  replaced  by  the  stronger  interfacial  region. 


9.6  Relevance  to  Air  Force  and  Future  Direction 

From  the  above  study,  it  is  clear  that  interfaces  play  a  crucial  role  not  only  in  MMCs  but 
in  high  temperature  heterogeneous  systems,  e.g.,  thermal  barrier  coatings.  Clearly  a  deeper 
understanding  of  the  interfacial  fracture  processes  is  needed  using  an  innovative  experimental 
measurements  and  advanced  compuational  schemes  at  micron  levels. 
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